AD'A154  292  MUTURL  COUPLING  RNRLVSIS  FOP  CONFORHRL  HICROSTRIP 

RNTENNRS(U)  OHIO  STATE  UNIV  COLUHBUS  ELECTROSCIENCE  LAB 
B  U  KURN  ET  RL.  DEC  84  712692-4  N0eei4-78-C-e049 
UNCLRSSIFIED  F/G  9/5  NL 


AD- A 154  292 


Kl  PRODUCrn  AT  GOVERNMENT  EXPENSE 


The  Ohio  State  University 


MUTUAL  COUPLING  ANALYSIS  FOR  CONFORMAL 
MICROSTRIP  ANTENNAS 


B.W.  Kwan  and  E.H.  Newman 


The  Ohio  Slot*  Univenity 

ElectraScieiice  Labaratwy 

Oaportmant  of  Eloctricol  Enginooring 
ColumbM,  0(«(0  43212 


Technical  Report  712692-4 
Contract  No.  N00014-78-C-0049 


December  1984 


Department  of  the  Navy 
Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  Virginia  22217 


IDIOT 
Api 


DISnUBUTlON  STATEMENT  A 
Approved  k«  pubUo  rolearoi  | 


DTIC 

ELECTE| 
MAY30  686 


Diatfibatton  UnHwtIed 


85  5  09  058 


I.,  J  ...y  a.  J  J  J  J  Iji.lijilji  piiwii,  IP  IP  J  ,1  JI.JIIHHIMHII 

lu  ('fionocro  Ai  (;ovrRNMENf  expense 


NOTICES 


When  Government  drawings,  specifications,  or  other  data  are 
used  for  any  purpose  other  than  in  connection  with  a  definitely 
related  Government  procurement  operation,  the  United  States 
Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever,  and  the  fact  that  the  Government  may  have  formulated, 
furnished,  or  in  any  way  supplied  the  said  drawings,  specifications, 
or  other  data,  is  not  to  be  regarded  by  implication  or  otherwise  as 
in  any  manner  licensing  the  holder  or  any  other  person  or  corporation, 
or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell 
any  patented  invention  that  may  in  any  way  be  related  thereto. 


S0772-I0I 


REPORT  DOCUMENTATION  1 »  "eport  no. 
PAGE  I 

4.  THI«  and  SuMill* 


2.y|  3.  Ra<')p* 


HI*  •Ad  SuMitIo  I  s.  Report  Dote 

MUTUAL  COUPLING  ANALYSIS  FOR  CONFORMAL  MICROSTRIP  ANTENNAS  I  December  1984 


7.  Author<«)  .  ^  . 

B.W.  Kwan  and  E.H.  Newman 

9.  Performing  Orgenization  Nome  and  Address 


The  Ohio  State  University  ElectroScience  Laboratory 
Department  of  Electrical  Engineering 
1320  Kinnear  Road 
Columbus,  Ohio  43221 


t.  Performing  Organization  Rept.  No. 

712692-4 

10.  Proioct/Taak/Worli  Unit  No. 


11.  Contract(C)  or  Grant(G)  No 

<C'  N00014-78-C-0049 


12.  Sponsorir^  Organization  Name  and  Address 

Department  of  the  Navy 
Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  Virginia  22217 

15.  Supplementary  Notes 


13.  Type  of  Report  A  Period  Covered 

Technical 


IS.  Abstract  (Limit:  2(X)  words) 


Integral  equations  based  upon  appropriate  Green's  functions  are  employed  to 
rigorously  treat  the  mutual  coupling  of  microstrip  patches  on  1)  an  infinite  grounded 
dielectric  slab,  and  2)  a  dielectric  coated  circular  cylinder.  A  moment  method  solution 
for  the  first  problem  using  the  dielectric  slab  Green's  function  is  presented.  Calculated 
values  of  the  mutual  impedance  are  found  to  be  in  good  agreement  with  measurements.  The 
second  problem  is  restricted  to  the  calculation  of  the  elements  of  the  impedance  and 
voltage  matrices  which  are  needed  to  carry  out  the  solution  by  the  moment  method.  These 
elements  are  derived  using  the  Green's  dyadic  for  the  coated  cylinder.  It  is  found  that 
the  self  impedance  of  a  dipole  mode  at  the  surface  of  a  coated  cylinder  converges  to  that 
of  the  infinite  slab  case  as  the  radius  of  the  cylinder  increases.  Two  examples  of  the 
mutual  coupling  between  two  dipole  modes  are  considered  and  numerical  results  are 
presented  for  cylinders  of  different  size.  In  addition,  alternative  representations 
suitable  for  asymptotic  evaluation  are  derived  from  the  eigenfunction  solution  obtained  by 
the  Green's  dyadic.  A  two-dimensional  example  is  presented  to  illustrate  the  usefulness 
of  asymptotic  solutions  when  the  dielectric  coated  cylinder  is  large  in  terms  of  a 
wavelength.  ^ 


17.  OoeumaM  Analy*!*  ■.  OMcrIptor* 

Microstrip  Antennas^ 
Mutual  Coupling 
Conformal  Arrays, 

b.  Idantlfl.rv/Op.n.Cnd.d  T.rmt 


Green's  Function, 
Cyl i nder 
Moment  Method  , 


Eigenfunction  Expansion 


c.  COSATl  Fi«ld/Group 
!■»  AvPilPbIbty  Slplpmpnt 


DISlWBUnON  STAreMENT-ft 


Approved  iw  public  releoeei 
Diatiibution  Uolimited 


It.  SMurily  Cl.ti  (Thi*  R.per1) 


leieoeei  UNCLASSIFIED.. _ 

ated  TO.  Sacurtty  Claat  (Thu  Ra(a) 

I 


Saa  InairwetJona  on  Ravaraa 


21.  No  of  Pacaa 


(Saa  ANSI-ZM.im 


OrriONAL  FORM  272  (4-77) 
(Formarly  NTIS-SS) 
Oaparlmant  of  Commana 


TABLE  OF  CONTENTS 


Page 

LIST  OF  TABLES  vi 

LIST  OF  FIGURES  vi i 

Chapter 

I.  INTRODUCTION  AND  PRELIMINARY  1 

A.  BRIEF  REVIEW  1 

R.  STATEMENT  OF  PROBLEM 

C.  OVERVIEW  OF  APPROACH  6 

1.  Integral  Equation  Formulation  6 

2,  Moment  Method  Solution  6 

3,  Feed  Modeling  H 

4.  Port  Impedances  10 

D.  ORGANIZATION  13 

II.  MICROSTRIP  PATCHES  ON  A  GROUNDED  PLANAR 

DIELECTRIC  SLAB  15 

A.  INTRODUCTION  15 

B.  THEORY  AND  GENERAL  SOLUTION  16 

C.  MUTUAL  COUPLING  ANALYSIS 

0.  NUMERICAL  EXAMPLES  36 

E.  FAR  ZONE  RADIATION  FIELD  46 


L. 

t  . 


Chapter 

Page 

III. 

DYADIC  GREEN'S  FUNCTIONS 

50 

A. 

INTRODUCTION 

50 

B. 

SOME  BASIC  RELATIONS  GOVERNING  E,  H.  Gg  AND  ^ 

53 

F. 

C. 

EIGENFUNCTION  EXPANSION  OF  THE  FREE  SPACE 

GREEN'S  DYADIC 

58 

0. 

GREEN'S  DYADIC  FOR  A  DIELECTRIC  COATED 

CONDUCTING  CYLINDER 

79 

E. 

EXPANDING  THE  GREEN'S  DYAOICS 

91 

1^ 

IV. 

PATCH  DIPOLES  ON  A  DIELECTRIC  COATED  CYLINDER 

98 

A. 

INTRODUCTION 

98 

B. 

A  SPECIALIZED  METHOD  OF  SOLUTION 

101 

r 

C. 

MUTUAL  COUPLING  BETWEEN  EXPANSION  DIPOLE  MOOES 

113 

L 

D. 

NUMERICAL  EXAMPLES 

m 

t-f 

E. 

FAR  FIELD  CALCULATION 

144 

V. 

ASYMPTOTIC  SOLUTIONS 

153 

A. 

INTRODUCTION 

153 

s 

8. 

ALTERNATIVE  REPRESENTATIONS  OF  THE  GREEN'S  DYADIC 

154 

< 

C. 

ASYMPTOTIC  EVALUATION  OF  Z12 

163 

D. 

A  NUMERICAL  EXAMPLE 

18? 

VI. 

SUMMARY 

185 

i  V 

•V 
•  • 

ft  ^ 

I 

\ ' 

APPENDICES 


A.  DETERMINATION  OF  SPECTRAL  FUNCTIONS 

B.  MUTUAL  IMPEDANCE  BETWEEN  TWO  EXPANSION 
(DIPOLE)  MOOES 

C.  EVALUATION  OF  AN  OSCILLATORY  INTEGRAL  BY  THE 
STATIONARY  PHASE  METHOD 

D.  ORTHOGONALITY  PROPERTIES  OF  CYLINDRICAL  VECTOR 
EIGENMODES 

E.  PROOF  OF  EQUATION  (3.73) 

F.  DETERMINATION  OF  THE  EXPANSION  COEFFICIENTS 

G.  ON  THE  FUNCTION 

H.  ON  THE  RATIOS  Pp/Sp  AND  Qp/Tp 

I.  CALCULATION  OF  CYLINDER  FUNCTIONS:  j'  H  h' 

n’  n*  n»  n 

J.  ASYMPTOTIC  FORM  OF  Q( v) 

K.  THIN  SUBSTRATE  APPROXIMATION  OF  Zv(h) 

L.  ROOTS  OF  P(a)  =  ywi(a)  -  wi(a)  =  0 


REFERENCES 


MAY  3  0  1985 


Aoc!.::.'.  von  For 

i- 

1  '  - 


IV.  n 

.need  Q 


Hy- 

b,!  •:  ;■<  ’'.V.vlcu/ 


Av  U.'ibillty  Codes 
Avfjill  and/or 
i  Special 


LIST  OF  TABLES 


Page 

SELF  IMPEDANCE  OF  A  SINGLE  EXPANSION  MODE  138 

MUTUAL  IMPEDANCE  BETWEEN  TWO  IDENTICAL  EXPANSIDN 

MODES  IN  THE  E-PLANE  140 

MUTUAL  IMPEDANCE  BETWEEN  TWO  IDENTICAL  EXPANSION 

MODES  IN  THE  H-PLANE  141 


LIST  OF  FIGURES 


Figure  Page 

1.1.  Microstrip  antenna  configuration.  3 

1.2.  Microstrip  antenna  feed  modeling.  9 

1.3.  Voltages  and  currents  of  an  M-port.  11 

2.1.  Geometry  of  a  microstrip  patch  printed  on  a  grounded 

dielectric  slab.  17 

2.2.  Proper  contours  of  integration  and  branch  cuts  in  the 

complex  k  plane.  26 

2.3.  Geometry  of  patches  and  expansion  dipole  modes  for  the 

evaluation  of  and  v^fa)  32 

mn  m'^'* 

2.4.  Current  distribution  on  an  isolated  microstrip  patch  at 

resonance.  38 

2.5.  E-plane  coupling  for  identical  microstrips.  40 

2.6.  H-plane  coupling  for  identical  microstrips.  43 

3.1.  Region  V,  bounded  by  a  perfect  conducting  surface  S  and 
a  surface  S„  at  infinity,  contains  an  electric  source  7. 

V  can  be  inhomogeneous.  56 

3.2.  A  general  orthogonal  curvilinear  system  with  coordinates 

^1’  ^2’  *^3  their  corresponding  unit  vectors  Uj,  u^, 
and  metric  parameters  hj,  h2,  h3.  59 

3.3.  Free  space  V^  is  partitioned  into  regions  V  and  V  with 
closed  surface  S  as  interface  which  contains  the  point 

•  III 

source  J  at  r'  =  (ui,u2,U3). 

vii 


3.4.  Configuration  for  determining  Dg  ,  Pg  and  Te 

0+0  0  ^  0^ 

-  region  V  is  bounded  by  surfaces  S  at  p=p  and 
S“  at  p=p“,  and  source  T  is  located  on  surface  S  at 

p=p'  ( p'''>p‘  >p“) . 

3.5.  Geometry  of  a  dielectric  coated  cylinder. 

4.1.  Two  z-polarized  dipole  modes  on  the  surface  of  a 
dielectric  coated  cylinder. 

4.2.  Proper  contours  of  integration  and  branch  cuts  in  the 
complex  h-plane. 

4.3.  E-plane  coupling  for  identical  expansion  modes  on  coated 
cylinders  of  radii  a  =  0.25Xo,  o.BXq,  I.Xq  and  10. Xr 

( Xq  =  0.474  m.) . 

4.4.  H-plane  coupling  for  identical  expansion  modes  on  coated 
cylinders  of  radii  a  =  0.25Xq,  O.BXq,  I.Xq  and  3.Xo 

( Xq  =  0.474  m.) . 

4.5.  Coupling  for  identical  expansion  modes  on  a  coated 
cylinder  of  radius  a  »  0.25Xo  (Xq  =  0.474  m.). 

4.6.  Coupling  for  identical  expansion  modes  on  a  coated 
cylinder  of  radius  a  =  O.BXq  (Xq  =  0.474  m.). 

4.7.  Coupling  for  identical  expansion  modes  on  a  coated 
cylinder  of  radius  a  =  I.Xq  (Xq  =  0.474  m.). 

4.8.  Contours  of  integration  in  the  h-plane  and  the  w-plane. 

4.9.  The  steepest  descent  path. 

5.1.  (a)  Contour  Cy  encloses  the  upper  half  v-plane  which 

contains  the  poles  of  gg  or  ggZ  vj,  V2, . . . , Vp. . . . 

(b)  The  deformed  contour  Cy'  also  encloses  all  the  poles 
of  ge  or  gg. 


Cross-section  of  the  dielectric  coated  cylinder  with  two 
infinite  conducting  strips. 


168 


Comparison  between  integral  representation,  residue  series 
and  eigenfunction  expansion  of  1^3 

Two  expansion  dipole  modes  on  a  grounded  dielectric  slab.  198 

Path  of  integration  r  the  h-plane  and  the  analytic 
properties  of  ^k^-h^  when  k''=0  (for  e“3^  time  dependence).  244 
Path  of  integration  r  in  the  h-plane  and  the  analytic 
properties  of  ^k^-h^  when  k"=0  (for  0^“^^  time  dependence). 


CHAPTER  I 


INTRODUCTION  AND  PRELIMINARY 


A.  BRIEF  REVIEW 

A  microstrip  device  basically  consists  of  a  sandwich  of  two 
parallel  conducting  layers  separted  by  a  single  thin  dielectric 
substrate.  The  lower  conductor  functions  as  a  ground  plane,  and  the 
upper  conductor  may  be  a  simple  resonant  patch  of  regular  shape,  a 
resonant  dipole,  or  a  monol ithically  printed  array  of  patches  or  dipoles 
and  the  associated  feed  network.  The  concept  of  microstrip  antennas  was 
first  proposed  in  this  country  a  by  Deschamps  [1]  and  in  France  by 
Gutton  and  Baissinot  [Zl.  However,  the  main  interest  in  the  microstrip 
antenna  concept  did  not  emerge  unitl  the  early  1970's.  It  was  probably 
the  need  for  conformal  missile  and  spacecraft  antennas  that  provided 
this  impetus  together  with  the  advance  of  printed  circuit  technology 
based  on  the  good  selection  of  microstrip  substrates  that  were  becoming 
available.  The  first  practical  antennas  were  developed  in  the  early 
1970's  by  Howell  [31  and  Munson  [4].  Since  then,  extensive  research 


and  development  of  microstrip  antennas  have  led  to  diversified 
applications  and  to  the  establishment  of  the  topic  as  a  separate  entity 
within  the  broad  field  of  microwave  antennas.  The  state  of  the  art,  in 
both  theory  and  experiment,  is  summarized  in  the  books  by  Bahl  and 
Bhartia  [51,  and  James,  et  al .  [61.  In  particular,  a  recent  special 
issue  of  IEEE  Transactions  on  Antennas  and  Propagation  [71  was  devoted 
to  microstrip  antennas  and  arrays.  In  this  special  issue.  Carver  and 
Mink  [81  presented  a  comprehensive  survey  of  the  state  of  microstrip 
antenna  element  technology;  while  Mailloux,  et  al .  [9]  discussed 
microstrip  array  design  techniques. 

In  its  simplest  configuration  a  microstrip  antenna  consists  of  a 
thin  metallic  radiating  patch  bonded  to  a  thin  grounded  dielectric 
substrate  (Figure  1.1).  The  patch  conductor  typically  has  some  regular 
shape,  for  example,  rectangular,  circular  or  elliptical.  The  feed  is 
often  a  coaxial  probe  or  a  microstrip  transmission  line.  Microstrip 
antennas  exhibit  all  the  properties  inherent  to  microstrip  devices:  a) 
they  are  light  weight,  small  size  and  low  profile  planar  configurations 
which  can  be  made  conformal;  b)  they  are  inexpensive  to  build  and 
ideally  suited  for  large  scale  production  by  printed  circuit  techniques; 
c)  they  are  compatible  with  modular  designs  (solid  state  devices  such  as 
oscillators,  amplifiers,  phase  shifters,  etc.,  can  be  added  directly  to 
the  antenna  substrate  board);  d)  their  feed  lines  and  matching  networks 
can  be  fabricated  simultaneously  with  the  antenna  structure  so  that 
discontinuities  due  to  connectors  can  be  eliminated.  All  these 
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RADIATING  PATCH 


Figure  1.1.  Microstrip  antenna  configuration. 

advantages  compensate,  at  least  In  part,  a  number  of  drawbacks:  a) 
simple  microstrip  antennas  have  narrow  bandwidths;  b)  their  gain  is 
low;  c)  they  have  a  small  power  handling  capability;  d)  their  dielectri 
losses  reduce  the  radiation  efficiency;  e)  unwanted  surface  waves  may 
cause  spurious  radiation  at  the  edges  of  the  microstrip  patch. 

Most  of  the  work  in  printed  antenna  theory  and  technology  during 
the  last  decade  was  for  antennas  operating  in  the  UHF  to  microwave 
frequency  bands  (300  MHz  to  10  GHz).  There  are  two  somewhat  successful 
methods  in  calculating  input  impedance  and  radiation.  One  is  the 
transmission  line  model  proposed  by  Derneryd  [10],  and  the  other  is 
the  cavity  model  used  by  Lo,  et  al ,  [11].  However,  both  methods  fail 


to  rigorously  account  for  the  surface  waves  on  the  antenna  substrate  and 
for  the  mutual  coupling  between  antenna  elements.  More  recently,  a 
moment  method  solution  to  the  microstrip  strip  antenna  problem  was 
proposed  by  Newman  and  Tulyathan  [12],  In  [12],  an  integral  equation 
was  formulated  with  the  aid  of  image  theory  which  accounted  for  the 
presence  of  the  ground  plane;  the  field  within  the  dielectric  slab  was 
treated  by  an  equivalent  current.  This  method  gives  good  results  for 
the  input  impedance,  but  requires  an  extremely  accurate  evaluation  of 
the  elements  of  the  impedance  matrix. 

The  above  studies  all  deal  with  planar  substrates.  From  a 
practical  viewpoint  it  is  also  important  to  consider  microstrip  antennas 
and  arrays  on  curved  surfaces,  especially  on  portions  of  cylinders, 
cones  or  spheres.  Microstrip  dipoles  have  indeed  been  considered  for 
cylindrical  substrates  and  some  preliminary  results  have  been  published 
by  Alexopoulos,  et  al .  [13],  who  utilized  the  dyadic  Green's  function  to 
give  the  electric  field  produced  by  an  electric  dipole  tangent  to  the 
outer  surface  of  the  coating  layer.  Their  results  were  restricted  to 
the  far  field  and  to  the  surface  field.  At  about  the  same  time,  Fonseca 
and  Attilio  [30]  reported  results  on  the  radiation  patterns  of 
microstrip  wraparound  antennas  using  a  theory  based  upon  the  dyadic 
Green's  function  for  a  dielectric  coated  cylinder.  In  both  papers, 
the  Green's  dyadic  was  constructed  using  the  principle  of  scattering 
superposition.  The  free  space  Green's  dyadic  employed  was  quoted  from 
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Tai  [19]  who  used  a  different  method  from  the  one  adopted  in  this  work 
for  constructing  the  free  space  Green's  dyadic.  However,  the  above 
papers  did  not  specifically  address  the  mutual  coupling  problem. 

B.  STATEMENT  OF  PROBLEM 

It  is  likely  that  there  will  be  increasing  interest  in  millimeter 
wave  systems  and  applications,  such  as  imaging  array  antennas  and 
ai rcraft-to-satel 1 ite  communications.  Also  of  interest  is  the 
development  of  complete  monolithic  systems  which  combine  antenna 
elements  or  arrays  on  the  same  substrate  as  the  integrated  RF/IF 
front-end  detector  and  amplifier  circuits.  In  these  applications, 
substrates  are  often  thicker  and  have  higher  dielectric  constants  than 
at  lower  frequencies.  Consequently  the  electi'ical  performance  of  these 
antennas  will  be  severely  degraded  due  to  surface  waves  or  mutual 
coupling.  It  is  then  evident  that  the  analysis  of  mutual  coupling 
between  microstrips  is  important  in  the  design  of  antenna  arrays, 
especially  if  tight  pattern  control  or  low  sidelobes  are  required. 

This  work  considers  the  mutual  coupling  of  two  types  of  conformal 
microstrip  antennas: 

1.  microstrip  patches  on  a  dielectric  slab 
(dielectric  slab  problem); 

microstrip  patches  on  a  dielectric  coated  cylinder 
(dielectric  coated  cylinder  problem). 


C.  OVERVIEW  OF  APPROACH 


1.  Integral  Equation  Formulation 

The  present  approach  is  based  on  the  Rreen's  function  method 
which  yields  the  total  electric  field  E  produced  by  the  electric  surface 
currents  .1$  on  the  microstrip  antenna  elements  (patches)  as 


E  =  i’  +  /  G 

S 


Jg  ds 


(1.1) 


where  E^  is  the  incident  field  excited  by  an  impressed  curent  source  3-; 
G  is  an  appropriate  dyadic  Green's  function  which  accounts  for  the 
grounded  substrate;  and  S  is  the  surface  of  the  microstrip  patches. 

Introducing  the  boundary  condition  on  the  microstrip  patches 


n  X  E  =  0  on  S 


(1.?) 


yields  an  electric  field  integral  equation  for  here  n  is  the  unit 
vector  normal  to  the  surface  of  the  microstrip  patches. 

Moment  Method  Solution 

The  integral  Equation  (1.2)  is  solved  using  the  fialerkin  form  of 
the  moment  method  where  both  the  basis  and  testing  functions  are  taken 
as  a  surface  patch  dipole  mode.  Thus  the  unknown  current  is  expanded 
in  a  set  of  N  basis  functions  or  modes 


6 


where  J,,  is  the  n^h  mode  and  In  is  its  unknown  amplitude.  Use  of  the 
same  set  of  functions  as  testing  function  leads  to  a  system  of  linear 


algebraic  equations  to  be  solved  for  the  unknown  Ip 


rzl[I]  =  [V] 


(1.4) 


where  [Z]  is  known  as  the  impedance  matrix  with  elements 


=  -  /  E  *0  ds 
mn  i  m  n 

^n 


(1.5) 


[Vl  is  the  voltage  vector  whose  elements  are  given  by 


m 


=  /  E, 


V. 


m 


dv 


(1.6) 


and  Ip  are  the  elements  of  the  current  vector  [I],  In  Equations  (1.5) 
and  (1.6),  Em  is  the  electric  field  due  to  the  mode  current  Jm  in  the 
presence  of  the  dielectric  substrate  and  ground  plane;  Ji  is  the 
impressed  (source)  current;  and  Sp,  Vi  denote  the  surface,  volume  where 
currents  Jp,  Ji  exist,  respectively.  The  expressions  in  (1.5)  and  (1.6) 
were  obtained  using  the  reciprocity  principle. 


3.  Feed  Modeling 


Referring  to  Figure  1.2(a),  a  coaxial  feed  with  terminal  current  li 
can  he  modeled  by 

Ji  =  zIi6(x-XQ)6(y-y^)  (I.7) 

A 

to  represent  a  z-directed  current  source  at  the  probe  position  (x  y  ). 
To  account  for  the  probe  self  inductance,  jXp  (for  0"^“^  time  dependence) 
can  be  added  to  the  input  impedance,  where  [8] 

h  — 

V  "  /i;  *  (1-8) 

with  Z  a  /n  /e  ,  k  =  u)/u  e  ,  u  and  e  are  the  constitutive 
0  00  0  00  0  0 

parameters  of  free  space  (air);  and  t  is  the  thickness  of  the  dielectric 
layer  with  dielectric  constant 

In  Figure  1.2(b),  a  microstrip  line  feed  with  terminal  current 
also  can  be  modeled  by  using  Equation  (1.7)  to  represent  an  equivalent 

A 

z-directed  current  source  at  the  point  (’^Qty^)  where  the  feed  line  joins 
the  microstrip  patch.  Such  a  model  furnishes  good  results  for  the 
narrow  feed  lines  that  are  in  common  use.  Alexopoulos  and  Rana  [141 
pointed  out  that  the  voltage  term  (1.6)  should  be  modified  by  the  factor 
•'W^/t  to  account  for  edge  effects  of  the  microstrip  line  of  width  W, 
where  is  the  effective  width  given  by  [81 
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_  -  1/2 
lot 

S  =  2  2  ^  ^  W  •  (1*10) 

The  feed  model  (Figure  1.2(c))  can  be  thought  of  as  a  port  with  terminal 
current  Ij  and  voltage  Vi  as  shown  in  Figure  1.2(d).  Vi  is  the  voltage 
due  to  surface  current  on  the  patch  excited  by  the  impressed  Ji , 

4.  Port  Impedances 

An  M-element  microstrip  array  can  be  modeled  as  an  M-port  whose 
currents  and  voltages  (defined  in  Figure  1.3)  are  related  in  matrix  form 
as 

[ZP]  [iP]  =  [yPj  .  (1.11) 

One  should  differentiate  between  these  port  quantities  (indicated  by  the 
superscript  p)  and  those  associated  with  the  moment  method  solution 
defined  in  Equations  (1.4)  -(1.6).  The  moment  method  quantities  depend 
on  the  choice  of  the  basis  functions.  The  port  quantities  are  the  ones 
of  interest  for  determining  input  impedance  and  mutual  coupling.  The 
relation  between  these  quantities  is  now  described. 


e  +  1  e  -  1 

r  r 
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p 

X* 


[z^] 


Figure  1.3.  Voltages  and  currents  of  an  M-port. 


The  input  impedance  at  port  m  is  given  by 


E  *0^  dv 


(1.12 


where  j'J'  is  the  impressed  source  current  at  port  m  and  e'^  is  the  total 
electric  field  due  to  the  N  expansion  modes  0^^,  t*l,2,...,N,  excited  by 
j'J’  with  the  other  ports  open  circuited,  e'”  can  be  expanded  in  terms 
of  E,  which  represent  the  electric  fields  due  to  expansion  modes  J.  as 
fol lows 


11 


(1.13) 


E  =  I  C  E, 


Substituting  Equation  (1.13)  into  (1.12)  and  applying  (1.6)  yields 


m 


(1.14) 


where  are  the  expansion  mode  current  amplitudes  found  from  Equation 
(1.4);  V*?  are  the  voltages  induced  at  port  m  due  to  E  .  and  V^(m) 

^  it  •  n 

denotes  the  total  voltage  at  port  m  due  to  e"’. 

The  mutual  impedance  between  port  m  and  port  n  can  be  written  as 


c"’  1"  ^ 

E  •  dv 


m  n 


(1.15) 


where  E  is  the  total  electric  field  due  to  N  expansion  modes, 
ilj,  a=l,2,...,N,  excited  by  j'!'  with  port  n  open  circuited.  Using 
Equations  (1.6)  and  (1.13)  in  (1.15)  produces 


I  C  v" 

m  n 


(1.16) 


where  V  are  the  voltages  induced  at  port  n  due  to  E  when  J.  i$  set 
to  zero  (open  circuited);  and  vP(m)  is  the  total  voltage  at  port  n  due 
to  e"*.  Via  reciprocity. 


,  for  all  m,n 


(1.17) 


It  is  common  that  the  mutual  coupling  is  measured  in  terms  of 
scattering  parameters  From  the  circuit  theory  for  waveguiding 
systems,  one  finds  that  the  scattering  matrix  [S]  is  given  by 


[S]  =  ([ZP]  +  [U])‘^  ([ZP]  -  [U])  ,  (1.18) 

where  [IJ]  denotes  the  unity  matrix. 

From  the  above  discussion  it  is  seen  that  the  port  impedances  are 
functions  of  the  impedance  matrix  [Z],  Therefore,  in  order  to  analyze 
the  mutual  coupling  between  microstrip  antennas,  it  is  essential  to  have 
full  knowledge  of  the  self  and  mutual  impedances  between  expansion  modes 
(i.e.,  the  elements  of  [Z]).  In  this  work,  the  major  effort  will  be 
focused  on  the  computation  of  these  parameters. 

D.  ORGANIZATION 

A  general  solution  to  the  dielectric  slab  problem  is  presented  in 


Chapter  II.  Numerical  results  for  specific  examples  are  obtained  and 
compared  with  measurements.  Expressions  for  the  far  zone  radiation 
fields  are  derived  and  details  of  the  calculations  are  discussed. 


In  Chapter  III,  the  structure  of  free  space  dyadic  Green's  function 
is  characterized  in  terms  of  solenoidal  and  irrotational  components.  A 
complete  eigenfunction  expansion  for  the  free  space  dyadic  Green's 
function  is  obtained.  This  expansion  is  used  in  determining  the  dyadic 
Green's  function  for  a  dielectric  coated  cylinder.  Certain 
orthogonality  properties  of  the  cylindrical  vector  wave  functions  [15] 
are  established  in  this  chapter. 

The  dielectric  coated  cylinder  problem  is  treated  in  Chapter  IV. 
Fields  due  to  current  sources  at  the  dielectric-air  interface  can  be 
obtained  in  two  ways.  One  approach  relies  on  the  technique  used  in 
Chapter  II.  The  other  is  based  on  the  dyadic  Green's  function 
determined  in  Chapter  III.  Use  of  these  fields  yields  the  impedance 
matrix  [Z].  The  method  of  steepest  descent  is  employed  to  derive  an 
expression  for  the  far  field  of  the  currents  on  a  patch. 

The  impedance  expressions  obtained  in  Chapter  IV  are  basically 
eigenfunction  solutions  which  converge  slowly  as  the  radius  of  cylinder 
increases.  In  Chapter  V,  the  Poisson  summation  formula  is  introduced  to 
convert  the  Green's  dyadic  function  to  an  alternative  representation 
which  is  more  rapidly  converging.  To  illustrate  this  approach,  a  two 
dimensional  coupling  problem  is  considered. 

Concluding  remarks  are  presented  in  Chapter  VI,  and  various 
analytical  details  are  given  in  the  appendices. 


CHAPTER  II 


MICROSTRIP  PATCHES  ON  A  fiROUNDED  PLANAR  DIELECTRIC  SLAB 

A.  INTRODUCTION 

This  chapter  deals  with  mutual  coupling  between  microstrip  patches 
printed  on  a  grounded  dielectric  slab.  For  the  calculation  of  impedance 
matrix  and  voltage  vector  elements,  the  electric  field  due  to  an 
expansion  patch  mode  is  needed.  A  general  solution  to  the  field 
equations  pertaining  to  the  microstrip  patch  on  a  grounded  dielectric 
slab  is  presented  in  Section  B.  The  surface  current  on  the  patch  is 
introduced  only  through  the  boundary  conditions;  this  simplifies  the 
calculation,  which  is  performed  in  the  Fourier  transform  domain.  The 
solution  is  essentially  of  the  Green's  function  type  even  though  the 
Green's  function  is  not  constructed  explicitly.  It  is  exact  in  the 
sense  that  both  the  dielectric  slab  and  the  ground  plane  are  taken  into 
account  rigorously.  Thus  surface  waves  and  coupling  to  adjacent  antenna 
elements  can  be  accurately  determined.  An  efficient  evaluation  of  the 
Green's  function  is  also  discussed.  A  moment  method  solution  for  mutual 
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coupling  between  rectangular  microstrip  antenna  elements  is  treated  in 
Section  C.  In  Section  D,  numerical  results  of  mutual  coupling  for  two 
antenna  geometries  are  presented.  These  results  are  compared  with 
measurements.  Finally,  in  Section  E  a  general  expression  for  the  far 
field  of  a  patch  mode  using  the  method  of  stationary  phase  is  derived. 

An  e-^*^  time  dependence  is  assumed  and  suppressed  throughout  this 
chapter. 

B.  THEORY  AND  GENERAL  SOLUTION 

The  geometry  under  consideration  is  shown  in  Figure  2.1.  The 
grounded  dielectric  slab  is  infinite  in  extent  in  the  x,y  directions 
with  uniform  thickness  t.  A  microstrip  patch  is  printed  on  the  slab  at 
the  dielectric-air  interface  with  current  density  J5.  Since  both 
regions  1  and  2  are  source  free,  an  arbitrary  field  that  satisfies 
Maxwell's  equations  can  be  constructed  from  two  scalar  functions  [16]: 
’Hni  which  generates  a  TM  field  and  '|>ei  which  generates  a  TE  field,  where 
i  =  1  for  the  region  inside  the  dielectric,  and  i  =  2  for  the  region 
outside  the  dielectric.  Both  scalar  functions  and  i|;g^  satisfy  the 
scalar  wave  equation: 


(7^  + 


'•'mi^'") 


=  0 


(2.1) 


"e  k2 
r  0 


in  region  1 


(i 


1) 


2 

where  k. 
1 


in  region  2  (i  =  2) 
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MICROSTRIP 

PATCH 


Figure  2.1.  fieometry  of  a  microstrip  patch  printed  on  a  grounded 
dielectric  slab, 

is  the  free  space  propagation  constant,  and  relative 

dielectric  constant  in  region  1. 

Upon  multiplying  and  1(1^^  by  the  unit  vector  z,  one  can 
construct  the  fields  as  [16] 

£^(0  =  -  7x(z4»ei-(r)]  +  7  x  7  x  [z^^^Cr))  , 

1 

H^(r)  =  7x(zi|»^^(r))  +  7  x  7  x  (zi|.g^{r)) 
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Ll- 


where  and  =  e^.  The  explicit  form  of  the  field  components 

can  be  found  in  Appendix  A. 

Since  the  structure  is  infinite  in  the  x-y  plane,  the  four  scalar 
wave  functions  can  be  represented  in  terms  of  their  2-D  Fourier 
transforms  (or  2-0  plane  wave  expansions)  as  follows: 

in  region  1 


_L  7, 

_'J^ei(r) 


in  region  2 


1  » 

.k  ) 
ml  X  y 

_*ez(?)_| 

•  00 

_'*'el(''x»%L 

•  ik  •  r 

e  dk^dky  , 


where 


‘zl 


V  r  0 


2  2 

k  -  k 
X  y 


'^z2  =  y  *^0  -  '^x  - 

M  A  A  A 

k  =  xk  +  yk  +  zk 


^  r> 

y  z?- 


and 


dk^dk^ 


(2.4) 


(2.5) 


,  Re  k  ,  >  0  ,  Im  k  ,  <  0  , 

zl  zl  ’ 


,  Re  kz2  >  0  .  '^z2  ^  ^  ’ 


(2.6) 


r  =  XX  +  yy  +  zz 


It  should  be  pointed  out  that  using  Equations  (2.4,5)  in  (2.2,3)  will 
automatically  satisfy: 


a.  the  boundary  condition  at  the  conducting  ground  plane: 

z  X  Ej  =  0  at  2  =  -t  ;  (2.7) 

b.  the  radiation  condition  as  r  =  |rl  >  «;  and 

c.  the  criterion  for  the  integrals  in  (2.5)  to  converge  as  z  > 

To  specify  the  fields  uniquely,  the  boundary  conditions  at  the 
dielectric-air  interface  must  be  satisfied: 


z  X  (E2  -  Ej)  =  0  at  2=0  , 

and 

^  •  m 

2  X  (H2  -  Hj)  =  at  z  =  0 


(2.8) 


(2.9) 


Without  loss  of  generality,  the  current  density  is  assumed  to  be  Fourier 
transformable  and  to  have  both  x  and  y  components.  Hence,  Jj  can  be 
written  as 


J5(x,y)  =  xJj^(x,y)  +  yJ,„(x,y) 


sy 


1  *  -j(kyX+k  v) 


dk  dk 
X  y 


(2.10) 
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Use  of  Equations  (2.8)  and  (2.9)  will  determine  completely  the 
spectral  functions  (i=l,?).  The  algebraic  details  are  carried 

out  in  Appendix  A.  Only  the  field  components  are  listed  below: 


region  1  (-«  <  x,y  <  ->  ,  -t  <  z  <  0) 
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-J 


^22^'“^  =  4ir2 


(k^+k^) 

X  y  ^ 


u>e. 


-j(k  x+k  v+k  _z) 
<ke2  ®  ^ 
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(2.20) 


^2^''^  "  47t2  [I 


k  k 

y  z2 


'^^9  +  k  't 


uUq  'qZ  ■'xW 


•j(k^x+k^+k^^z) 
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(2.21) 
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je  k  - 
r  z2 


’^"'1  "  ■  ,^Z 
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k  -Sink  ,t 
zl  zl 


^>"2  =  •  ^Vsx  ^  Vsy^ 


(kW)o  ''S’'  y  sy 
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(2.26 


He  =  k^^cosk^^t  + 


(2.27 


zr""'zi 


(2.28 


It  should  be  pointed  out  that  there  are  two  dyadic  Green's 


functions  and  5^ ,  associated  with  the  grounded  dielectric  slab. 


corresponds  to  the  case  where  the  field  points  are  inside  the  substrate 


(region  1),  and  corresponds  to  the  case  where  the  field  points  are 


outside  the  substrate  (region  2).  Only  six  components  of  each  dyadic 
Green's  function  are  considered  since  the  surface  current  on  a 
microstrip  patch  has  no  z-component.  These  components  can  be  identified 
from  Equations  (2.11)  through  (2.13),  and  (2.17)  through  (2.1<1)  as 
fol lows; 


OJU 


0 


jk*r' 


dk  dk 
X  y 


(2.29) 


where 

i  =  i,2 

A  A 

r'  =  xx'  +  yy' 

AAA 

r  =  XX  +  yy  +  zz 
a  =  x,y,z 
0  =  x.y 

and  g^g  are  obtained  from  the  electric  field  components  which  typically 
have  the  following  form; 


E.i(r)  =  JI 

•  00 

(2.30) 

It  follows  from  Equation  (2.30)  that  the  calculations  of  either 


electric  fields  or  mutual  impedances  will  invariably  involve  the 
numerical  evaluation  of  an  integral  of  the  form; 


r 


// 


F(k,.k^) 

“A 


dk  dk 
X  y 


(2.31) 


which,  however,  can  be  facilitated  by  changing  to  polar  coordinates  k,(^ 
where 


^  •  (2.32) 

=  •'Sin*  .  (2.33) 


Thus, 


2n  F(k,'t) 

r  =  J  kdk  J  -r-i  d* 

0  ‘ e  m 


(2.34) 


The  contour  for  the  k  integration  is  shown  in  Figure  2.2.  The  branch 
cuts  for  the  branch  points  k  =  ±ko  are  defined  by  the  analytic 
properties  that 


2  2 

a)  Im  k^p  =  im  ^k^  -  k  <  0  on  the  entire  top  Riemann  sheet; 


b)  Re  k^2  >  0  in  the  first  and  third  quadrants;  and 

c)  Re  k^2  <  fl  in  the  second  and  fourth  quadrants. 

However,  k  =  ±  /e^  k^  are  not  branch  points  since  the  integrand  is  a 

/  2  2 

single-valued  function  of  k^^  V^r*^o  branch  cuts  are  also 


shown  in  Figure  2.2. 


The  rienominator  in  the  integrand  of  Equation  {2.?A)  defines  the 
surface  wave  modes.  These  modes  are  determined  by  the  roots  of 


^  (for  the  TE  waves) 


D  =  0  (for  the  TM  waves) 
m 


(2.35 


(2.36 


J2  2 

e^k^-k  into  Equations  (2.35)  and 
(2.36)  and  rearranging  terms  yields,  respectively. 


/  2  2 

'a  -a  +  ocoto  =  0 


Cp  w  a  -a  +  otana  =  0 


(2.37 


(2.38 


It  can  he  demonstrated  that  if  e  is  real  and  e  1,  the  roots  of 

r  r 

Equations  (2.37)  and  (2.38)  are  real  and  located  inside  the  segemnt 
k  <k</ek.  IfN,N  represents  the  number  of  roots  for  the  TE,  TM 
case,  respectively,  then  it  can  be  shown  that 


0  ,  for  a  <  ir/2 

n  ,  for  (n-l/2)ir  <  a  <  (n+l/2)n,  n  =  1,2,3,.... 


N  =  n+1  ,  for  nit  <  a  <  (n+l)it,  n  =  0,1,2,.. 


It  is  noted  that  the  dominant  TM  mode  has  a  zero  cutoff  frequency.  Fo 
lossy  dielectric  (Im  t  0),  the  roots  just  move  off  the  real  k-axis 
with  the  form  k  =  k^-jk.,  k.  >  0.  The  exact  root  locations  can  be 
determined  by  using  the  Newton-Rhapson  procedure. 

For  efficient  numerical  evaluation  of  the  integral  in  Equation 
(2.34),  the  new  contour  is  adopted  by  deforming  (as  shown  in 
Figure  2.2).  The  integration  along  is  performed  by  computing  the 
Cauchy  value  of  the  integrals  around  the  surface  wave  poles.  Assuming 
the  surface  wave  poles  constitute  an  ordered  set  as  {p^,  *,=1 ,2,. . .  ,n} , 
Equation  (2.34)  can  be  written  as 


r  =  ?,im 
5->-o+ 


Po-'S 

p  -S 

p^  +  5 

Pn-l"' 

~F(k,*)“ 

d^i  Residue 

0  D 

e  m 
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n  2Tr 

I  / 

il  =  l  0 


In  the  case  of  lossy  dielectric,  the  integrations  from  p^_5  to 
Pj^+6,  11=1 ,2,3 ,...  ,n ,  can  be  evaluated  analytically  without  indenting 
the  contour  C^.  This  is  done  by  using  two  terms  of  a  Taylor  series 
expansion  of  D  D  about  p  ,  and  by  taking  the  value  of  the  numerator 
F(k,!(i)  at  k=p^  throughout  the  interval. 

In  actual  numerical  evaluation  of  (2.39),  6  ~  O.OOl/kg,  and  the 
infinite  integral  is  terminated  at  k  ~  150  kg. 


C.  MUTUAL  COUPLING  ANALYSIS 


In  this  section,  the  moment  method  solution  for  coupled  microstrip 
antennas  is  described.  Without  loss  of  generality,  an  array  of  M 
arbitrarily  oriented  rectangular  microstrip  antenna  patches  on  a 
grounded  dielectric  slab  is  considered.  When  patch  q  is  excited  by  a 
vertical  filament  of  constant  current  density  (a  rough  approximation 
to  a  coaxial  or  microstrip  line  feed),  an  incident  field  E?  is  radiated 
Let  the  sum  of  the  surface  currents  induced  on  each  of  the  patches  by  J 
be  denoted  as  this  3^  radiates  to  produce  the  scattered  field  E^. 

On  the  surface  of  the  patches,  the  total  tangential  electric  field 
vanishes  so  that 

I  jj  *  ^^i  ^  »  (2.41 

where  is  an  arbitrary  non-zero  test  current  (usually  called  the  test 
mode)  located  on  the  surface  S  of  the  microstrip  patches  (where 
exists).  Denoting  the  electric  field  due  to  Jj  by  and  applying 
the  reciprocity  theorem  to  the  pair  of  fields  (E^;  E^-)  and  sources  (J^: 
J^)  yields 

/  3t  •  E^  ds  =  /  L  •  3^  ds 

S  ^  5  S  ^  s 

Applying  the  reciprocity  theorem  once  again  to  the  pair  of  fields 
(E?;  E^)  and  sources  (J^;  J^)  yields 


It  then  follows  from  the  above  reciprocity  relations  and  Equation  (2.4i) 
that 


/  E,  •  ds  =  /  •  J?  dv 

^  T  s  ^  T  1 


(2.42) 


The  next  step  in  the  moment  method  solution  is  to  expand  in 
terms  of  known  basis  functions  from  a  finite  ordered  set 

"  d  ^  Cl  ^01 

“  ~  ”  ^^nx  ^  '^'^ny*^  ~  =  1,2,. ..,M,  +...+  =  N| . 


M  N, 


I  I  i!(q)  j! 


a=l  n=l 


(2.43) 


where  J^,  known  as  an  expansion  dipole  mode,  is  the  current  density  of 
mode  n  on  patch  a,  and  I^(q)  is  its  unknown  current  amplitude.  In  the 

“  “A 

Galerkin  scheme  J,  is  chosen  to  be  an  element  J  eH  which  excites  the 
T  m 

electric  field  E^^.  Hence,  inserting  Equation  (2,43)  into  (2.42)  for  all 
m  and  S  yields 


Na  a  a 

J:  J: 


(2.44) 


where  m  =  1,2,. ...Ns,  3  =  1,2,., 


,  ,M,  5)  N3  =  N,  and 
S 


3a  (  "3, -8a  -a. 

=-JE(r  +r)*J(r  )ds  , 
mn  ^  m'  mn  n'  n'  n' 


’'.*('1)  =  /  E!(r!)  •  j?(r!)d' 


(2.45; 


(2.46) 


E 


-3a  8a  Sa 

^mn  "  dTsplacement  vector  from  the  center  of 

mode  m  on  patch  3  to  the  center  of  mode  n  on 
patch  a,  and 

“6  S3 

r*^  =  (x^,y^,z)  is  the  postion  vector  referenced  to  the  center  of 
mode  m  on  patch  3. 

In  Equation  (2.46),  the  excitation  3^  at  port  q  is  iDOdeied  as 
oJ(r®)  =  zlj6(x-x®j)5(y-y®j),  -«  <  x,y  <  ».  -t  <  z  <  0,  where 
(X|^^,y^^,0)  is  the  displacement  vector  from  the  center  of  mode  m  on 
patch  3  to  the  feed  location  on  patch  q. 

The  geometry  of  modes  m  and  n  pertaining  to  the  evaluation  of 
Sex  S 

^mn  depicted  in  Figure  2.3.  In  terms  of  matrix  notations, 

(2.44)  can  be  written  as 


[Z]Cl(q)]  =  [V(q)] 


where 


CZ]  = 


a^*’  column  block 


'  ,3a 

,3a  ' 

1  11  *  * 

•  "in  1 

1  • 

a  1 

• 

• 

• 

1  • 

.  •  1 

• 

.  • 

1 

.  1 

N-1 

N-N 

1  ^ 

3  a 

(2.47) 


3^*^  row  block 


(2.48) 


-.vs 


[I(q)j  =  [i;{q)...i;^(q)  ...  IjCq) . . .  l"  jq)  ...  I  ^  (q) . . .  I^j^Cq)  ] 


[V(q)J  =  [vj(q)...vj^(q)  ...  Vj(q)...v“  (q)  ...  Vj(q) . . .vJJ^(q) ]  . 


In  Equations  (2.49)  and  (2.50),  T  denotes  taking  the  transpose  of  the 
row  vector.  Thus,  the  current  amplitudes  [I(q)]  can  be  solved 
algebraically  by  matrix  inversion: 


[I(q)l  =  [Z]"^[V(q)] 


As  indicated  in  Equation  (2.51),  the  current  amplitudes  will  depend  on 
the  feed  location  on  patch  q. 

According  to  Equation  (1.14),  the  input  impedance  at  port  q  is 


given  by 


M  N. 


I  1 


a=l  n=l 


Similarly,  it  follows  from  Equation  (1.16)  that  the  mutual  impedance 
between  port  q  and  port  r  is  given  by 


M  Nq 

I  I 

a=l  n=l  ^ 


In  Equations  (2.52)  and  (2.53),  it  is  understood  that  and  L  are  the 
terminal  currents  at  port  q  and  port  r,  respectively;  the  Ip(q)  are 
obtained  from  (k:.51);  v“(q)  and  v^(f*)  the  voltages  at  port  q  and 
port  r,  respectively,  due  to  mode  n  on  patch  a  excited  by  source  when 
all  other  ports  are  open  circuited. 

To  compute  the  moment  method  solution  for  mutual  coupling  between 
microstrip  patch  antennas,  one  needs  to  evaluate  Equations  (2.45)  and 
(2.46).  As  presented  in  Appendix  B,  the  exact  expressions  for  and 
v^(q)  are  given  by 


'^x''z2 


73  .  !  70 
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01  a 


(2.62) 

is  the  angle  between  the  x-axis  of  mode  m  and  the  x'-axis  of  mode  n 
3  S 

(see  Figure  ?.3);  2x  and  2y  are  the  widths  of  trode  m  on  patch  0  in  the 

X  and  y  directions,  respectively,  and  likewise,  2x®,  2y“  are  the  widths 

n  n 

of  mode  n  on  patch  ti  in  the  x',y'  directions. 

D.  NUMERICAL  EXAMPLES 

Calculation  of  mutual  coupling  ($12  parameter)  between  two 
identical  coax-fed  rectangular  microstrip  antennas  is  presented  in 
this  section.  Two  antenna  geometries  (E-plane  and  H-plane  coupling)  are 
considered.  In  both  examples,  the  expansion  set  is  chosen  to  be 

^  "  ^^n  "  ^^nx*  Nj=N2=1  }  ,  (2.63) 

where 

sin<(x“-|x| ) 

2yois'jnKX®  *  ^  ^  ■^1’  0-1,2 

0  ,  otherwise  .  (2.64) 

In  Equation  (2.64),  the  coordinates  x,  y  are  referenced  to  the  center  of 
the  expansion  dipole  mode  of  length  2x“  (along  the  direction  of  the 
current)  and  width  2y^.  Furthermore, 
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1*  V 


<  = 

0  eq 


^eq  ^0 


e  +  1  e  -  1 
r  r  lot 
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(2.65) 


(2.66) 


The  choice  of  dipole  mode  currents  and  constant  ic  is  suggested  by 
[17].  It  should  be  remarked  that  a  current  expansion  mode  in  the  set  H 
is  employed  to  approximate  the  even  (dominant)  current  mode  (Figure 
2.4(a))  on  a  isolated  microstrip  patch  at  its  first  resonance.  As 
indicated  in  Equation  (2.64),  such  an  expansion  mode  is  uniform  in  the 
y-direction  and  piecewise  sinusoidal  in  the  x-direction.  It  is 
understood,  nevertheless,  that  at  least  one  odd  (attachment)  mode 
(Figure  2.4(b))  also  exists  simultaneously  at  resonance  to  account  for 
the  discontinuity  in  current  due  to  the  feed.  Summing  the  even  and  odd 
modes  will  reasonably  represent  the  true  current  distribution  on  the 
patch  (Figure  2.4(c)).  For  the  sake  of  simplicity,  however,  only  the 
even  dominant  mode  is  considered  in  the  calculations  with  the 
understanding  that  at  resonance  only  the  even  mode  will  dominate  and 

closely  resemble  the  true  current  distribution. 

12  12  12 
Since  x^  =  x^,  y^  =  y^  (identical  patches),  and  =  0, 

Equations  (2.54)  and  (2.55)  reduce  to,  respectively, 

_  k  k 
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(2.67) 


2  2  2 


(a)  EVEN  (DOMINANT)  (b)  ODD  (ATTACHMENT)  (c)TRUE  CURRENT 
MODE  MODE 
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-j[x(k^cos*jj+k^sin^jj)  +  yC-k^sin^^j+lc^cos^j)] 
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(2.72) 


(2.73) 


It  should  be  noted  that  (x”*,  is  determined  from  Equation  (2.68)  by 

enforcing  the  condition  that  Vj(a)/lf  is  a  real  quantity.  This 
corresponds  to  the  condition  that  the  isolated  microstrip  antenna  is  at 
resonance. 

Example  1:  E-plane  coupling 

The  antenna  geometry  is  shown  in  Figure  2.5.  In  this  special  case, 
1  2 

Xj  s  =  3.275  cm 

yj  =  y^  =  5.285  cm 

11  22  ,  ... 

’^If  ”  ~  ” 

yj‘  -yjf  -0  cm 

♦ll  ■  “ 


y^^“  =  0  ,  S,o=1.2 

frequency  =  1417  MHz 
t  =  0.1575  cm,  and 


|s,j|  dB 


ij'.i  u  ■■  jp  ym  V  1  'j  •  'j  ■  V  •  'J  ■  .'  m  4  i  u  J'  «'■  u  ^ 


Also,  the  even  and  odd  properties  of  the  integrands  can  be  used  to 
reduce  the  domain  of  the  a  integration  from  [O.Zn]  to  [0,n/2]. 
Equations  (2.74)  and  (2.75)  can  then  be  written  as 
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(2.80) 


^  «  ir/2  ‘'x‘'z2 

v,(a)  =  ^ZTiT  i  !  t  n  Sink  t  sink  x  ,J  kdadk  , a, 3=1, 2  . 

^  0  0  0  '^zr’m 

(2.81) 


Equations  (2.80)  and  (2.81)  are  the  final  forms  used  for  computing  the 
impedance  matrix  and  voltage  vector  elements. 


Example  2;  H-plane  coupling 


Figure  2.6  shows  an  antenna  configuration  for  studying  the  H-plane 
coupling.  In  this  case. 


1  2 

Xj  =  x^  =  3.275  cm 
yj  =  =  5.285  cm 


1.115  cm 


Figure  2.6.  H-plane  coupling  for  identical  microstrips 


frequency  =  1417  MHz 


t  =  0.1575  cm,  and 


e 

r 


2.5  . 


As  noted  in  Example  1,  a=l,2.  Moreover,  it  follows  from 

Equation  (2.73)  that  (Xj“,  y^“)  =  (x^,  )  *  a,B=l,?.  Then  Equations 

(2.67)  and  (2.68)  can  be  simplified  to 
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Vj(a)  =  - 
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a, 6=1, 2 


(2.83 


Introducing  polar  coordinates  k  and  a  in  Equations  (2.82)  and  (2.83)  to 
obtain 
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(2.R5) 


From  Figures  2.5  and  2.6,  one  can  observe  that  the  E-plane  coupling 
is  larger  than  the  H-plane  coupling.  This  is  due  to  a  stronger  surface 
wave  which  is  excited  for  the  E-plane  configuration.  In  general,  the 
magnitude  of  the  scattering  parameter  sj^  is  seen  to  have  good  agreement 
with  measurements  by  Jedlicka,  et  al .  [18].  For  the  data,  the  isolated 
microstrips  were  resonant  with  50-1  resistance.  This  occured  at 
1417  MHz  for  the  computations,  and  at  1410  MHz  for  the  measurements.  It 
is  noted  that  the  percentage  tolerance  on  the  substrate  permittivity  and 
patch  size  are  of  the  same  order  of  magnitude  as  the  antenna  bandwidth 
(typically  a  few  percent).  Such  an  error  causes  a  shift  in  resonant 
frequency,  hut  it  is  found  that  the  change  in  the  calculated  terminal 
impedance  is  negligible. 


E.  FAR  ZONE  RADIATION  FIELD 


I 

The  radiation  field  refers  to  the  electric  field  in  region  2  when 
the  field  point  is  significantly  removed  from  the  source.  From 
Equations  (2.17)  to  (2.19),  it  can  be  shown  that  the  electric  field  may 
be  expressed  by  the  2-0  Fourier  transform  as  follows: 


-  _  1  “  =  -jk*r 


(2.86) 


where 


k  =  xk^  *  yky  ♦  zk^2'  "22  '  '^22  "  ''22  " 


r  =  XX  +  yy  +  22  ,  -«  <  x,y  <  «  ,  0  <  z  <  «  (region  2), 


E,(k^.ky)  =  xE^j(k^,ky)  *  yE^2(^.k^)  *  ^E^^(k^.k,). 
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Since  v  •  E2(r)  =  0  in  region  2,  the  divergence  of  Equation  (2.86) 


gives 


3  -j^*r 


-jk»r 


-jk»r 


7  •  E^e  =  E^  •  7e  =  -jk  •  E^e  =  0 


(2.87) 


Thus  only  two  components  of  E^  are  independent.  Then  E^  can  be 
expressed  by 


^2  '  *  ^^2  -  ^  k,j  "'k^2  *  Vy2> 


(2.88) 


As  shown  in  Appendix  C,  the  integral  in  Equation  (2.86)  can  be 
evaluated  asymptotically  by  the  method  of  stationary  phase  as  r+*. 

Using  Equation  (2.88)  and  (C.15)  in  (2.86)  yields  the  far  zone  radiation 
field  as  follows: 
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(2.89) 


where  =  (k^sin9cos<^,  k^s1n9sin4>)  is  the  stationary  phase 
point,  k^2  =  k^cose,  and  9  and  <j>  are  the  usual  angle  variables  in 
spherical  coordinates. 


In  terms  of  spherical  components.  Equation  (2.89)  becomes 


“■''‘o  2»r  l*[cos*  EyjC'xo-V’  ■ 

*  e(cos»  E^^C'xo’So'  *  'yj'^xo-^yo'))  •  (E-’n) 


Once  the  current  density  on  the  microstrip  antenna  is  determined,  the 
radiation  field  can  be  obtained  from  Equation  (2.90). 
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CHAPTER  III 


DYADIC  GREEN'S  FUNCTIONS 


A.  INTRODUCTION 

For  many  electromagnetic  problems,  the  eigenfunction  expansion 
method  provides  the  most  useful  solution  to  the  (electric  or  magnetic) 
field;  the  boundary  conditions  of  the  problem  are  readily  incorporated 
into  the  solution.  For  waveguides  and  cavities,  in  particular,  the 
eigenfunction  expansion  method  is  usually  the  preferable  way  to  solve 
problems  involving  current  sources  and  wave  scattering.  The  whole 
procedure  for  finding  the  eigenfunction  expansion  of  the  electric  field 
can  be  organized  in  a  systematic  way  by  introducing  the  electric  dyadic 
Green's  function  or  simply  Green's  dyadic  (Gg)  which  is  discussed  in 
this  chapter. 

s 

The  eigenfunction  expansion  of  Gg  had  been  a  subject  of 
misconception  and  controversy  in  the  past.  The  difficulty  centered 
about  the  completeness  of  the  expansion.  Initially,  Tai  developed 
expansions  for  the  dyadic  Green's  functions  in  his  book  [19]  which 


Included  only  the  M  and  N  functions  which  are  the  solenolda! 

(transverse)  Hansen  vector  wave  functions.  After  the  publication  of  his 
book,  Tal  noted  that  his  expansions  were  Incomplete  In  the  source  region 
and  discovered  the  need  for  the  Irrotatlonal  (longitudinal)  Hansen 
vector  wave  function  L.  In  his  subsequent  amendments  [20,21],  Tal 
Indicated  how  the  complete  expansions  of  the  dyadic  Green's  functions 

M  « 

could  be  obtained  via  the  method  of  Gm,  where  G^,  Is  the  magnetic  dyadic 
Green's  function  with  zero  divergence  and  may  be  constructed  In  terms  of 

s  a 

only  the  solenoldal  eigenfunctions.  Ge  Is  then  obtained  from  v  x  G„, 
together  with  a  dyadic  delta  function  term  at  the  source  point.  Thus, 
he  avoided  having  to  deal  with  the  L  functions  explicitly.  About  the 
same  time,  Collin  [22]  discovered  Independently  the  lack  of  completeness 
for  the  E  and  H  modes  In  the  same  waveguide  problem  considered  by  Tal. 

a 

Later,  Rahmat-Samll  [23]  modified  Tal's  method  of  Gm  via  the  use  of 
symbolic  functions  (distributions).  More  recently,  Collin  [24] 
successfully  unified  the  various  results  on  the  dyadic  Green's  functions 
as  presented  by  different  authors.  Various  representations  for  the 
Green's  dyadic  were  presented  and  their  Inter-relationships  pointed  out. 
Most  recently,  Pathak  [25]  described  a  relatively  simple  method  for 

a 

constructing  a  complete  expansion  of  Gg  In  a  compact  and  useful  form 
which  contains  only  the  solenoldal  type  eigenfunctions  plus  an  explicit 
dyadic  delta  function  term  that  accounts  for  completeness  at  the  source 


A  main  feature  of  this  chapter  is  on  the  eigenfunction  expansion  of 
Ge  in  free  space.  The  present  approach  follows  that  of  [25]  in 
obtaining  the  source  conditions  that  govern  the  Green's  dyadic  in  free 
space.  It  then  conjectures  that  the  free  space  Green's  dyadic  can  be 
expressed  as  the  sum  of  a  solenoidal  component  and  an  irrotational 
component  which  is  allowed  to  be  a  symbolic  function.  The  solenoidal 
component,  as  usual,  is  expanded  in  terms  of  the  M  and  N  functions;  and 
the  irrotational  component  is  deduced  from  the  differential  equation 
governing  G0.  Of  course,  the  irrotational  component  turns  out  to  be  in 
the  form  of  a  dyadic  delta  function  as  in  [20,21]  and  [25].  Making  use 
of  the  free  space  Green's  dyadic  and  the  principle  of  scattering 
superposition,  the  dyadic  Green's  function  for  an  infinite  dielectric 
coated  conducting  cylinder  can  be  determined  without  additional 
analytical  difficulty. 

The  remainder  of  this  chapter  is  organized  as  follows.  Section  B 
reviews  some  basic  relations  governing  the  electric  field  (E),  the 
magnetic  field  (H),  Gg  and  G^.  Only  E  and  Gg,  however,  will  be  used 
in  the  chapters  that  follow.  The  free  space  Green's  dyadic  is 
characterized  and  expanded  in  Section  C  which  also  establishes  the 
orthogonality  of  the  M  and  N  functions  in  cylindrical  coordinates.  In 
Section  D,  the  Green's  dyadic  for  a  dielectric  coated  cylinder  is 
constructed.  Finally,  Section  E  expands  the  Green's  dyadic  into 
component  form  from  which  the  electric  field  is  easily  identified  when 
the  source  is  known. 
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In  the  remaining  portion  of  this  report,  all  time-harmonic 
quantities  are  assumed  to  have  an  e“j“^  time  dependence  which  will  be 
suppressed  throughout. 

B.  SOME  BASIC  RELATIONS  GOVERNING  E,  H,  Gg  AND  G^, 

The  electromagnetic  fields  E  and  H  excited  by  the  electric  current 
source  J  in  a  region  V  of  linear,  isotropic  and  time-invariant  medium 
satisfy  the  Maxwell's  equations: 


ja)uH(r)  =  '7  X  E(r)  , 

(3.1) 

k2E(r)  =  ja)u(7  X  H(r)  -  J(r))  , 

(3.2) 

7  •  H(r)  =  0  , 

(3.3) 

and 

-  _  p(r)  ju)M  _ 

7  •  E(r)  =  =  -  7  .  [-p-  j(f))  ,  (3,4) 

where  r  e  V,  u  and  e  denote  the  constitutive  parameters  of  region  V; 

and  p  is  the  charge  density  associated  with  J.  One  can 
verify  from  Equations  (3.1)  and  (3.2)  that  E  and  H  satisfy  the  following 
differential  equations,  respectively: 


7  X  7  X  E{r)  -  k2E(r)  »  j(DuJ(r) 


f 


(3.5) 


and 

7  X  7  X  H(r)  -  k2H(r)  =  7  x  J(r)  .  (3.6)  ^ 

I 

The  electric  Green's  dyadic  Gg  for  region  V  is  defined  as  the 
solution  of  the  vector  wave  equation: 


7  X  7  X  Gg(r.r')  -  =*  !5(r-r')  ,  (3.7) 

and  the  boundary  conditions  pertaining  to  V,  where  r,r'  c  V  denote, 
respectively,  the  field  observation  and  the  point  source  locations;  I  is 
the  unit  dyad;  and  6(r-r')  is  the  Dirac  delta  function.  It  can  be 
proved  that,  for  region  V,  there  exists  another  dyadic  function  known  as 
the  magnetic  Green's  dyadic,  denoted  by  Gp,,  which  satisfies  the 
differential  equation: 

7  X  7  X  6^(r,r')  -  k2G^(F,F')  =  7  x  I6(r-r')  (3^8) 

together  with  the  associated  boundary  conditions,  and  is  related  to  Gg 
as 


G„(r,r')  =  7  X  §  (r,r') 


(3.9) 


From  Equations  (3,7)  and  (3.9)  it  can  be  readily  verified  that 

k^Sg(r,r')  =  7  X  G|j|(r,r')  -  I6(r-r')  .  (3.10) 

It  is  also  clear  from  (3.9)  and  (3,10)  that 

,3.H, 

and 

7  *  Sg(r,r‘)  =  -7  •  [I6(r-r')/k2]  *  -  76(r-r')/l(^  .  (3.12) 

Comparing  Equations  (3.1)  -  (3.4)  with  (3^.9)  -  (3.12),  it  is  seen  that 
^e.  G|n/ju>u  and  I6(r-r')/juu  represent  the  dyadic  analogs  of  E,  H  and  J, 
respectively.  It  may  be  noted  that  relations  (3,9)  and  (3.10)  must  be 
modified  should  r  and  r'  be  located  in  different  regions  each  with 
distinct  permeability.  This  will  be  discussed  further  in  Section  D, 

It  is  of  interest  to  consider  that  E  and  Gg  are  solutions  in  region 
V  to  the  differential  Equations  (3,5)  and  (3,7),  respectively, 
pertaining  to  the  geometry  depicted  in  Figure  3.1.  Thus  E  and  Gg  must 
satisfy  the  following  boundary  conditions: 


on  the  perfect  conducting  surface  S, 


Figure  3.1,  Region  V,  bounded  by  a  perfect  conducting  surface  S  an^ 
a  surface  at  infinity,  contains  an  electric  source  J. 
V  can  be  inhomogeneous. 


on  the  surface  S«,  at  infinity. 


tim  r  7  X  E(r)  -  jkr  x  E(r)  *  0  ,  (3^1 


tim  r 


r+<» 


^  ^  6p{r,r')  -  jkr  x  6^(r,r') 


(3.1 


The  latter  set  of  conditions  is  commonly  known  as  the  radiation 
condition.  The  region  V,  bounded  by  surfaces  S  and  S*,  contains  the 
source  J  and  can  be  inhomogeneous.  Next  employing  the  Green's  theorem 
for  the  pair  E  and  Ge  over  V,  and  enforcing  the  appropriate  boundary 
conditions  on  E  and  Gp  yield 


(3.17) 


B 


I 


E(r)  =  juy  /  G  (r.r")  •  J(r")dv",  r,r"  e  V 


From  (3.1)  and  (3.9),  it  can  be  verified  that 


H(f)  =  /  Gn,(r.r'‘)  •  3(?")dv".  r,r“  e  V 


(3.18) 


In  particular,  if  J(r")  is  an  arbitrarily  oriented  electric  point  source 
of  strength  pg  located  at  r"=r',  i.e.. 


J(r")  =  Pg6(r”-r') 


(3.19) 


then  (3.17)  and  (3.18)  will  reduce  to 


E(r)  =  juju  Gg(r,r')  •  p^ 


(3.20) 


and 


H(f)  *  G„(r,r')  •  5 


(3.21) 


C.  EIGENFUNCTION  EXPANSION  OF  THE  FREE  SPACE  GREEN'S  DYADIC 


The  analysis  of  mutual  coupling  between  microstrip  patch  antennas 
involves  only  the  electric  field  and  electric  current.  It  is  sufficient 
to  consider  the  electric  Green's  dyadic  (GgCr.r'))  because,  as  observed 
in  the  previous  section,  Ge(r,r')  is  directly  related  to  the  electric 
field  due  to  an  electric  current  point  source.  The  mangetic  Green's 
dyadic  Gp,(r,r')),  if  needed,  can  be  obtained  by  taking  the  curl  of 
Ge(r,r'). 

This  section  is  devoted  to  the  eigenfunction  expansion  of  the  free 
space  Greeen's  dyadic  which  will  be  denoted  by  Gjj(r,r').  Without  loss 
of  generality,  it  is  assumed  that  the  electromagnetic  problem  under 
consideration  can  be  formulated  in  a  general  orthogonal  curvilinear 
coordinate  system  (u^,  u^,  u^)  with  corresponding  metric  parameters, 

AAA 

^1*  ^2’  ^3’  associated  unit  vectors  u^,  u^,  (Figure  3.2).  This 
requires  that  the  scalar  wave  equation  is  separable  in  this  system. 

For  the  sake  of  being  specific,  it  is  assumed  that  coordinate  uj 
satisfies  the  following  two  criteria: 


C.l.  Any  set  S  =  {(ui, 02,03):  uj  »  pg,  pg  fixed,  03,03  arbitrary} 
forms  a  closed  surface,  i.e.,  surface  S  partitions  the  space 
into  two  distinct  regions,  V+  *  {(01,02,03):  oi>Pg, 

02,03  arbitrary}  and  V-  =  {(01,03,03):  oi<Pg,  03,03  arbitrary}. 
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Figure  3.2.  A  general  orthogonal  curvilinear  system  with  coordinates 

AAA 

Un  Ug,  Uj  and  their  corresponding  unit  vectors  Uj,  u^, 
and  metric  parameters  hj,  h^,  h^. 

C.2.  The  ui'propagating  solutions  of  the  scalar  wave  equation  can  be 
obtained  so  that  the  construction  of  the  u^-propagatlng  vector 
solenoldal  eigenmodes  or  wave  functions,  denoted  by  Mq  and  Nq, 
Is  possible. 

It  should  be  remarked  that  at  any  point  on  S,  the  outward  unit  vector 

A 

normal  to  S  coincides  with  Uj.  Also,  the  shorthand  notation  q  stands 
for  a  particular  combination  of  various  mode  characteristics  such  as 
evenness,  oddness  and  order. 


vv  ■  .  ^  — T  wj~  -j-  ^  ' ’v  "v  ,  "v  I  TV  .  n,  .  /  :rv  "v.  ^  - 


It  is  of  significance  to  consider  an  arbitarily  oriented  point 

•  III 

source  J  of  strength  p^,  located  at  r'  *  (u^,  u^),  radiating  in  free 

space  Vo.  According  to  criterion  C.l,  is  partitioned  into  two 
distinct  regions  (Figure  3.3): 

^  _  • 

V  =  {r  =  (Uj,U2,U2):  Uj>Uj^,  af*^>^trary}  , 


and 

I 

V"  =  (r  =  (Uj.u^.Uj):  Uj<Uj^,  U2»‘^3  arbitrary} 

The  interface  between  V*"  and  V  forms  a  dosed  surface  defined  by 

I 

S  =  {f  =  (u^.Ug.Ug):  Uj=Uj,  u^.Uj  arbitrary}  , 

which  contains  the  source  represented  by 

6(u^-uJ)6(u2-U2)6(u3-U2) 

0(f)  =  PeS(?-;')  »  .  (3.22) 

From  (3.12)  it  is  clear  that 


-76(r-r')/ko 

0 


,  r  e  S 

,  r  e  or  F  e  V 


(3.23) 


2 

where  k©  =  (^^uo^o*  ^'o*  ^o  are  the  constitutive  parameters  of  free 

space  Vo.  Equation  (3.23)  simply  says  that  Go(r,r')  is  solenoidal 
in  both  V'*'and  V",  but  irrotational  on  S.  Thus  one  can  conjecture  that 


i 

i| 


•‘1 
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POINT  SOURCE  J (r)sp  6 (r-  ?) 


/  0  / 

I  REGION  V“:  j  REGION  U,  >  u', 


SURFACE  S  :  u,«u', 


Figure  3.3.  Free  space  V  Is  partitioned  Into  regions  v'*’  and  v‘  with 


closed  surface  S  as  Interface  which  contains  the  point 
•  •  •  •  • 
source  J  at  -  (ui.u-.u,). 


%{r,r')  »  U(uj-Uj)gQ(r,r')  +  U(uj-Uj)gjj(r.r' )  -  ♦(r.r*) 


(3.24 


where 


go(r,r')  =  G^(r,r') 


for  r  c  V 


(3.25 


9o(r,r')  =  6  (f.r')  for  r  e  v‘ 


(3.26 


u(0  * 


1 

0 


.  5  >  0 

(unit  step  function)  » 

,  5  <  0 


and  '|'(^,F')  is  a  dyadic  symbolic  function  (distribution)  which  is  used 
as  if  it  was  an  ordinary  function  with  the  understanding  that  it  will 
produce  valid  results  under  the  integral  sign  (in  the  present  case,  the 
volume  integral  over  a  subset  of  Vq), 

It  should  be  recognized  that,  by  definitions,  g^  and  g”  are 
solenoidal  (divergenceless)  in  V"*"  and  V",  respectively,  while  ^  is 
irrotational  on  S.  Mathematically  these  mean 

7  •  9”(r,r')  »  0  for  r  e  V*  ,  (3^27 

and 

V  •  ^(r,r‘)  *  0  for  r  e  S  . 


It  may  be  mentioned  that  (3.28)  Implies  Is  non-zero  In  the  sense  of 
distributions.  In  addition,  g*  satisfies  the  homogeneous  vector  wave 
equation: 

7xvxg*(r,r')  -  kogo(r,r')  »  0  r  e  V*  .  (3.29) 

and  the  radiation  condition: 


(3.30) 


Now  from  (3.20)  and  (3.21),  It  is  seen  that  the  electromagnetic 
fields  apart  from  the  surface  S  can  be  defined  as 


E“(r)  »  j%gQ(r,r')  •  p^ 


,  r  e  V 


(3.31) 


H*(r)  *  7  X  gQ(r,r')  •  p^ 


,  r  e  V 


(3.32) 


Since  the  tangential  electric  field  Is  continuous  at  S,  It  requires  that 


-+  - 


-  «(ui-ui) 

Uj  X  [E  (r)  -  E  (r)]  - -  -  0 


or 


.•+  - 


,  .  ^  *-  -  -  ,  «(ui-ui) 

“l  X  [go(r,r')  -  g-(r,r')l  fT 


»  0 


(3.33) 


Equation  (3.33)  Is  obtained  by  using  (3.31)  and  the  fact  that  "pg  Is 
arbitrary.  Likewise,  by  enforcing  the  boundary  condition  on  the 
magnetic  field  at  the  Interface  S  one  gets 

I 

-+  -  .  fi(ui-ui) 

Uj  X  [H  (r)  -  H  (r)]  *  ^t  *  ^  » 

(3.34) 

s 

where  denotes  the  transverse  part  of  the  unit  dyad  with  respect  to 

A 

the  propagation  direction  u^,  and  Is  defined  by 

S  S  AA  S  AA  AA 

I  =  It  +  UjUj  ;  =  UgUg  +  U3U3  .  (3.35) 

Employing  (3.32)  In  (3.34)  and  the  fact  that  p^  Is  arbitrary  lead  to 

I 

»+  -  -  « -  «(ui-ui)  -  _  . 

X  [V  X  gQ(r.r’)  -  V  x  gQ(r.r’)  - -  »  It6(r-r')  .  (3,35) 

Equations  (3.33)  and  (3.36)  constitute  the  boundary  conditions  on 
surface  S. 

Next  the  Irrotatlonal  component  |  will  be  determined  by  enforcing 
^o(r,r')  to  satisfy  the  Inhomogeneous  vector  wave  equation  described  In 
(3.11)  (with  replaced  by  k^).  Taking  the  curl  of  6^  In  (3.24)  gives 

^  ^  GQ(r.r')  =  U(Uj-Uj)7  X  gQ(r,r’)  +  U(Uj-Uj)7  x  gQ(r.r') 

I  "  I 

=+  .  .  = -  ,  «(ui-Ul)  «  .  _  6(ui-ui) 

+  ^^1  X  [gQ(r,r')  -  gQ(r,r')]  jq  -  7  x  i|»(r,r') - — - 

(3.37) 
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■i-'.  -r. 


The  third  term  on  the  R.H.S.  of  (3.37)  vanishes  because  of  (3.33)  and  It 

m 

Is  permissible  to  choose  ti)  such  that 


*0  ,  r  e  S  ,  (3,38) 

since  up  to  this  point  <|i  is  still  arbitrary  other  than  the  requirement 

I 

■  =  _  _  ^(U1“U]^) 

that  v*,|»  ^  0  on  S.  It  may  be  noted  that  ♦(r,r') - r - is  a 

symbolic  function,  hence  (3.38)  means  Us  curl  vanishes  on  S  in  the 
sense  of  distributions.  Consequently,  one  will  find  (3,37)  reduces  to 

^  X  S^(r,r')  *  U(Uj-Uj)7  x  ^^(r.r')  +  U(Uj-Uj)7  x  g’(r,r')  . 

(3.39) 

Taking  curl  of  both  sides  of  (3.39)  yields 


7  X 


-.,«(ui-ui) 
'l»(r,r') - r - 


7x7xG^(?,?')  =  U(Uj-Uj)7x7xg^(r.r')  +  U(Uj-Uj)7x7xg^(r,r' ) 

I 

+  Uix[vxgjj(r,r’)  -  7xgg(r,r’)]  .  (3.40) 

From  (3.29)  and  (3.36)  it  can  be  verified  that  (3.40)  will  become 

A  A 

=_>  2  '»+  _  '  «---  “1^1  -  - 

7x7xG^(r,r')  “  kg[u(uj-Uj)gQ(r,r' )+U(uj-Uj)gQ(r,r' )  “  "j^ 


+  I5(r-r') 


(3.41) 


Comparing  (3.24)  and  (3.41)  suggests  that  choosing 


i|»(r,r')  = 


6(u2-u')6(u3-u;)  ^  ^ 

Vl 

^0  ^2^3 


will  make 


=  __  •«+__  •  s _  ujui  _  _ 

Go(r.r*)  =  U(Uj-Uj)gQ(r,r')  +  U(u^-Uj)g^(r,r' )  -  -j^«(r-r') 


a  valid  representation  of  the  free  space  Green's  dyadic  provided 
conditions  (3.23),  (3.28)  and  (3.38)  are  satisfied.  To  verify  these 
conditions  one  may  Integrate  7»\|i  over  region  V-us  and  then  apply  the 
Gauss'  theorem: 


/  7»4)(r,r')dv  = 

V-US 


/  'l»(r,r')»Ujh2h3du2dU3  »  ^ 

5  ^ 


Equation  (3.44)  indicates  that  V*i|i(r,r')  does  not  vanish  on  S. 

I 

,=  4(U1-Ul) 

Likewise,  Integrating  V  x  [<|>(r,r')  - - ]  over  region  V  US  and 

applying  Stokes'  theorem  show 


*  5(ui-ui) 

/  7  X  [t|»(r,r') - ^ - ]dv 

V-US  ^ 


6(ui-ui) 

/  Uj  X  [♦(r,r’)  T  ]h2h3du2du 


I  ^  [~TT~  4(f*-?')]  h,h,du,du,  »  0 


-■  f 


If  m  , .  I ,  ,  .i|  >,1  >1^  v>  ■  ^!  ■  I  ■  L.  •  I.  ■  ij  ■  ■■  I  u  ■  I  ■  1 1  1 w  ".vw  «.'■  r 


Q 


=  «(ui-ui) 

Equation  (3.45)  Implies  that  7  x  - ]  indeed  vanishes  on  S. 

s 

Thus  the  validity  of  the  representation  of  Gq  by  (3,43)  has  been 
established.  Moreover,  this  representation  proves  to  be  most  useful  and 
compact  for  the  eigenfunction  expansion  of  Gq  which  is  considered  next. 

s 

Anticipating  the  Gq  will  be  needed  in  the  construction  of  the 
Green's  dyadic  for  a  dielectric-coated  circular  cylinder,  one  would  look 
for  the  eigenfunction  expansion  of  Gq  in  cylindrical  coordinates  (p,4,z) 
which  are  natural  for  cylindrical  structures.  To  do  so,  one  of  the 
coordinates  p,  <^,  z  must  be  designated  as  the  "u^  coordinate"  that 
satisfies  the  criteria  C.l  and  C.2  stated  earlier. 

It  is  relatively  easy  to  verify  that 

1.  the  set 

S  =  {r=(p,(fr,z):p*p* ,  p'  fixed,  forms 

an  infinite  circular  closed  surface  that  divides  the  entire 
(free)  space  into  two  distinct  regions  denoted  by 

V'*’  =  {r  =  (p,t,z):  p>p',  0<4<2x,  -od<z<«»}  and 

V-  =  {r  =  (p,ifr,z):  p<p',  0<(><2ii,  -•<z<»}  ; 

2,  the  "p-propagating"  solenoidal  vector  eigenmodes  or 
wave  functions  which  satisfy  the  homogeneous  vector  wave 
equation  are  found  be  [15,19] 
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Memx'(h',r)  =  ^  p  f,  1^{X  p){cos'"*}  -  *  p) icin’"*! 


3p  m' 


Nemx'ih'.r)  =  7  X  Me^^.(h'.F) 

0  0 

1  jh'  3  +  rCOS  .  “  +  cin 

■  !  1“  ko  3p  yk‘»Hs,„i«*!  '  tkoP  I‘(k'p)l„5«i»} 


where 


=  0,1,2,...;  -c»<h'<«o;  X'  =\ 


*^o  "  ^  0 


z-(x'p)  - 


hI^^X'p),  an  order  Hankel  function  of  the 
m  ' 

first  kind 

.  an  order  Bessel  function; 
m'  ’ 


and  e,o  refers  to  even,  odd  mode,  respectively. 


Evidently,  it  is  permissible  to  assign 


and,  correspondingly,  the  metric  parameters  are  given  by 


=  (1,P,1)  . 


To  be  more  specific,  it  is  assumed  that  an  arbitrary  electric  point 
source  of  strength  Pe  is  located  at  r*  »  e  S,  and  the 

observation  point  is  at  r  *  (p,t,2)  e  V+  or  V",  Then,  according  to 
Equation  (3,43),  Gq  has  the  specific  form  below: 

A  A 

=+--  .  PP_. 

Go(r,r')  =  U(p-p')g^(r,r*)  +  U{p'-p)g^(r,r' )  -  «(r-r')  . 

(3.48) 

It  follows  that  the  resulting  electric  field  is  given  by 

“i  —  si  —  »  —  —  i 

E“(r)  =  j<^UQgQ(r,r')  •  p^  ,  r  e  V  '  .  (3.49) 

-±.  =+___  -± 

Since  E  (r)  or  gQ(r,r')»p^  is  solenoidal  for  r  e  V  ,  they  can  be 

expanded  in  terms  of  the  solenoidal  eigenmodes  as  follows: 


=  /  dh' 
^  *00 


*  rq(h',r') 
">=0  Cq(h',r') 


M-,.(h-,r)  + 


Tq(h’,r') 
,Dq(h'  ,r') 


N-,.(h',r) 

(3.50) 


where  and  denote  the  expansion  coefficients;  q  is  a 

shorhand  notation  for  and 
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^  0  0 


•  "’^O  •  (3.51 


It  may  be  noted  that  (3.51)  will  also  hold  if  M“j^,(h',r)  is  replaced  by 

“±  -  "±  -  “±  - 
since  both 

The  determination  of  the  expansion  coefficients  in  (3.50)  requires 
a  set  of  orthogonal  properties  of  the  vector  wave  functions  which  are 
stated  in  the  form  of  lemmas  below.  The  proofs  of  these  lemmas  can  be 
found  in  Appendix  D. 

Lemma  1. 


/  dh'  /  dz  /  Pd»  lAe„  ;  X  7  X  fi^Sx(d,f))-[re„M5;;x.(h',?)J 
0  0  0  0  0 

—  Sn  3  Sn, 

"  ^n  V  3p  ^m  *mn  • 


Lemma  2. 


/  dh'  /  dz  /  pd<^  [Ae,  p  X  7  X  Ne;;x(h.F)j.[re^  Cx' •") ] 


•  OO  .00  Q 


0  0 


^  3  S  Sm 

7“  S  Tex,  ^  Z  "(ip)Z;7(Xp)  5„„  . 


Lernma  3. 


-S  -S 

/  dh*  /  dz  /  Pd^  [Ae^  p  X  V  X  Me;;;^(h.F)].[re^  ^  • 

•  eo  .00  0  0  0  0  0 

Letnma  4. 

”  *  2tt  »  5  c 

/  dh'  /  dz  /  pd<^  [Ae^  p  X  V  X  °  • 

-o»  -oa  0  0  0  0  o' 

In  these  Lemmas,  the  shorthand  notation  Ag^  Pg  has  the  same  meaning 

0  0^ 

defined  by  (3.51);  e  -jj  i  fg^  f,=o  and  2  for  n>l;  is  the  Kronecker 
delta  function  which  is  1  when  m=n  and  0  otherwise;  and  stands  for 
•'+"  or  sign,  n=0,l,2 . 

It  is  worth  pointing  out  that  the  orthogonal  properties  of  the 
cylindrical  vector  wave  functions  have  been  established  in  a  mixed 
domain  of  (<|i,z,h)  which  is  partly  spatial  and  partly  spectral.  This 
implies  that  these  functions  are  orthogonal  on  any  closed  surface  at 
which  the  p  coordinate  is  constant.  It  is  the  very  feature  that  makes 
the  eigenfunction  expansion  process  simple  and  elegant.  It  should  be 
mentioned  that  different  forms  of  orthogonal  properties  of  these  vector 
wave  functions  have  previously  been  investigated  by  Stratton  [15]  who 
used  a  mixed  domain  of  (p,iti,X);  and  by  Tai  [19]  who  used  the  entire 
spatial  domain  of  (p,i^,z).  Since  coordinate  p  is  involved,  both  [15] 
and  [19]  spent  extra  effort  in  dealing  with  the  orthogonal  properties 
of  the  cylinder  functions. 


To  see  how  the  expansion  coefficients  (Ce^,  Tg^  and  Tg^)  are 

0  0  o"  0^ 

determined,  it  is  helpful  to  consider  a  region  V  containing  the 
point  source  J(r)  =  p  6(r-r')  and  bounded  by  a  surface  in  denoted  by 


=  {r  =  (p^,<|i,z):  p^>p'  ,  -"•<z<«l  , 


and  another  surface  in  V"  denoted  by 


S~  =  {r  =  (p“,(|),z):  p"<p'  ,  0<(j><2n^  _a><z<»} 

The  geometry  is  shown  in  Figure  3.4,  Now  Green's  theorem  (vector 
version)  is  applied  to  the  pair  (Mg^^j^Ch.r) ,  G  (r,r')»p  )  over  region  V: 


I  [Menx(^*f')*'^><'7x(Vr,r')«p  )  -  VxVxMe^j^(h,r)  *(6  (r,r' ) ‘p Jjdv 
Vo  ®  0  ® 


=  I  [(GQ(r,r')«Pg)xVxMe^j^(h,r)-Me„;^(h,r)x7x(Gp(r,r')»Pg)]-nds  . 


S'*’US" 


(3.52 


Since  GQ(r,r')  and  Mgf,j^(h,r)  are  solutions  of  (3.11)  and  (3.29), 

0 

respectively,  the  L.H.S.  of  (3.52)  can  be  reduced  significantly  as 
f ol 1 ows : 


L.H.S.  =«  / 


V  0 


=  /  Me„;^(h,r)*Pg6(r-r')dv  =  Me„j^(h,r')-pg 


From  (3.48)  and  the  fact  that  n=p  at  S  and  n=-p"  at  S",  it  is  clear 


R.H.S.  =  /  [(gQ(r,r')‘P  )x7xMe„x(h.r) 

S+  0 


-  f^nx(^»^)x^’'(9o(f'.'*')*Pp)]*P^ds 
0  u  e 


+  /  [(go(f'.r')-pg]x7xMenx(*i»r) 


■  Menx(f’.r)x7x[g"{r,r')*p  )]*(-p")ds 


Employing  the  expanded  form  of  (9o»pg)  in  (3.50),  the  two  surface 
integrals  of  (3.54)  can  be  manipulated  as  follows: 


/  [  ]»P  ds 


Cn  .  • 

/  dz  /  p  d«fr  /  dh'  I 


+  Te^(h'.r')Nen,x'(^’'.^)]x7xMe  j^(h.r)*p'^ 
0  0  0 


=  -  I  /  dh'  /  dz  /  p'^d^  P^x7xMe^j^(h.r) 


m=o  - 


0  0  0  0 


-  [ren,{h',F')7xMj^;^.(h'.?)  +  Te^(h' ,F' )7xNe^;^.  (h* .?)  ] 


•  Menx(H.f') 


Using  Lemmas  1,  3  and  4  one  obtains 


/[>P*<I5  =  -  I  |^re™pVHy>(Xp*) 

5+  m=o  n  0 

lii  ■•• ,  (1)  ♦  ’  (1)  ♦  ~ 

-  p  W  (xp  )^hJ  '(Xp  )S^  .  0 


/  [  ]*(-p~)ds 
S' 

•  2ti  0.  «o 
=  /  dz  /  p"dt  /  dh'  I 


m=o 


[C, 


0  0 


+  Den,(h'  ,r' )Ne„j^.  (h' ,r)  ]xVxMe^;^(h,r) -(-p") 
0  0  0 


Menx(h.r)xVx[CeJh',?‘)M;^;,.{h'.F)  +  Dejh' .F*  (h*  ,F)  ] -{-p") 

o  0  0  0  0 


2x 


I  /  dh'  /  dz  /  p‘d^ 

cn— 0  "o#  *00  0 


p"x7xMe„;^(h.r) 

0 


0  0  0  0 


-  p'x7x  [Cen,(h'.f')MemX'^^'»^^  ^  ^^emx' i 


•  Menx(^r) 


Again,  making  use  of  Lemmas  1,  3  and  4, 


/  [  ]*(-p")ds  *  I 

j+  I11“0 


4it^  -  -«(1)  -  S  • 

—  Cem(-^r*)p  x2h;"'(Xp  )— JJXp  )6„ 


n  0 


3p-  '''mr 


4x^ 


-  -  -  9  (li  - 

Ce„(-h.r')»  X20„(X.  )^h;,  '(Xp  )S„ 
mo 


4x2  _  _  _,  _  -(IV 

7~Cen(-h.r')p  x2.x[h;'^Xp  )J^(Xp  )-J^(Xp  )h;;  (Xp  )] 


The  Wronskian 


WlH<'>(z),0„(z)l  .  h'‘>(z)0>)  -  H<‘>'(z)J„(z)  »  -  ^  . 

Hence, 

-jSir 

/  [  ]  •  (-P  )ds  *  -j—  Cen(-h.r') 

S-  no 

It  is  obvious  from  (3.55)  and  (3,56)  that 

-jSit 

R.H.S.  =  /  [  ]  •  p  ds  +  /  C  ]  •  (-P  )ds  =  — —  x2  Cej,(-h,r') 

S+  S-  no 

Equating  L.H.S.  and  R.H.S.  gives 

-+  -  _  -jSirX^ 

Menx(^.f'‘)*Pe  *  ~  Cen(-h.r') 

0  n  0 

Reversing  the  sign  of  h  and  solving  for  Ce^  gives 

0 

j  ^n  -+  __ 

Cen(^''  )  =  8^Menx(-h.r*).p  . 

0  0 

If  one  applies  Green's  theorem  to  the  following  pairs: 
(NJnx(^’.0»  S  (r,r')*p  ),  (Menx(^»^)»  ^ *Pe^ *  (Nenx(^»^^» 

)  *Pe^  •  follows  the  preceeding  steps  closely,  then  one 

obtains 


VJ  V"' 


_  J  “+  -  > 

“jn*’’-'''*  "  8^x2  "enx(-''''''l-Pe 


J  -_  _  - 

ren(h.r')  =  Menxt-^.f*' ) ‘Pe 


(3.58) 


(3.59) 


a^n  -  - 

Ten(^»^')  =  3^x2  ) *?© 


(3.60) 


respectively. 


Substituting  (3.57)  -  (3.60)  into  (3.50)  results  in  the  complete 
eigenfunction  expansions  of  the  solenoidal  components  of  Gq  given  below: 


/.'‘''.Li?  ^(-h.f) 


.«  n®o 


+  N|nx(^’»^)Ne-x(-h.r')  . 

0  0  (3.61) 


Hence,  the  formal  eigenfunction  expansion  of  Gg  has  been  completed.  It 
is,  however,  of  interest  to  mention  that  Tai  [19]  has  shown  that 


»  for  all  r,r'  e  Vjj  ,  (3,62) 


where 


(3.63) 


9n(r.r')  = 


jkolr-r' 


4ii|r-r' 


Since,  for  r^F',  ^^(r.r')  Is  identically  equal  to  either  iQ(r,r')  or 
gQ(F,r')  (depending  on  the  location  of  r),  it  is  straightforward  to 
conclude  from  Equations  (3.61)  and  (3.62)  that 


”  ^n  -±  _  -:f  _  -♦  _  -X 

8ir  /  I  \2  Menx(h.r)Menx(-h.r')  +  Nenx{h.r)Nenx(-^’.'*') 
-•  n=o  0  0  0  0 


(3.64) 


=  (I  +  1(2  77)  g  (r,r')  .  r  ^  r* 

0 


0.  GREEN'S  DYADIC  FOR  A  DIELECTRIC  COATED  CONDUCTING  CYLINDER 


The  geometry  of  a  dielectric  coated  cylinder  is  illustrated  in 
Figure  3.5.  It  consists  of  an  infinitely  long,  perfectly  conducting 
circular  cylinder  of  radius  a  coated  by  a  uniform  dielectric  layer  of 
thickness  t  =  b-a,  permittivity  =  epeo  and  permeability  m  =  uq* 
is  in  general  complex  to  account  for  lossy  dielectric.  The  z-axis  of 
the  cylindrical  coordinate  system  (p,t,z)  coincides  with  the  axis  of  the 
cylinder. 

In  order  to  determine  the  electric  Green's  dyadic  for  this 
geometry,  it  is  convenient  to  divide  the  entire  space  surrounding  the 
conducting  cylinder,  denoted  by  V,  into  two  distinct  regions  defined 
below: 


Figure  3.5.  Geometry  of  a  dielectric  coated  cylinder 


Vj  »  {p  ■  (p,(|»,z):  a<p<b,  0<^<2ir,  -»<z<-}  , 


(3.65) 


^^2  “  t**  “  (p«^»z):  b<p<«,  0<i><2n,  -«kz<«} 


(3.66) 


As  a  consequence,  Vj  1s  the  dielectric  layer,  while  V2  is  the  air  region 
surrounding  the  coated  cylinder,  and  V  =  V1UV2. 

In  analyzing  conformal  microstrip  antennas.  It  Is  sufficient  to 
consider  the  electric  current  source  that  Is  located  In  region  V2  only. 

B 

Consequently  the  electric  Green's  dyadic  Gg  for  region  V  will  consist 

s  s 

of  two  components,  denoted  by  Gel  ^^d  Ge2'»  such  that 


Gg(?,?')  =  U(p-b)Gg2(^»^')  +  U(b-p)Ggj(?.?') 


(3.-7) 


where  Gel  and  Ge2  satisfy  the  following  vector  wave  equations: 


VxVxGg^(r,r')  -  ki(5gj(r.?’)  *  0  ,  ?  e  Vj,  F'  e  V2  ,  (3.68) 


Vx7xGg2(r,r')  -  k2Gg2(r,r')  =  i6(r-r')  ,  r,r'  e 


(3.69) 


2  2  2  “  “ 
kj  =  and  k2  »  <a^V2^2  *  “^^0^0  “  *^o*  general  terms,  ^^^(6^2) 

can  be  thought  as  the  electric  field  In  Vi(V2)  due  to  a  point  source  In 

V2.  Correspondingly,  there  exists  two  mangetic  Green's  dyadics,  denoted 

B  S 

by  Gmi  and  Gni2,  such  that 


7xVxGi^l(r,r‘)  -  kiGn,i(r.r')  »  0  .  r  e  V^,  r'  c  ^1^  .  (3.70) 

*  7xi6(r-r')  .  r.r'  e  ,  (3.71) 

and  they  are  related  to  their  electric  counterparts  as 

^  “  11^  )  ,  r  e  V^,  r'  e  ,  (3,72) 

Sn,2(r.r')  =  VxGe2(r.r')  ,  r.r'  e  V2  .  (3.73) 


Equation  (3.73)  follows  directly  from  (3.9),  while  (3,72)  is  established 
In  Appendix  E  which  also  shows  that 

Ei(r)  =  juug  /  Ggj(r,r')  •  J(r')dv'  ,  f  e  (3.74) 

*2 

“  I  •  j(r')dv'  ,  r  e  Vj  (3.75) 

''2 

where  Ej  and  Hj  refer  to  the  electric  and  magnetic  fields  in  Vj  due  to 
source  J  in  V2.  And,  of  course,  it  follows  from  (3.17)  and  (3.18)  that 

E2(r)  »  j‘^U2  /  G  2(r,r')  •  J(r')dv’  ,  r  e 


(3.76) 


^2^^^  “  ^  •  j(r')dv'  ,  r  e  .  (3.77) 

''2 

where  and  represent  the  electric  and  magnetic  fields  in  V^. 
may  be  noted  that  the  magnetic  Green's  dyadic  for  region  V  can  also  be 
cast  into  a  form  similar  to  (3.67)  as  follows: 

Gn,(^r')  =  U(p-b)G^2(^*^‘)  +  U(b-p)G^l(?.r' ) 

=  U(p-b)7xGg2(r.r')  +  U(b-p)—  vxGgj(r,r' )  .  r'  t  .  (3.78) 

However,  G^  will  not  be  involved  in  the  following  analysis,  and  it  is 
being  intr'^duced  for  the  sake  of  completeness. 

s  s  s 

To  specify  Gg  uniquely  for  region  V,  it  requires  .Ggj  and  Ge2  to 
satisfy  the  boundary  conditions  pertaining  to  V.  These  conditions  can 
be  readily  derived  by  considering  the  corresponding  boundary  conditions 
imposed  upon  the  electromagnetic  fields  in  V.  First  it  is  observed  that 
the  tangential  electric  field  vanishes  at  the  conducting  surface  of  the 
cylinder,  namely, 

*  w 

p  X  Ej(F)  -  0  .  (3.79) 

p=a 

Secondly,  at  the  interface  between  regions  Vj  and  Vg,  the  tangential 
electric  and  magnetic  fields  are  continuous,  i.e. 


83 


p  X  [E2(F)  -  =  0 


P  X  [H2(F)  -  H^(F)]p^jj  =  0 

Next  substituting  the  integral  representations  of  the  electric  and 
magnetic  fields  defined  by  (3.74)  -(3.77)  Into  (3.79)  -(3.81)  one 
obtains  in  the  following  set  of  boundary  conditions  on  the  dyadic 
fields: 

px6gj(F,F')  =0 

p=a 

P  «  [''>'Se2(P.P')  -  ~  VxGjj(?,?')]„.i,  =  0  . 

In  addition,  Ge2  must  also  satisfy  the  radiation  condition: 

Aim  r  [Vx6g2(r,r')  -  jk2rxGg2(r,r' ) ]  =  0 

The  simplest  way  of  finding  Ge2  is  to  apply  the  principle  of 
scattering  superposition.  Thus  one  will  consider 

Ge2(r.r')  =  ^  G  e  . 


i 


In  Equation  (3.86),  Gg  is  the  free  space  Green's  dyadic  found  in  Section 
C;  and  is  expressed  explicitly  below  for  convenience: 


GgCr.r')  =  U(p-p')g^(r,r')  +  U(p'-p)g^(r,r' )  -  6(r-r') 


=  U( p-p' ) 


“  c. 


n 


8ir  /  dh  I  ^2 

n=0  0  0 


0  0 


1  J  *  *  Cp 

I  -«  n=0  0  0 


i 


+  Nenx(^’.^)Nenx(-^’*^')] 
0  0 


PP 


■  S(r-r')  ,  r,r'  c  V^, 


(3.87) 


and  Gs2  satisfies  the  homogeneous  vector  wave  equation  and  therefore 
must  have  the  form 


<"  Cl 


n 


8s2('-.''')  ■  8.  L'"’  „l„  1? 


lAe„Se„x(h,f)  *  (h.F) ) 

0  0  e  e 


”enx(-f’.''' )  +  [CenNenX^^’*'*)  +  ^  Nenx(-h.r') 


0  0 


^2  •  (3.88) 
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UNCLRS5IFIE0 


NL 


It  is  clear  that  the  sum  of  and  '♦ill  automatically  satisfy  (3.69) 
that  governs  8^2.  Loosely  speaking,  6^  may  be  interpreted  as  the 
incident  or  direct  field,  and  G^2  as  the  scattered  field.  As  for  the 
form  of  Gg2.  the  choice  of  anterior 

elements  is  dictated  by  the  radiation  condition  (3.85),  and  the  choice 

of  Menx(-*'»f*' )  and  Ne^;^(-h,r')  as  the  posterior  elements  is  guided 
0  0  _ 

by  the  expression  for  g“  because  various  boundary  conditions  such  as 
(3.83)  and  (3.84)  can  be  satisfied  only  if  the  posterior  elements  of 

s 

and  G^2  are  the  same. 

Once  the  form  of  Gs2  is  determined,  it  is  relatively  easy  to 

a 

formulate  the  expansion  for  Ggi  given  below: 


“iT  /  ■"'.L  $ 


n«0 


0  0 


^enMenM(^».r)  +  TonNony(h.r)  ] 

0  0  e  e 


”enx(-H.r')  +  [aenNenn(h.r)  +  bon^  (h,f) 

0  0  0  e 


^  ‘‘enN;„^(h,F)  +  foX^(h,r-)] 

A  A  A  A  A  I 


toil 


In  the  case  of  6^^,  the  posterior  elements  Me„;^(-h,r')  and  Nej,j^(-h,r' ) 

B  0  0 

must  be  chosen  to  match  those  of  6^2  such  that  the  boundary  conditions 

at  the  interface  (p=b)  can  be  satisfied,  and  the  anterior  elements 
-  -±  - 

Meny(li»f*)  and  (h,r)  are  chosen  because  they  are  the  solenoidal 
eigenmodes  and  will  constitute  the  bouncing  waves  inside  the  dielectric 
(region  V^), 

It  is  important  to  note  that  the  parameters  X  and  u  which  appear  in 
(3.87)  -  (3.89)  are  defined  as  follows: 


‘V 


-  h2  ,  Im  X  >  0  ,  Re  X  >  0 


u  kf  -  h2  ,  Im  u  >  0  ,  Re  u  >  0 


(3.90) 


(3.91) 


There  are  totally  24  unknown  coefficients  in  the  expansions  of  Ggl 
and  Ge2.  Not  too  surprisingly,  however,  there  are  exactly  24  linear 
equations  that  govern  these  coefficients  after  expanding  the  boundary 
conditions  (3.82)  -  (3.84).  This  24  x  24  system  of  simultaneous  linear 
equations  can  be  solved  analytically.  The  algebraic  details  are 
presented  in  Apendix  F  and  only  the  final  results  are  listed  here: 


^en  =  A  (3.92) 

0 

®en  “  (3.93) 
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21 

^  ^kfHl  -  l2  )  H„(Xb)Y, 


Jn(Xb)  vi2  Tfj 

■  H„(xb)  M  x2  h;^ 


^ki)  x2  H„(xb) 


Jn( k2  Sq 

■  H„(Xb)  +  UJ  X2  Hn(Xb) 


( 


and 


■  i;  L Jo 

*  'n  (Nonafi.r)  ♦  P„  "o„„(l'.r))] 

Menx(-h.r')  +  [d^  (Neny(h,r)  +  Ng^yC^.r)) 

"  ^(Mon,(h.r)  ^  q„  M;„^(h.r-))]  Knx^-K^')  . 

(3.121) 

This  completes  the  construction  of  the  dyadic  Green's  function  for  the 
dielectric  coated  cylinder, 

E.  EXPANDING  THE  GREEN'S  DYAD ICS 

One  may  have  already  observed  that  dyadic  functions  are 
notationally  so  compact  that  they  greatly  simplify  many  formal 
manipulations.  Ironically  it  is  this  very  compactness  feature  that 
tends  to  conceal  Information.  Often  enough  deeper  insight  can  only  be 
gained  after  the  dyadic  functions  have  been  expanded.  Anticipating  that 
the  Green's  dyadic  for  the  coated  cylinder  will  be  used  in  the  analysis 

s 

in  the  chapters  that  follow,  it  may  be  worthwhile  to  expand  Ge  into 
component  form  such  that  the  electric  field  can  be  easily  read  off  once 
excitation  is  known. 


Various  components  of  Gg  in  regions  Vj  and  V2  are  presented  in  the 
following. 


Region  V2  (air):  b<p,p'<-,  0  -«kz,2'<» 


It  follows  from  Equations  (3,67),  (3.86)  and  (3.87)  that 


,  =+ _  _  pp  = 

=  U(P-P  )go(r.r')  +  U(p'-p)g^(r,r')  -  «(r-r')  +  G^2(F,F'). 


O'-  »■  '  k? 

Making  use  of  (3.46)  and  (3,47)  in  (3,87)  yields 


(3.122) 


g;;(r,r')  = 


j  «  “  en 

8n  /  I  ^ 

-«  n®0 


PP 


-n2 


hX  ^  ■ 

')  +  fr")  z  ixp)z  fxp* 


^Z„(Xp)Z„(X,’)  *  (-)  Z„  (Xp)Z„  (Xp') 


•  cosn(  )  +  p(>' 


nX  + 


nXh2 


;z„(xp)z„  (Xp-)  *  J|7Z„  (Xp)Z„(Xp') 


sinn(4-^') 


A  A 

jhx3  - 

A  A 

+  PZ 

(Xp)Z„(Xp') 

cosn( ♦-<>')  -  (^p* 

nx  +•  5 

;rz-  (xp)z„(xp') 


nxh^ 


^  k^p 

nh  ^  1 

*  (IT 

4 

PP 

A  A 

“jhx3 

-zp' 

1 

1^2 

A  A  1 

+  zz 

sinn(  (>-♦')  + 

A  A 

cosn((fr-4')  -  4z 


X2Z*'(Xp)Z^'(Xp') 


^nX^h  t  + 
k?p 


cosn(^-^')  - 


jnhx^  +  j 
‘Wz-(xp)z„(xp’) 


sinn(^-^' ) 

sinn( ) 


’X- 


1(2  Z^(Xp)Zj^(Xp')  cosn((^-^')  e 


(1.123) 


where. 


^n(0  ■  Hy’(C):  Z„'(6)  •  0„(C)  and  Z^'(t)  -  ^Z*(e) 


It  Is  of  interest  to  point  out  that  formally  carrying  out  the 
=  1 

operation  (I  +  ^2  vv)  g.Cr.r')  in  (3.64)  will  lead  to  the  closed  forms 


s  X 

of  g~  given  below: 


V''*'''*  =  4^1 


jk2R  jk2R 


pp'  |_(k2 


+  P  1^3  )  COS(^-*') 

jk2R 


+  Q(p-p‘cos(<|»-<t.'))  •  (pcos(<fr-t')-p*) 


R3  J 


+  P(P' 


0k2R  jk2R 

e  e  “I 

^  k^ — ^ — +  [p+p(p-p'cos(^-(j>'))Q]  • -p  jsin(t-^') 


jk2R 


+  PZ 


|_(z-z')(p-p'cos{i>-t‘))Q  “p — _ 

_  jk2R  jk2R 

"  L  *^2  *  “rS — _|sin((|.-(i>') 

jk2R 

lee  < 


jk2R  jk2R 
e  e  e  “ 

(k^  g  +  P  p“  )  cos(t-^')  +  pp'Qsin2(t-i>' )  p“_ 


+  ♦z 


2  R 

jk2R  jk2R 

^  6  ^  AA  I**  0  ™ 

_p'(z-z')s1n(<|i-(ti')Q— p-_  +  zp'|U-z')[pcos((|)-<i>')-p')0— p“ 


jk2R 


+  z^' 


jk2R  jk2R 

e  e 


_p(z-z')sin((»-(»')Q  p"_  +  zz  |_k2  y-  +  [P+(z-z')^Q)  p" 


(3.12 
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It  o 


|r-r'|  =  [p2  +  p'2  -  2pp'cos((|.-4>’)  +  {z-z')^] 
jkaR-l  , 


Q  = 


33^2  3 


Next  substituting  (3.46)  and  (3.47)  into  (3.120)  gives 


J 


-«  n=ij 


pp 


n2 


jnh  1  '  1  * 

♦  ®r,  [7H„(Xo)H„(i»')  t  7h„(1p)H„(Xp')] 


h  =>  . 


cosn(^-i>' ) 


+  P<>' 


*  ®n  kT  lH„(MH„(Xp') 


nh^ 


^  X2pp'  H,(Xp)H^(Xp')]  +  Cj;;^  H„(Xp)H^(Xp')] 


sinn((fr-(j>') 


+  PZ 


rJhX 


‘^nlTT  '  ®n  1717 


cosn(4»-(^' ) 


•♦p' 


n  I 


MT7""pl^‘''"n<*‘’'>J  *  ®n  ir  17^'  H„(Xp)H„(Xp') 


nh^ 


♦  H„(Xp)H„(Xp')]  *C„  [7|I7H„(Xp)H„(Xp’)] 


sinn(  ) 


I  I 


*  1  a„[h„(xp)h„(xp’)]  *  b„  ^7  I  7  h„(xp)h‘{xp')  ♦ 


*  p'  H„(»<>)H„(Xp')l  ♦  ^H„(Xp)H„(lp')l  cosn(».t') 

'**  inh  X  • 

,  n  jhX  I 

^  ®nt(i^')Hn(Xp)H^(Xp‘)]  -  H^( Xp)H^( Xp' ) ]  cosn(*-*') 

^  '  ,  jnh 

*  "♦  '  <:n[kI7 

**  ^  ^  ihf 2-7M 

+  “  CJ(— )  H^(Xp)H^(Xp')]  •cosn(*.^')  , 

-  _  (3.12 

where 

H„(5)  -  h‘‘’(?)  and  H^(5) 

Region  Vj  (dielectric):  a<p<b,  b<p'<<»»,  <2tj,  -«»<z,z'<<»» 

It  is  clear  from  Equation  (3.67)  that 

K^r,r)  =  6gj(r,r') 

Using  Equations  (3.46)  and  (3.47)  In  (3.121)  yields 


,  jnhu  •  “  I 

^  ^n  UiX^p'  [Jn(wp)+PnH^(up)] 

h^M  '  <  I  jnh  i 

"  "nlkjkji  ”„(»p')1[0„(upKp„H^(,.p)1  +  [—  H„(Xp')] 

[Jn(up)+q„Hn(up)]  cosn(k-*')  ♦  p«'  t„  [i;^  H^{ Xp' ) ](j„(up)tq„H^(up)) 

rilm  '  ,  '  •  ,  nh2u 

^n  [kiX  Hn(^p')]  [Jn(up)+PnH^(up)]  +  H^(Xp')] 

r,'  '  ,  jn^h  “I 

Nn(pp)+Pn^n(^P^J  +  [k2X2pp'  Hp^  ^  J  ^  ^ 


,  jnjj  •  •  n 

«  “ali^  H„(X»')1[0„(u»)*p„H„(p»)]  -  f„  [—  H„(Xp’)] 


r  1  1  ny 

LJn{np)+qnH^(iip)]  cosn((^-*')  -  (>p’  l'j^Hp(^p')] 

r,'  '  jhn2 

lJ„(»PKq„H„(up)]  +  t„  [|<,j2„'  H„(Xp')][jJuo)kp„H„{up)] 

nh2  I  , 

*  Mk.k^Xp  H„(Xp')l  [0„(ppKp„H„(„»))  k  f„  [i^H„(Xp')] 

•  lJ„(»P)kq„H„(up)]  s1nr(p-*-)  *  H^( Xp' ) ][J^(up)+q^H^(up) ] 

■jnh  >  p2h2 

'p  [k,Xp  H„(Xp'))[J„(pp)*P„H„(pp))  *  <l„[k,k,x2»p'  «n('»'>l  • 


Kl 


r  jnnu  I  I 

*  tk2X2p'  ^p' )  ] [J„{wp)+q„H^(wp) ]  jCOsn(*-(fr' ) 


I  r  jnri  \x 


r,  '  nu^ 

LJn(iip)+qnH^{pp)]  sinn(^-.fr')  +  zp'  UiX^p' 

r  jhu^  *  I 


'  jP  hy2 

*  'n  [m  H„(x»'1[J„(,o)»p„H„(„p)]  -  d„  (t;,^  H„(Xp’)] 

•  [Jn(pp)+PnHn(pp)]|sinn(t-<t.')  +  zz  i  H^( Xp' ) ] 


•  [Jn(up)+p^H^(pp)] jcosn(t-(^')  i  ^  ^ 


(3.126) 


where 


H„(5) 


Hy’(5)  «nd  H>) 


Thus,  the  dyadic  Green's  function  for  the  coated  cylinder  has  been 
expanded  into  a  form  from  which  the  electric  field  components  can  be 
readily  identified  once  the  current  density  is  specified. 


97 


CHAPTER  IV 


PATCH  DIPOLES  ON  A  DIELECTRIC  COATED  CYLINDER 

A.  INTRODUCTION 

The  main  goal  of  this  chapter  Is  to  obtain  expressions  which  will 
enable  one  to  compute  either  directly  or  Indirectly  the  self  and  mutual 
impedances  between  expansion  dipole  modes^  and  to  evaluate  the  elements 
of  the  voltage  vector  for  the  case  of  a  dielectric  coated  cylinder. 

These  impedance  and  voltage  elements  are  crucial  to  implementing  the 
method  of  moments  to  analyze  the  coupling  between  microstrip  antennas. 
The  single  most  important  quantity  in  these  computations  is  the  electric 
field  excited  by  an  expansion  dipole  mode  at  the  air-dielectric 
interface  associated  with  the  dielectric  coated  conducting  cylinder. 

This  is  due  to  the  fact  that  an  appropriate  combination  of  expansion 
modes  is  used  in  the  moment  method  solution  to  model  or  approximate  the 
true  current  on  a  printed  antenna  patch. 

Extensive  use  is  made  of  the  electric  Green's  dyadic  constructed 
in  Chapter  III  to  obtain  expressions  for  the  electric  field  or  its 
components.  To  ensure  continuity  and  consistency  relevant  symbols 


or  notations  will  retain  their  original  meanings  defined  in  Chapter  III 
unless  stated  otherwise.  For  convenience,  some  pertinent  notations  are 
summarized  below. 

V  refers  to  the  space  surrounding  the  circular  conducting  cylinder 
of  radius  a.  It  consists  of  two  distinct  regions  Vj  and  V2  such  that 
V  =  V1UV2,  where 

Vj  =  (f  *  (p.i^.z):  a<p<b  =  t+a  ,  0<^<2ir  , 

represents  the  uniform  layer  of  dielectric  coating  of  thickness  t,  and 
V2  =  (r  =  (p.<l».z):  b<p<-  .  0<«<2it  ,  -oKzc-} 

stands  for  the  (free)  space  or  air  surrounding  the  coated  cylinder.  Th 
constitutive  parameters  for  Vj  and  V2  are,  respectively,  (m,ei)  and 
(u2.e2)  corresponding  propagation  constants  ki^u/ii'iVi  and 

_  m 

k2=(j)/iJ2e2«  The  Green's  dyadic  for  V  is  denoted  by  Ge(F,r'),  where 
r  e  V  designates  the  field  point  and  r*  e  V2  refers  to  the  source  point 
which  is  restricted  to  V2.  The  Green's  dyadic  can  be  conveniently 
represented  as 

Ge(r,r')  =  U(b-p)Ggj(r,r' )  +  U(p-b)Gg2(r,r' )  ,  r  c  V,  r*  e  V2  , 

(4.1) 

with  Gel  (Equation  (3.126))  being  the  Green's  dyadic  for  Vi,  and  Ge2 
(Equation  (3.122))  being  the  Green's  dyadic  for  V2.  Likewise,  the 
electric  field  is  written  as 


E(f)  =  U(b-p)E^(F)  +  U(p-b)E2(r)  .  ?  e  V 


(4.2) 


where  E^(r)  =  pEjp(r)  +  ♦Ej^(r)  +  zEj^Cr)  1s  the  electric  field  in 
given  by 

El(r)  =  j‘^U2  /  Ggj(r,r')  •  J(r')dv'  ,  r  e  .  (4.3) 

''2 

a  A  A  A 

and  E2(r)  =  PE2p^^^  electric  field  in 

given  by 

£2(0  =■  j 011*2  I  Gg2(r,r')  •  J(r')dv'  ,  r  e  V2  .  (4.4) 

^2 

J  denotes  the  electric  current  source  immersed  in  V2.  The  corresponding 
magnetic  field  H  takes  the  form 

H(r)  =  U(b-p)Hj{r)  +  U(p-b)H2(r),  r  e  V  ,  (4.5) 

with  Hj(r)  =  l/(jo)Pj)  VxEj(r),  r  e  Vj,  and  H2(r)  =•  l/(jo>U2^  ^xE2(r), 
r  e  V2.  Finally,  all  the  expansion  coefficients  associated  with  the 
Green's  dyadics  are  defined  by  Equation  (3.92)  through  Equation  (3.119) 
which  in  general  cannot  be  readily  expressed  because  of  their 
formidable  complexity. 
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It  should  be  mentioned  that  there  Is  an  attractive  alternative  to 
the  Green's  function  approach  for  calculating  fields  due  to  the  surface 
current  source  located  right  at  the  air-dielectric  interface  of  the 
dielectric  coated  cylinder.  Such  an  approach  is  discussed  in  Section  B. 
It  is  analogous  to  the  planar  dielectric  slab  case  treated  in  Chapter 
II.  Section  C  utilizes  the  Green's  dyadic  to  obtain  the  electric  field 
in  regions  Vj  and  V2.  The  excitation  is  taken  to  be  the  current  of  an 
expansion  dipole  mode  placed  at  the  interface,  A  general  expression  for 
the  mutual  impedance  between  two  arbitrary  expansion  modes  is  derived. 
The  computation  of  these  mutual  impedances  is  an  extremely  important 
step  in  the  moment  method  solution  to  the  problem  of  analyzing  the 
coupling  between  microstrip  antennas.  In  Section  D»  examples  of 
coupling  between  two  expansion  dipole  modes,  both  with  z-polarized 
currents,  are  studied.  It  reveals  certain  interesting  but  important 
numerical  aspects  of  the  Green's  function  solution.  In  the  last 
section,  an  expression  useful  for  far  field  calculation  is  derived  via 
the  method  of  steepest  descent. 


B.  A  SPECIALIZED  METHOD  OF  SOLUTION 


This  section  presents  a  method  particularly  suited  for  analyzing 
patch  antennas  printed  on  a  dielectric  coated  conducting  circular 
cylinder.  This  method  is  conceptually  simple  and  is  analogous  to  the 
approach  described  earlier  in  Chapter  II  for  the  planar  dielectric  slab 
case.  Thus  the  presentation  will  be  brief.  The  inclusion  of  this 


method  is  due  to  the  fact  that  It  does  provide  an  attractive  alternative 
to  the  Green's  function  method  on  which  the  analysis  in  later  sections 
is  based.  In  the  present  method,  the  electric  surface  current  on  the 
printed  conductor  that  is  situated  right  at  the  air-dielectric 
interface  is  introduced  only  through  the  boundary  conditions.  The 
solution  is  rigorous  since  the  effects  of  the  dielectric  coating  and  the 
cylindrical  conductor  are  both  accounted  for.  It  inclui^  the  radiation 
field  and  surface  waves.  Only  formal  expressions  for  the  total  fields 
(electric  and  magnetic)  are  derived  here. 

Being  a  surface  current,  the  microstrip  patch  current  is  introduced 
through  the  boundary  conditions.  As  a  result  Maxwell's  equations  have 
to  be  solved  for  homogeneous  source-free  regions  Vj  and  V2.  An 
arbitrary  field  in  ,  i=l,2,  can  be  constructed  from  two  scalar 
functions:  which  generates  a  TE  field  and  <|i^^  which  generates  a  TM 

field.  Both  \i)gi  and  satisfy  the  scalar  wave  equation: 


■  0  . 


where  ,  1=1,2. 


(4.6) 


The  electric  and  magnetic  fields  can  be  constructed  as 


follows  [16]: 


E.(0  =  pE4„(r)  +  .(r)  +  zE.Jr) 


(4.7) 


(4.8) 


Hi(r)  -  pH^p(r)  +  ♦H.^(r)  + 


The  explicit  forms  of  the  field  components  are  given  below  for 


convenience: 

Eip(r)  * 

13  1  32 

"  p  3(>  ’•'ei  "  juie^  3p3z  ♦mi  » 

(4.9) 

3  1  1  32 

Ei^(r)  = 

3p  '*’ei  "  p  3^32  'I'mi  • 

(4.10) 

Ei,(F)  . 

1  32  2 

"  ju»c^  Uz2  *  ♦mi  » 

(4.11) 

13  1  32 

Hip(r*)  = 

p  3(^  "^mi  “  j(i)u^  3p3z  ’^ei  * 

(4.12) 

3  1  32 

«i^(r)  * 

3p  ♦mi  ”  jidu^p  3^3z  ♦ei  * 

(4.13) 

1  32  2 

ju)u^  (322  *^i)  ♦ei 

(4.14) 

In  order  to  specify  the  fields  uniquely,  the  wave  Equation  (4.6)  must  be 
solved  together  with  the  following  boundary  conditions: 


on  the  surface  of  the  conducting  cylinder. 


p  X  Ej(r)  *  0 


(4.15) 


on  the  air-dielectric  interface. 


P  X  [E2(F)  -  ijCF)]  b  =  0  . 


(4.16) 


P  X  [H2(F)  -  Hi(F)]p=b  = 


(4.17) 


where 


(4.18) 


is  the  surface  current  density  on  the  interface  at  p»b.  Of  course,  the 
solution  must  be  chosen  so  that  the  radiation  condition  is  satisfied  as 
|r|  r  e  V^. 

Due  to  the  cylindrical  symmetry  of  the  structure  under 
consideration,  the  scalar  functions,  and  >lnni ,  can  be  expressed  as  a 
Fourier  series  -  Fourier  integral: 


1  -  - 


n=-c.  .«  i|;j^^{p,n,h) 


(4.19) 


where 


♦e^Cp.n.h) 


g-jn^  g-jhz 


dz  . 


(4 


Similarly,  Js{F)  can  be  formally  represented  as 


2Tt  I  {♦JcA(p.n.h)  +  zJ  (p,n.h)} 


dh 


» 


(4 


with 


2n 

J 


0 


c» 

/ 


Jsz(J') 


-jn*  -jhz 
e  e 


dz  d(^ 


(4 


Substituting  (4.19)  into  the  wave  Equation  (4.16)  leads  to  the 
Besssel's  equation  of  order  n: 


♦e^(p,n,h) 

?^^(p,n,h) 


=  0  .  i-1,2  , 


(4 


where  k  .  *  -  h^  .  The  solution  for  region  V,  can  be  written  as 


^el(p*'’*^) 

Ag(n.h) 

Hp'\k  p)  + 

Bg(n,h) 

\l{p,n,h) 

^m(n,h) 

n 

®m(n,h) 

(4.24) 


where  and  are  the  Hankel  functions  of  the  first  and 

second  kind  which  represent,  respectively,  the  outgoing  and  ingoing 
waves  in  V^.  As  for  region  ^2*  solution  takes  the  following  form: 


'I'e2(p*'^*^) 

^(P.n.h) 

= 

rg(n,h) 

^m(n,h) 

(4.25) 


since  there  are  only  outgoing  waves  in  V2. 

Making  use  of  (4.24),  (4.25)  and  (4.19)  in  (4.9)  to  (4.14)  one  can 
express  the  field  components  more  explicitly.  In  particular,  the 
transverse  field  components  with  respect  to  the  p-di recti  on  can  be 
rewritten  as  follows: 

in  region  Vj  (a<p<b,  0<(^<2n,  -<»<z<<>»). 


1 


nh 


-  Ill'  f?!' 

l  /  dh{k„  [aX  Mk  p)  +  '  (k  p)]+^ - 

19'  p-_oo  -00  ^  ^  ^  ^  e  n  '  pr jue^i 


(2) 


jn9  jhz 


(4.26) 
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=  2Tr  ^  J 

^  f^S.OD  mOO 


P2 


I-  «'"’e 


jni^  jhz 


■  2*  „i-  L'’''*'f'‘»2''">''i^’'"‘p2'’’! 


^  j«jy2P  ®  . 


*  2.  I  /.  dh 

fls^QO  «CO 


1- 


jojM^  ^e”n  ®  ® 


(1  r  d  f  i  ^ 

where  (0  =  — hJ^^(5),  1=1.2. 


The  unknown  coefficients  can  be  found  from  the  boundary 


(4.31 


(4.32 


(4.33 


conditions 


The  solution  procedure  is  rather  straightforward,  hence  the  algebraic 
details  will  be  omitted.  Only  the  key  steps  and  final  results  are 
presented  here.  From  boundary  condition  (4.15),  it  is  evident  that 


Next  it  follows  from  boundary  condition  (4.16)  that 

(1)'  (2)' 

Vp,t"n  "'pi'*')  ■ 

nl'  (11  (i>l 

*  ajh„  (k„,b)  -  vip  (Ap.M) 

(1)'  nh  (1) 

^e*^P2^n  *  jucjb  ^m^n  ^*'p2^^  * 

and 

Pi  fi\  (o\ 

-  IKp,*')  -  »»<  (Kp,‘’)l 


Finally,  boundary  condition  (4.17)  implies 


1 


(4.53) 


Ve 


zm  ~  j“Pi 


(4.54) 


(4.55) 


Thus  all  the  coefficients  have  been  completely  determined.  These 
in  turn  uniquely  specify  the  scalar  functions  >)»ei  and  ,  i=l,2,  from 
which  the  electric  and  magnetic  fields  in  both  regions  Vi  and  V2  are 
readily  obtained  via  Equations  (4.9)  to  (4.14).  One  may  have  already 
observed  that  this  solution  is  very  appealing  for  analyzing  problems 
that  involve  printed  dipoles  or  patch  antennas  provided  the  S'-.  "ace 
current  density  has  a  closed  form  transform  in  the  sense  of  (4.22).  For 
general  volumetric  electric  sources,  the  present  solution  scheme  will 
certainly  fail.  In  such  event,  the  Green's  function  approach  can  be 
used  although  finding  the  Green's  function  itself  often  presents  the 
most  difficulty;  however,  once  the  Green's  function  is  found  its  formal 
use  in  the  solution  is  straightforward. 


C.  MUTUAL  COUPLING  BETWEEN  EXPANSION  DIPOLE  MODES 


In  this  section,  the  field  evaluation  is  based  on  the  Green's 
dyadic  Gg  for  the  dielectric  coated  cylinder  constructed  in  Chapter  III, 
Gg  essentially  describes  a  linear  relation  between  the  electric  field 
and  an  elementary  current  source  as  indicated  by  (3.20).  Hence  the 
total  electric  field  E  produced  by  a  distributed  current  source  J  is 
obtained  through  an  integration  (summation)  process  described  by  (4,3) 
or  (4.4). 

In  the  analysis  of  mutual  coupling  between  conformal  microstrip 
antennas  on  the  dielectric  coated  cylinder  formulated  by  the  moment 
method,  the  chosen  set  of  expansion  dipole  modes  plays  a  very  important 
role.  The  efficiency  as  well  as  accuracy  of  the  analysis  depend 
intimately  on  how  closely  an  individual  mode  resembles  the  actual 
current  distribution  on  a  microstrip  patch,  and  on  the  transformability 
of  a  mode  in  the  Fourier  sense  (i.e,,  the  sense  of  Fourier  series  and 
Fourier  integral).  Furthermore,  one's  ability  to  accurately  evaluate 
the  field  due  to  a  mode  will  also  affect  the  quality  of  analysis  in  a 
grand  manner.  The  electric  field  obtained  via  the  Green's  dyadic  will 
be  an  exact  solution  since  the  dyadic  function  accounts  for  the 
conducting  cylindrical  ground  and  its  dielectric  coating  rigorously. 

A  microstrip  antenna  on  the  dielectric  coated  cylinder  consists  of 
a  printed  conductor  that  conforms  with  the  coated  cylinder.  The  printed 
conductor  is  typically  very  thin  so  that  the  current  flowing  on  it  can 


be  regarded  as  a  surface  current.  Hence  It  Is  appropriate  to  choose  a 
set  of  expansion  dipole  modes  as  described  below: 

H  »  {3^(r);  t  “  1,2,...,  L  <  •}  , 


where 


0  .  -r  t  S, 


(4.56) 


=  {r=(p.<>.z):  p=b‘'  ,  * 

The  symbol  b'*’  denotes  a  real  number  which  Is  Infinitesimally  greater 
than  b.  This  essentially  says  the  basis  set  s  consists  of  L  (finite) 
expansion  dipole  modes.  A  typical  mode  Jji  represents  a  surface  current 
density  which  vanishes  outside  a  rectangular  surface  patch  on  the 
alr-dlelectric  Interface  with  widths  In  the  ^-direction,  Zzi  In  the 
z-directlon,  and  centered  at  (b+,  tjjc,  z^c).  One  may  note  that  It  Is 
not  necessary  to  think  of  Sj^  as  being  a  conducting  surface  patch.  The 
reason  for  requiring  to  be  located  at  p=b+  Instead  of  p=b  Is 

consistent  with  the  use  of  the  Green's  dyadic  and  also  Is  closer  to  the 
actual  structure  being  modeled. 

Denoting  the  electric  field  In  V  due  to  mode  J)i,  t=l,2,...,L,  by 


if(?)  “  pEfp(f)  +  iEf^(?)  +  kf^(F)  ,  ?  e 


^*1.2  .  (4.57) 


V  , 
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and  employing  (4.3)  and  (3.126),  one  will  find 


--  n=0 

r  '  jnh  ,  hp  •  i 

•  '„Hn<“»l  ♦  r;!  1J„(|‘I>)+P„H„(UI>)1 

r  '  nh 

•  U'„H„(Xb) 

,  jhu  •  • 

n  "  «  / 2—2  *  ^ 

•  M  Kl""!  *  <l„Hn<“P>lW«'>)«'"(*-*'>Jtz('''>le 

( 

(l)U 

4  2  •  •  li  > 

-  Sir  I  ‘‘s'  /  I  sUx 
^4  -•  n=0 

r  •  jnh  nh 

•  l’'-."n(«>l  *  kjXb  Vn(«>)]  *7;i;  iJnlPb'^XnlPP') 

r  •  nh 

•  WVn(‘'>>  *  k;xbVn<‘'>> ))“’"(♦-♦' >V'> 

,  jnh 

'  -  jh(z-z') 

"  k2  *  qn‘‘n^‘‘P^]Vn^^P^)sinn(t-*')Jj^2(''')l® 


I- . 


^S,  “* 


jnh 


[Jim  ,  -  » 


2 

+  (  k^  l‘Jn(np)+PnH^(up)]d^Hj^(Xb))cosn(^-*')J^^(r')}e^^^^  ^  ^ 

(4.6 


and  similarly,  using  an  appropriate  form  of  Ge2  in  (4.4)  yields 


■  -  iT  I  “S'  /  «'>  i/nKi;;  H„(x»)[j^(xb) 

-•  n®0 

jnh  h  •  nh 

,  '  nh 

+  +]^C„H„(Xb))l)sinn(^.t')J^^(r’) 

JhX  . 

+  (]^H^(Xp)[j^(Xb)+C^H^(Xb)] 

n  _  jh{Z“Z') 

■  k  p  BpH  (Xp)H  (Xb))cosn(^-(f')J^^(r')}e  .  (4j 


.1 


'  jnh  nh  nh 


nh 


+  UB„H^(Xb)  +77^C„H^(Xb))])cosn(*-^*)J^^(r') 


,jnh 


■  te«n^^'^HJn(Xb)+C„H^(Xb)] 


(4. 


and 


X  • 

+  k7  Vn^  ^p)Hn(  ^  ♦-♦'  ) 

*  jh(z*z* ) 

^  ^kf  Hn(Xp)[J„(Xb)+CnH„(xb)])cosn(*-*')J^^(r')}e 

(4. 

#  J)  • 

It  should  be  remarked  that  the  expressions  for  £2^,  and  £2^  are 
derived  with  the  understanding  that  b*<p<».  Furthermore,  in  carrying 
out  the  calculation  of  the  electric  fields  using  the  Green's  function 
according  to  Equations  (4.3)  and  (4,4),  the  Jy(r')dv'  term  has  been 
replaced  by  the  J  (r')ds'  term,  where  and  J  denote  the  volumetric 
and  the  surface  current  densities,  respectively.  This  is  necessary 
because  the  expansion  modes  in  H  are  surface  currents.  Clearly, 


(4.64) 


Gg2(f‘»^' ) 


GQ(r,r')  +  =  9Q(r.r')  +  • 


where  and  G^2  expressed  by  Equations  (3.123)  and  (3.125), 
respectively. 

The  fundamental  definition  of  mutual  impedance  between  expansion 
modes  and  is  given  by 

~  ~  I  E2(r)  •  J^(r)ds  ,  t,  m=l,2,...,L,  (4.65) 

^m 

where  the  limiting  process  p*b*  is  implied.  It  follows  from  (4.56)  and 
(4.57)  that  can  be  rewritten  as 

hm  =  -  I  ^  •  (^*56) 

^m 

Letting  b+-*-b  and  defining 


~c 

Ja0(t,n';n,h) 


/ds  /ds'0,,(?M0„,(?) 


sinn( ) 


gjh(z-2') 


cosn( ) 


^mc+Zm  <l'mc+'t»m 

/  dz  /  bd4i 

Zmc'Zfii  '^mc“'^T\ 


Zbc'^Zj  , 

W  ^  bd*'J^Jb.t',z')J„g(b.4,z) 
^U~^l  ^IC~^Z 


sinn(  t-')'' 


)  gjh(z-z') 


cosn(^-^' ) 


(4.67] 


where  a  =  (t>,z  and  e  »  ♦.z,  one  will  arrive  at  a  general  expression 
for  Zjtni  hy  substituting  (4.62)  and  (4.63)  Into  (4.66)  as  follows: 


*  8.  1  '"'Xs  l(H„(»)lJ„(Xb) 

-*  n=u 

'  jnh  nh  nh 

*  *  k,Xb  ♦  i;ib  «„(“>)  li^ 

Jnh  X  • 

■  ^kfb  B^H^(Xb)H^(Xb))J^^(t.m;n.h) 

.  jnh  X  * 

+  l-  Ic2b  H„(Xb)[j„(Xb)+C„H„(Xb)]  +rB„H„(Xb)H^(Xb))J^2(t.n.;n,h) 


/■_  ~c 

+  (^2  H„(Xb)[j^(Xb)+C„H^(Xb)])J2^(t.m;n.h)l  .  4.m=U2....L. 


(4.68 


From  (4.68)  one  can  readily  obtain  the  mutual  Impedance  between  two 
♦-polarized  modes,  two  z-polarized  modes,  or  two  modes  each  with 
different  polarization. 

The  remaining  quantities  to  be  calculated  In  order  to  set  up  the 
moment  method  solution  are  the  elements  of  the  voltage  vector  which.  In 
the  present  case,  are  defined  as  follows: 


[f 


where  Ji  denotes  the  excitation  current  source.  Assuming  the  constant 
current  filament  model  described  in  Chapter  I  for  the  planar  slab  case 
still  holds  for  both  the  coaxial  feed  and  the  microstrip  line  feed,  each 
of  which  joins  a  microstrip  patch  at  (b.iis.Zs)  with  terminal  current  1^, 
one  can  immediately  write 


-  .  -  6(4'-<l>s)5(z-Zs^ 

JT(r)  =  pl^  - p  ,  r  e  Vj 


(4.70) 


Making  use  of  (4.58)  and  (4.70)  in  (4.69)  will  lead  to 

«  00  ^ 

"Mo 

•  !i(q„.o-i)  -  r(q„.r)| 

nh  ~5  jh 

^  k^Xb  ^kik2  ^n'^n^ 

1  ~c  jhZj 

-k  ^H^(^b)[l(q^,n.l)  -  r(qn,n)])j^(!l;4.5,n,h)}e  .t=l ,2. . . . ,L, 


(4.71) 


where 


I’(Pn,n) 

iMq^.n) 


b  <  Pn  ' 

I  [Jr,(up)  +  M^(up)]  dp 


LVI 


[Jf,(up)  +  ” 

-V  ^ 


(4.72) 


=  /  W  [Jn(up)  + 


H^(pp)]  dp  , 


(4.73 


JJt;(|),n,h) 

=  /  ds‘J^Jb.<>*,2') 

cosn(  ) 

St 

sinn((>-(j»' ) 

cosn((|.-i|.')  g-jhz' 
s1nn(<|>-(|>' ) 


(4.74 

with  a  =  <{1,2.  In  Equation  (4.71),  the  integral  I(qn,n-1)  does  not  have 
a  closed  form  and  hence  must  be  evaluated  by  numerical  means.  However, 
in  the  case  of  thin  coating,  where  t«a,  I(qn,n-1)  can  be  approximated 
as 


nq^,n-l)  =«  u(b-a)[J^_j(ua)  + 

u^(b-a)2  I  • 


(4.75 


which  is  obtained  by  integrating  the  two  leading  terms  of  the  Taylor 
series  of  the  integrand. 
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D.  NUMERICAL  EXAMPLES 

The  mutual  coupling  between  two  expansion  dipole  modes  is  studied 
in  this  section.  Figure  4.1  shows  two  expansion  dipole  modes  on  the 
surface  of  the  dielectric  coated  cylinder.  Mode  1,  centered  at  (b,0,0), 
has  width  2wi=2b  tii  and  length  2hi=2zi;  mode  2,  centered  at  (b,(|)2c.Z2c^ » 
has  width  2w2=2b'}i2  and  length  2h2=2z2.  Both  modes  have  only  z-polarized 
currents.  Past  experience  with  the  planar  dielectric  slab  problem  in 
Chapter  II  suggests  that  a  reasonable  choice  of  expansion  modes  (basis 
functions)  is  a  set  of  piecewise  sinusoids.  Hence  the  currents  of  modes 
1  and  2  are  defined  as 


^  2b(^  sin(<  z  )  *  ^  ^ 

0  ,  elsewhere  , 

(4.76) 


where  Sj^={r=(b,i»,z):  ,  \z-z^^\<.z^} ,  t=l,2,  and  (referring  to 


It  then  follows  from  Equation  (4.67)  that,  for  t=l ,  m=2. 


J^^(l,2;n,h)  =  ?^(1.2;n.h)  =  jJ^(1.2;n.h)  =  0 

J2z(1.2;n.h)  = 

dz  ^2cj‘*’2  Jl 

^2c'^2  ‘•’2c''*’2  '^1  ■'^1 

•  cosn(  ^ 


4k^k2 


^i4'2sin(<^z^)sin(K2Z2) 


sin(n<|i^)  sinCncji^) 
n  n 


> 


where 

cos(hz  )  -  cos(k  z  ) 
~  i  i  i 

=  {<2  -  h^l 


Substituting  (4,79)  and  (4.80)  into  (4.68),  with  t 

j^Z2c 

recognizing  that  the  integrand,  excluding  e  , 
of  h  yield 


•  (4.79) 

^2(b,<t)'  ,z'  )J22(b,<t>,z) 

■  cos(nt2c)Zi(h)Z2(h) 

(4.80) 

(4.8i: 

^1,  m=2,  and 
s  an  even  function 


J  dh  e„«(n.h)  -  - 


sin(n(t>i)  sin(n|i2) 


CJ 


•  cos(n*2c)  cos(hz2j  Zl(h)Z2(h) 


(4.82) 


where 


^^^2  <1  <2 


12  ■ 


(4.83) 


W(n.h)  =  x2H^(xb)[j^{Xb)  +  C^H^(Xb)] 


(4.84) 


Recalling  the  definition  of  Cf,  in  (3.110),  one  can  rewrite  W(n,h)  as 


"  ^ffb)  (k^)  [X2(n,h).X3(n,h)/Xi(n,h)J  * 


(4.85) 


where 


X  (n.h) 
1 


On  U2 

T„  *  ~ 


(4.86) 


X  (n.h) 
2 


^  Pn  k2  2  m2 

^  ^  ^k  J  A  '"n  ’ 


(4.87) 


X  (n.h) 

3 


,^2  nh  2  y  2  2 

[1  -  (7)  ]  . 


(4.88) 
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(4.89) 


r_  = 


Hn(Ab) 

Hn(Xb) 


= 


Jn(iib) 


J  (ua) 
n 

Hn(iia)  "n 


H„(ub) 


(4.90) 


J'(;ia) 

n 


(4.91) 


JpCwa) 


■  H  fwa) 


(4.92) 


Jr,(pa) 


and 


-  H*(iia) 


X  =  /k|-h^  ,  M  =  /k^-h^  »  with  Re 


*■  — 

_ 

X 

>  0»  Im 

1 

X  1 

1 

u 

Jj 

(4.93) 


>  0  . 


(4.94) 


It  may  be  assumed  that  Mj=U2  without  much  to  the  loss  of 
generality.  In  the  evaluation  of  (^*82),  there  are  several 

critical  points  that  need  special  attention,  namely: 


a.  It  is  computationally  more  efficient  to  perform  summation  over 
n  before  integration  over  the  h  variable.  This  is  due  to  the  fact  that 
the  sequence  of  cylinder  functions 

|J„(0,V?),h„(5).h„U))5 

can  be  generated  efficiently  for  5=ua,  ub,  and  Xb.  N  is  the  maximum 
index  number  which,  as  Indicated  in  Appendix  I,  is  dictated  by  the 
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dynamic  range  of  the  computer  being  used  to  compute  the  sequence.  The 
number  of  terms  summed  in  the  series  of  (4.82)  can  be  determined  by 
considering  the  well-known  series: 


(4.95) 


Supposing  6>0  is  an  error  threshold,  there  will  exist  an  integer  Nn,aix 
such  that 


N. 


S  -  I 


max  1 


n=l 


n 


2  <  6 


(4.96) 


It  is  clear  that  using  N^ax  as  the  terminating  index  value  for  the 
series  in  (4.82)  will  also  satisfy  the  same  error  criterion  in  the  sense 

m 

of  (4.96)  since  the  the  tail  end  of  the  series  is  dominated  by  I  l/n^ 

^ax 

In  general,  N  is  no  greater  than  N^ax.  The  evaluation  of  W(n,h), 
for  0<n<N,  is  straight-forward  provided  the  sequence 

has  been  calculated.  The  calculation  of  this  sequence  of  cylinder 
functions  is  discussed  in  Appendix  I  which  also  deals  with  the  case  as 
ImC>“. 

As  for  N<n<Nn,ax,  evaluation  of  W(n,h)  requires  two  separate 
treatments: 


1.  When  U-n|>n''3,  the  Debye's  asymptotic  approximations  of 

I  I 

Jn(5).  valid  and  can  be  used  to 

reduce  the  burden  to  computing  W(n,h).  More  specifically, 
one  can  write 
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-■  .'v 


.v>: 


SE 


2Csinh2YY- ^  n  tanhy  j  |  7  tanhy  (7  tanhY)  | 

r  J.  J_  1  21 

+  1  -  4n  (cothY  -  Scoth^Y)  +  3^2  (4  coth^Y  -  3  coth‘*Y 


385 

+  36  coth^Y) 


(4.105) 


The  asymptotic  approximations  of  J^(0  and  Jp(5)  are  naturally 

1  (11  (21 

derived  from  the  relation  that  0^(5)  =  [Vi  +  (5)]. 


When  l5-nl<n^''3^  so-called  transition  region,  the  asymptotic 
approxmations  of  0n(5),  8^(5)  and  their  corresponding  derivatives 
take  on  different  forms  which  are  described  in  [26]  and  will  not  be 
repeated  here.  It  may  be  noted  that  in  this  region,  both  n  and  5 
are  sufficiently  large.  Again  employing  these  approximations  to 
compute  W(n,h)  will  significantly  improve  the  efficiency. 
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As  a  consequence,  it  is  convinced  that  the  mutual  impedance  between 
mode  1  and  mode  2  is  reasonably  approximated  by 


«  max  sin{n(|)i)  sin(n(|i2) 

"12  '  ‘l2  /  Jo  ~  " - 


•  cos(n4.2c)  cos(hz2c)  •  Zj(h)Z2(h) 


(4.106 


b.  Appendix  G  suggests  that  condition  (4.94)  will  restrict  the 
contour  of  integration  to  the  fourth  quadrant  in  the  complex  h-plane. 

In  Figure  4.2,  the  broken  lines  denote  the  appropriate  branch  cuts 
starting  from  the  branch  point  h=k2(-k2)  going  to  h=0  on  the  real  axis 
and  then  to  h=j»(-j«).  Quite  significantly.  Appendix  H  indicates  that 
both  Pp/Sn  and  Qp/Tp  are  odd  functions  of  u.  Hence,  it  follows  from 
(4.85)-(4.87)  that  W(n,h)  is  an  even  function  of  u.  Thus  no  branch  cut 
is  required  for  the  square  root  function  u,  i.e.,  k^  is  not  a  branch 
point. 

c.  There  is  at  least  one  pole  for  each  n  on  or  above  the  real 

h-axis  (depending  on  lossless  or  lossy  dielectric)  between  h=k2  and 
h=Rek^.  The  original  integral  for  (4.106)  can  be  computed 

along  the  contour  r  which  runs  along  the  real  h-axis  and  overpasses  the 
the  poles  with  infinitesimal  half  circles  as  shown  in  Figure  4.2.  Such 
contour  requires  the  knowledge  of  the  pole  for  each  n,  0<n<Nn,aj<*  To 
avoid  the  numerical  burden  of  locating  these  singular  points,  it  is 
computationally  more  suitable,  and  justifiable  by  Cauchy's  theory  as 


well,  to  use  the  deformed  contour  Fq  which  starts  from  h=0  going  down 
the  negative  imaginary  axis  to  some  point  h=hi  sufficiently  away  from 
the  singularities.  Next  Fq  moves  parallel  to  the  real  axis  to  some 
point  h=h2,  with  Reh2>Reki,  then  travels  vertically  upward  to  meet  the 
real  axis  at  h=h3.  The  final  segment  of  Fq  extends  from  h3  to  h=».  The 
contour  Fg  is  also  shown  in  Figure  4.2.  It  should  be  cautioned  that  one 
cannot  let  Imh^  become  too  negative  since  then  the  cosine  functions  in 
(4.106)  (e.g.,  cos(hzj^)  and  cos(hzjjc),  4=1,2)  become  exponentially 
large.  This  will  unfortunately  cause  numerical  difficulty  in 
accurately  evaluating  zi2  which  is  typically  very  small. 


d.  The  integral  in  (4.106)  can  be  evaluated  in  three  stages  as 
symbolically  written  below: 


H 


■12 


^  Q  dh  +  Q  dh  +  /  0  dh  , 


(4.107) 


where  H2>Hi>Reki  and 
N 


max  sin(n(|»i)  sin(n()>2) 

‘X"’ '  *12  ~  ;; 


•  cos(ni.2^)  cos(hZ2c)Zj(h)Z2(h) 


(4.108) 


Stage  1  refers  to  the  integration  from  h=0  to  h=Hi.  The  upper  limit 
Hi  is  chosen  such  that 


The  integration  is  done  numerically  using  Simpson's  quadrature  rule. 

Stage  2  refers  to  the  integration  from  h=Hi  to  h=H2.  Because  of 
condition  (4.109),  for  all  h>Hi,  the  function  W(n,h)  can  be  greatly 
simplified  by  replacing  the  cylinder  functions  with  their  Debye 
asymptotic  forms.  Taking  the  leading  term  in  (4.99)  and  substituting  y 
by  jy,  one  can  deduce  from  (4.97)  and  (4.103)  that 


(2) 

(5) 


ir^siny 


±j[5(siny-ycosy)-rr/4] 


and 


1  '  1  i 

(2)  d  (2)  (2) 

^  ^^^“d5”n  (5)  ~  (5)  , 


(4.110 


(4.111 


where  y(n,5)  is  defined  by  cosy  =  n/5.  It  readily  follows  that 


j2[5(siny-ycosy)-ir/4] 

e 


(4.112 


2 


j2[C(siny-ycosy)-ii/4] 


hj 


n  I  n  ^ 

r(n,0  =  -jin  [  .  +  j  0  -  (  )  ] 


(4.114) 


Employing  (4.112)  -  (4.114)  will  enable  one  to  obtain 


^  eM 

5  ~  jsinY(n,iib)  — 

^n  e^-l 


(4.115) 


'n 


e^-1 


T  ~  jsinY(n,vib) 
n  e^+1 


(4.116) 


where 


X  =  j2{wb[sinY(n,ub)-Y(n,nb)cosY(n,Mb)] 

-  ua[sinY(n,ua)-Y(n,ua)cosY(n,ua)]}  .  (4.117) 

It  is  now  quite  apparent  that  these  asymptotic  forms  will  reduce  the 
complexity  of  W(n,h)  many  folds. 

Finally,  Stage  3  refers  to  the  integration  from  h=H2  to  h='».  H2 
can  be  chosen  so  that  X=u=jh  and  jc|>>n,  with  ?=ua,  ub  or  Xb. 
Consequently  one  sees  that  Y(n,5)=Ti/2,  cosY(n,5)=0,  and  sinY(n,5)=l. 

One  can  further  approximate  that  Pn/Sn^-j,  On/Tn=-j  and  rpsj.  All  these 
approximations  greatly  simplify  W(n,h)  to 

~  j2h 

*  irb(l+G  )  » 


(4.118) 


which  is  independent  of  n.  In  addition,  noting  that 


cos(hZi)-cos(<iZi) 

Zl(h)  =-  - ^2 - 


(4.119) 


~  C0S(hZ2 )-COS( <2^2) 

Z2(h)  =  -  ^ 


(4.120) 


and  making  use  of  trigonometric  identities,  one  can  write 

"  j2Ai2  Nmax  sin(n(t.i)  sin(n()»2)  9 

/ ■  ;S(TTr)  !  s  - r-  “'(''♦20'  J  >,  • 

^2  ^  n**u  *  A 

(4.121) 


where 


'2  “  /  cos(b^h)  ^3  >  t*l,2,...,9 

^2 


(4.122) 


^1  =  4  ,  bi  =  |Zj-Z2!  -  Z2^ 


(4.123) 


^2  =  4  »  ^2  ^ 


l^r^2'  ^2c 


(4.124) 


‘3  =  4  •  1^3  =  |Zl^Z2-^2cl 


(4.125) 


(4.126) 


“4  '  4  •  '>4  -  "i*42*Z2c 


COS(kiZi) 


=  - - T 


.  bg  =  . 


(4.127) 


^6  = 


COS(KiZi) 


•  ““e  =  ^2*^2c 


(4.128) 


cos  ( *^2^2  ) 


»7  =  -  - ^ 


.  ^  =  U1-22J 


(4.129) 


COS(<2Z2) 

“8  •  -  I  .  >>8  ■  2j.22^ 

ag  =  cos(<iZi)cos(<2Z2)  ,  bg  =  Z2c  . 


(4.130) 


(4.131) 


Finally,  employing  an  integration  by  parts,  one  can  verify  that,  for 
^~1 ,2, ...,9, 


/  cos(b  h) 

Ho  ^ 


h3  = 


=  0 


cos(b|iH2) 


bj^sin(b^H2) 


2 

Hz 

+  2  ci(bj^H2), 

\ .  0 


(4.132) 
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(4.133 


is  known  as  the  cosine  integral. 

For  illustration  purposes,  two  numerical  examples  of  coupling 
between  identical  expansion  modes  are  considered  using  the  theory 
developed  in  this  section. 

Example  1:  E-plane  coupling 


The  geometry  of  this  example  consists  of  two  identical  modes,  of 
length  L  =  0.15  m.  and  width  W  =  0.075  m.,  on  a  coated  cylinder  of 
radius  a  and  dielectric  coating  of  thickness  t  =  0.003175  m.  with 
He  =  <t>2c  “  0  ^Ic  ~  separation,  S  =  Z2c,  between  the  modes 

varies  in  the  z-direction  only.  The  operating  frequency  is  633  MHz, 
with  Ep  =  2.56  and  tan5  =  0.0015. 

In  Table  4.1,  results  are  given  for  the  self  impedance  of  a  single 
expansion  mode  (zn)  on  cylinders  of  various  radii.  It  is  evident  that 
zji  converges  to  the  infinite  dielectric  slab  value  as  the  cylinder 
becomes  large.  This  is  a  good  indication  that  the  theory  and  analysis 
techniques  presented  here  are  valid.  Table  4.1  also  shows  the 
convergence  of  the  summing  series  (Nn,ax)t  contour  of  integration 
Tq  (hj  and  h3'  and  the  limits  of  the  three  stages  of  integration 
(Hi  and  H2). 


3 


TABLE  4.1 


Lti 


o 


SELF  IMPEDANCE  OF  A  SINGLE  EXPANSION  MODE 


Radius  a 
(^o) 

— 

Self 

Impedance  zn 
(volt-amp. ) 

^max 

hi 

h3 

Hi 

H2 

0.25 

0.145-j0.429 

62 

-jO.2 

3.0 

50.0 

70.0 

0.5 

0.118-j0.422 

121 

-jO.l 

3.0 

50.0 

70.0 

1.0 

0.104-j0.418 

240 

-j0.15 

3.0 

50.0 

70.0 

3.0 

0.0957-j0.415 

717 

-jO.l 

3.0 

50.0 

80.0 

10.0 

0.092-j0.419 

1998 

-j0.05 

3.0 

55.0 

70.0 

00 

0.092-j0.399 

- 

- 

- 

- 

Mutual  impedance  (Z12)  of  the  expansion  modes  is  plotted  in  Figure 
4.3  as  a  function  of  separation  between  the  modes  for  cylinders  of  radii 
a  =  0.25Xo,  0.5X0.  l.^o  3nd  IOXq,  where  \q  denotes  the  free  space 
wavelength.  It  is  seen  that  the  coupling  rapidly  weakens  for 
S>0.3Xq.  Generally  this  is  in  agreement  with  the  planar  dielectric  slab 
case  studied  in  Chapter  II.  It  is  interesting  to  observe  that  the 
degree  of  coupling  between  two  expansion  modes  in  the  E-plane 
configuration  varies  slightly  among  cylinders  of  different  radii.  This 
leads  one  to  conclude  that  the  curvature  of  the  cylindrical  surface  has 
little  effect  on  the  energy  propagating  in  the  z-direction  as  surface 
wave  modes.  The  values  of  2^2  are  also  listed  in  Table  4.2. 

Example  2:  H-plane  coupling 

This  example  repeats  the  data  of  Example  1  except  that  zic=2:2c®0» 
^ilc=0,  and  S=b'ti2c  is  separation  between  the  expansion  modes  which 
varies  in  the  (|i-di rection  only. 

The  results  are  plotted  in  Figure  4.4  and  tabulated  in  Table  4.3. 

It  is  apparent  that  the  coupling  in  the  H-plane  configuration  gradually 
reduces  as  the  radius  of  the  cylinder  increases.  This  can  be  explained 
by  the  fact  that  waves  propagating  circumferentially  shed  more  energy 
into  the  surrounding  space  when  the  curvature  is  large.  One  may  notice 
that,  in  the  a=0.25XQ  case,  coupling  is  reinforced  at  S=1Xq  which  is 
close  to  the  back  side  of  the  cylinder  relative  to  expansion  mode  1. 

This  is  due  to  the  fact  that  waves  which  originate  from  mode  1  and 


s 

MUTUAL  IMPEDANCE  (volt-amp.)  j 

^  0-' 

a=0.25XQ 

3i“0  •  5  Xq 

a=1.0XQ 

a=10.0XQ 

0.1 

0.119-j7.528 

0.095-j7.491 

0.082-j7.472 

- 

0.25 

0.068-j5.578 

0.055-j5.555 

0.048-j5.542 

0.041-35.523 

0.32 

0. 0504- j. 1981 

0.04331-3.2028 

0.03894-3.2037 

- 

0.4 

0.023+j0.016 

0.023+j0.007 

0.022+30.002 

- 

0.5 

-0.005+j0.026 

O.OOl+j.020 

0.005+j0.015 

0.006+30.014 

0.6 

-0.021+j0.014 

-0.013-j0.013 

-0.009+30.014 

- 

0.75 

-0.02n-j0.010 

-0.016-j0.004 

-0.014-30.002 

-0.011+30.001 

0.9 

-0.002-j0.021 

-0,004-j0.014 

-0.004-30.010 

- 

1.0 

0.011-j0.016 

0.005-j0,013 

0.001-30.009 

0.0-30.008 

Note:  Xq  =  free  space  wavelength  =  0.474  (m) 
S  =  separation. 


MUTUAL  IMPEUANCE  BETWEEN  TWO  IDENTICAL  EXPANSION  MOOES  IN  THE  H-PLANE 


Figure  4.3.  E*plane  coupling  for  Identical  expansion  modes  on  coated 
cylinders  of  radii  e  •  O.ESXq,  0.5U,  l.L.  and  10. U 
(V)  -  0.474  M.). 


>-v 

L-  ■ 


142 


(VOLT-AMPERE) 


0.0  0.2  0.^  0.6  O.e  1.0  l.i 

S(Xo) 

Figure  4,4.  H-plane  coupling  for  Identical  expansion  modes  on  coated 
cylinders  of  radii  a  •  O.ZSXq.  O.BXq,  I.Xq  and  S.Xq 
(\q  ■  0.474  m.). 


propagate  circumferentially  in  opposite  directions  interact  more 
strongly  when  the  paths  of  propagation  are  short  (as  for  the  case  of 
small  cylinders),  and  weakly  when  the  paths  of  propagation  are  long  (as 
for  the  case  of  large  cylinders)  because  of  the  loss  of  energy  along  the 
paths. 

Figures  4,5  to  4.7  show  that  E-plane  coupling  is  in  general 
stronger  than  that  of  H-plane.  This  is  also  attributed  to  the  shedding 
of  energy  away  from  the  surface  of  the  coated  cylinder  by  the 
circumferentially  propagating  waves.  Even  though  there  is  no  curvature 
effect,  it  is  interesting  to  recall  that  the  H-plane  coupling  is  also 
weaker  than  the  E-plane  coupling  in  the  case  of  the  planar  structure. 

E.  FAR  FIEin  CALCULATION 

M  A  A  ^ 

If  the  current  distribution  J(r')  »  (j'J^i(r')  +  zJ^,(r')  is  known 
for  a  microstrip  antenna,  then  the  field  can  be  computed  for  an 
arbitrary  r  e  V^.  In  practice,  it  is  of  great  interest  to  know  the 
radiation  pattern  of  the  antenna.  To  this  end,  this  section  will  derive 
a  general  expression  via  the  method  of  steepest  descent  to  facilitate 
the  far  field  calculation. 

By  studying  Equations  (4,61)  to  (4,63),  it  becomes  apparent  that 
any  electric  field  component,  due  to  J,  can  be  expressed  in  the 
following  form: 


I  T(F,n)  , 


AMPERE) 


o.e 


1.0 


0. 

Figure  4.] 


— .  •  •  •  I  I  •  •  •  I  ■  I  >  I  I  I  ■  I  I 

0  0.2  O.M  0.6 

SlXg) 

Coupling  for  Identical  expansion  modes  on  a  coated 
cylinder  of  radius  a  •  I.Xq  {X^  ■  0.474  m.). 
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T(r,n)  =  /  f(4',n,h) 

r 

f(<l',n,h)  =  g(n,h)  /  JJr') 

S 

a  =  (j)',  2*,  S  denotes  the  surface  of  the  antenna,  and  g(n,h)  involves 
the  expansion  coefficients  of  the  Green's  dyadic.  Furthermore,  f(()),n,h) 
is  an  analytic  function  of  the  complex  variable  h  along  the  path  of 
integration  r  (Figure  4,8)  with  end  points  at  infinity.  The  choice  of 
r  is  discussed  in  Appendix  G. 

In  the  far  zone,  |xpl  =  ] p/k^-h^  |  *,  one  can  employ  the  large 

argument  approximations  for  the  Hankel  functions: 


cos 


sin 


e'^^^  bdi,'dz'  ,  (4,135) 


H^^^(Xp) 

h1^^‘(Xp) 


e-^^^  dh 


(4.134) 


(4.136) 


and 


X  .  j(^P  -  nir/2  -  n/4) 

^n  ~  j  V  irXp  ®  •  (4.13 

Using  these  approximations  in  (4.134)  reveals  that  it  is  sufficient  to 
consider  integral  of  the  type: 
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T(f  ,n) 


1  f(<t',n,h) 
r 


gj(Xp  -  nir/Z  -  iT/4)+jhz 


dh 


(4.138) 


which  can  be  evaluated  asymptotically  using  the  method  of  steepest 
descent  as  |  Xp|  -►  *. 

Introducing  the  transformation 

h  =  -kjcosw  (4.139) 

and  realizing  that  z  =  rcos0,  p  =  rs1n9,  where  r  and  9  are  the  usual 
spherical  coordinates,  one  finds  that 

P  \/ k|-h^  +  hz  =  -  k2rcos(w+9)  .  (4.140) 

Hence,  (4.138)  can  be  rewritten  as 


•'2rq(w) 

T(r,n)  =  J  F((|»,n,w)  e  dw  ,  n=0,l,2,...  ,  (4.141) 

r 

w 

where  the  contour  r  in  the  h-plane  is  mapped  onto  the  contour  in  the 
w-plane  as  shown  in  Figure  4.8, 

q(w)  =  -jcos(w+9)  ,  (4.142) 


and 


F(K-..w)  =  *  ’'“V  f(».n.-kacosw) 


(4.143) 


.  ) 
v-'i 

.  o 


1.4 


*  1 


J 
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The  saddle  point  is  given  by  solving 


q'(w)  =  jsin(w+0)  =  0  .  (4.144) 

w=w 

s 

It  is  clear  that  Wj  =  n-9  is  the  appropriate  saddle  point.  Observing 
that 

q(ws)  =  -jcos(n-0+9)  =  j  (4.145) 

and 

q"(ws)  =  jcos(Tr.9+9)  =  -  j  ^  0;  ,  (4.146) 

one  realizes  Imq(w)  =  1  along  the  steepest  descent  path  (SDP),  and  Wg  is 
a  first-order  saddle  point  since  q''(Ws)?‘0.  Letting  w=Wp+jwi ,  finds, 
along  the  SDP, 

Imq(w)  =  Im{-jcos(wp+9+jwi ) }  =  -  cosh(wi )cos(wp+9)  =  1  .  (4.147) 

As  cosh (wi  )■*■<»,  thus  cos(Wp+9)-»-0  in  order  to  satisfy  (4.147). 

This  suggests  that  the  SDP  is  asymptotic  to  Wp=tt/2-9  and  Wp=3ii/2-0.  In 
addition,  introducing  the  real  variable  B  via  the  transformation 

q(w)  =  q(ws)-0^  ,  (4.148) 

leads  to 

dw  2S 

^  =-q'(w)  „  •  (4.149) 

3=0  ^  3=0 
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Upon  applying  L'Hopital's  rule  to  the  indeterminate  form  for  dw/dS  in 
(4.149)  when  0=0  (i.e.,  w=w^),  one  obtains 


dw 

d0 


-2 


0=0 


q"(w5) 


=  /r  & 


-jn/4 


(4.150; 


Equation  (4.150)  implies  the  SOP  makes  an  angle  of  45°  with  the  real 
w-axis  at  w=Ws.  The  SOP  is  depicted  in  Figure  4.9. 

Finally,  a  general  expression  for  the  first-order  approximation  of 
T(7,n)  is  given  by  [28]: 


T(r,n) 


-2w 


k2rq"(w^) 


f(4'.n,w^)  e 


k2rq(w  ) 


,  k2r 


•  (4.151 


Combining  (4.136),  (4.137),  (4.138)  and  (4.151)  yields 


T(r,n)  = 


/  f(<^.n,h)  j^y 


h!^\xp) 


h;*'  (Xp) 


e^^"  dh 


jk2r 


,  k2r 


•  (4.152 


This  completes  the  fundamental  portion  of  the  far  field  calculation 


CHAPTER  V 


ASYMPTOTIC  SOLUTIONS 


A.  INTRODUCTION 

In  the  last  chapter,  all  electric  fields  and  impedances  are  exact 
eigenfunction  solutions  since  they  are  derived  by  the  method  of  Green's 
function  which  itself  is  expressed  in  terms  of  eigenmodes  of  the 
dielectric  coated  cylinder.  These  eigenfunction  solutions  are  employed 
in  Chapter  IV  to  calculate  the  mutual  coupling  associated  with 
microstrip  arrays.  In  many  situations,  however,  microstrip  antennas  may 
be  used  in  the  construction  of  phased  arrays  that  are  flush-mountable  on 
large  cylindrical  objects  such  as  an  aircraft  fuselage  or  space 
vehicles.  The  coupling  between  antenna  elements  plays  a  significant 
role  in  the  design  of  these  arrays.  In  such  situations  where  the 
cylinder  is  large  in  terms  of  wavelength,  it  is  best  to  use  asymptotic 
high  frequency  methods,  as  opposed  to  an  eigenfunction  series  solution, 
to  analyze  the  coupling  effect  because  the  eigenfunction  solution  is 
poorly  convergent.  Moreover,  the  radiation  mechanism  is  often  obscured 


in  an  eigenfunction  solution.  On  the  other  hand,  the  asymptotic 
solutions  are  far  less  complicated  from  the  viewpoint  of  computational 
effort  and  cost  involved.  More  attractively,  the  asymptotic  solutions 
may  be  interpreted  in  terms  of  ray-optics  thereby  revealing  the 
radiation  mechanism  involved. 

Asymptotic  high  frequency  solutions  usually  result  from  asymptotic 
evaluations  of  the  integral  representations  of  Green's  functions  of  the 
circumferentially-propagating  type.  Section  B  will  discuss  the 
alternative  representations  of  the  dyadic  Green's  function  for  the 
dielectric  coated  cylinder.  It  includes  the  integral  representation  of 
the  circumferentially-propagating  Green's  function  from  which  one  can 
derive  the  residue  series  representation  that  provides  nice  ray-optical 
interpretation.  Section  C  devotes  to  the  derivation  of  alternative 
representations  of  the  mutual  impedance  zi2  between  mode  1  and  mode  2 
considered  in  the  numerical  examples  of  Chapter  IV,  For  simplicity,  a 
two-dimensional  (2-D)  coupling  problem  involving  two  infinite  conducting 
strips  is  formulated  and  analyzed  asymptotically.  A  numerical  example 
of  the  2-0  problem  is  treated  in  Section  D.  It  serves  to  illustrate  the 
validity  and  usefulness  of  asymptotic  solutions. 

B.  ALTERNATIVE  REPRESENTATIONS  OF  THE  GREEN'S  DYADIC 

Integral  transforms,  such  as  the  Laplace,  Fourier,  and  Mellin 
transforms,  can  be  found  useful  in  many  cases  in  converting  relatively 
slowly  converging  series  into  much  more  rapidly  converging  ones.  Upon 
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examining  the  series  associated  with  the  dyadic  Green's  function  which 
converge  slowly  as  the  cylinder  radius  increases,  it  appears  that  the 
Poisson  summation  formula,  which  is  based  on  the  Fourier  transform,  can 
convert  these  poorly  converging  series  into  fast  converging  ones,  at 
least  in  principle.  One  form  of  the  Poisson  summation  formula  can  be 
stated  as  follows: 

I  I  ,  (5.1) 

n=-<»  n=-<» 

where 

oO 

F(0  =  /  f(v)eJ^^dv  ,  (5.2) 

•  eo 

and 

1  ** 

=  YU  j  F(5)e"'^^'’d5  (5.3) 

•  os 

form  a  Fourier  transform  pair,  and  a  is  a  scalar. 

The  utility  of  the  Poisson  summation  formula  in  converting  a 
slowly  converging  series  into  a  rapidly  converging  one  may  be 
appreciated  from  the  properties  of  the  Fourier  transforms.  If  f(v)  is  a 
function  which  is  concentrated  near  v=0,  then  F(t)  is  spread  out  over 
the  C-axis,  or  has  an  appreciable  value  for  a  large  value  of  5;  for 
instance,  if  f(v)=6(v),  then  F(f)=l,  Conversely,  if  f(v)  is  spread  out 


along  the  v-axis,  then  F(c)  is  localized  near  5=0.  A  slowly 
converging  series  with  a  typical  term  f(an)  is  then  converted  to  a 
rapdily  converging  series  with  a  typical  term  F(2mr/a),  since  F  will 
be  small  for  large  values  of  n  while  f(an)  decreases  slowly  with 
increasing  n. 

From  Equations  (3-123),  (3-125)  and  (3-126)  one  observes  that  every 
component  of  the  Green's  dyadic  is  composed  of  functions  of  the  form 


«  00 


where 


“  g  (p,p';n,h)cosn((t.-().') 

S(  p,p'  ;h)  *  1  e  ® 

"=0  gg(p,p';n,h)sinn((|)-(|)') 


(5.4) 


with  g^  and  g^  being  even  and  odd  functions  of  n,  respectively.  Thus 
one  only  needs  to  consider  S  in  looking  for  the  alternative 
representations  of  the  Green's  dyadic. 

The  series  S  in  (5.4)  can  be  rearranged  such  that 


S(p,p'  ;h)  =  I 


ge(p.p';n,h)  I  jn(4,-4,') 


n=-<»  -jgg(p,p' ;n,h) 


(5.5) 


In  addition,  employing  the  large  order  approximations  of  the  cylinder 
functions : 


(5.6) 


(2) 

(z)  ~ 


ez 


V 


(5.7) 


(5.8) 


Hy  (z)  ~  -  (2)  .  (5.9) 

with  z,y  being  complex  numbers,  one  can  easily  verify  that  both  g^  and 
9g  are  square  integrable  over  the  interval  (-•,<•),  i.e.. 


/ 


l9g( p,p' 
l9o( PtP' 


dv  <  • 


(5.10) 


Clearly  (5.10)  also  implies  l9g|,  Iggl^o  as  v+-«».  These  properties  of 
boundedness  warrant  the  existence  of  the  Fourier  transforms  of  g^  and 
9q.  Hence  the  Poisson  summation  formula  can  be  applied  to  the  series  in 
(5.5)  to  yield 


ge(p.p' ;n,h) 
-jgo(p»p' 


gg(p.p'  ;'>.h) 
'jgQ(p,p';v,h) 


gj'>(*-<t>')  gj2nnv 


=  I  I 

ps.oo  .00 


ge(p.p' ;v,h) 
-jsgn((|»-<|i'  )g^(p,p‘  ;v,h) 


gj'^l  1  +  j2mrv 


where 


sgn(x) 


1  ,  X  >  0 
1  ,  X  <  0 


Rearranging  (5.11)  will  result  in  the  following  form; 


gg(p.p' ;n,h) 
-jgQ(p.p' ;n,h) 


gjn(t-<t>') 


gg(p.p' ;v,h) 

-jsgn((|.-<^')g^(p,p';v,h) 


|+2nn)  jv(2it-|(^-y  |+2n7T) 

I  e  16  n  \j 


The  original  representation  which  appears  on  the  L.H.S.  of  (5,13)  is 
associated  with  the  "radially  propagating"  Green's  function.  Whereas 
the  alternative  representation  which  appears  on  the  R.H.S.  of  (5.13)  is 
found  in  the  "circumferentially  propagating"  Green's  function  which  is 
particularly  suited  to  emphasize  the  ray-optic  behavior  of  radiation  in 
the  shadow  zone  of  the  cylinder. 

Since  the  integrand  of  (5.13)  vanishes  as  in  the  upper  half 
plane,  the  contour  of  integration  (along  the  real  u-axis)  can  be  closed 
in  this  half  plane  to  form  a  new  (closed)  contour  Cy,  as  shown  in  Figure 
5.1(a),  without  affecting  the  value  of  the  integral,  or  symbolically, 

00  00  00 

I  /  dv  =  I  /  dv.  (5^14) 

n=0  -«  n=0  Cy 

It  is  known  that  the  simple  poles  of  g  or  g  form  an  ordered  set 
{''ptRe Vp,ImVp  >  0,p=l,2,3,  in  the  first  quadrant  of  the  v-plane. 

One  can  use  Cauchy's  theory  to  evaluate  the  R.H.S.  of  (5.14)  to  obtain 
the  residue  series: 

Res[g^(p,p' ;Vp,b)] 
-jsgn(t-'fr’)ResCg^(p,p';Vp,h)] 


I  /  du  =  j2n  I 

n=0  Cy  n=0  p=l 
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This  residue  series  representation  is  Important  for  analyzing  the 
diffracted  fields  in  the  deep  shadow  region  of  the  cylinder. 


Physically,  each  term  in  the  series,  corresponding  to  p=l,2,...,  can  be 

interpreted  as  a  creeping  wave  mode  that  creeps  around  the  cylinder  n 

1 4)-4>'  |+2mi) 

times,  n=0,l,2,....  The  exponential  terms  e  ^  ,  and 

j  Vp ( 2 n I  ((>- (ji  |+2nir)  ,  .  . 

e  (which  corresponds  to  waves  creeping  around  the 

cylinder  in  the  direction  opposite  to  that  of  the  former)  become 

exponentially  decaying  terms  which  attenuate  faster  as  | <!>-♦' |  become 

larger.  Since  Imv^  increases  with  p,  the  higher  order  modes  are 

attenuated  very  rapidly.  So  consequently  one  would  probably  need  to 

retain  only  the  first  couple  of  terms.  In  addition,  for  large 

cylinders,  the  n=0  encirclement  is  sufficient  because  contributions  from 

higher  encirclements  are  negligible  as  the  wave  attenuates  while 

creeping  around  the  cylinder.  The  exponential  decay  of  the  creeping 

wave  modes  on  the  cylinder  suggests  that  energy  leaks  away  from  the 

cylinder  surface  into  the  surrounding  space. 


Analytically,  one  may  note  that  Imv^,  p=l,2,3...,  is  significantly 
greater  than  zero.  Hence,  contour  can  be  deformed  in  such  a  way 
that,  along  the  new  contour  C^,(p-jgyre  5.1(b)),  the  integration  variable 
V  always  has  lmv»0  thereby  rapidly  damping  the  encirclement  terms  with 
n>0.  Consequently  one  can  write 


Figure  5.1.  (a)  Contour  Cv  encloses  the  upper  half  v-plane  which 

contains  the  poles  of  or  g^:  V2,...,vp.. 

(b)  The  deformed  contour  also  encloses  all  the  p 

of  ge  or  g^. 


s(  p.p'  ;h)  =  I 

n=-no 


gg(p»p';n.h) 
-jgg(p»p' ;h,h) 


gg(p.p' 

-jsgn((t.-4.')g^(p,p';v,h) 


gjv|(|)-(fr'l  ^  gjv(2ir-l4.-({.'|)| 


~  j2n  I 
p=l 


Res  Cgg(p,p' ;Vp,h)] 

-jsgn(  ^-4.'  )Res[g^{p,p'  ;Vp,h)] 


jvpl^-i^'l  jvp(2if-(j)-!|>' 1 ) 
e  i  e 


(5.16) 


The  three  representations  of  S  In  (5.16)  are  of  different  utility  in 

evaluating  S.  The  eigenmode  solution,  the  first  form,  is  exact  and 

suitable  for  small  cylinders  numerically.  The  remaining  two  forms  are 

valid  for  large  cylinders.  In  particular,  the  integral  representation 

can  be  used  when  |4»-4i‘l  is  small  while  the  use  of  the  residue  series  is 

^0 

most  appropriate  when  U-ffr'l  is  large  (in  general,  |4'-i>'|>Ib.  is  the 
free  space  wavelength). 
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C.  ASYMPTOTIC  EVALUATION  OF  212 


As  the  radius  of  the  coated  cylinder  increases,  the  asymptotic 
expansions  of  fields  and  impedances  are  deemed  appropriate.  The  mutual 
impedance  212  between  mode  1  and  mode  2,  both  2-polari2ed  as  shown  in 
Figure  4.1,  obtained  in  Chapter  IV  is  an  eigenfunction  solution  which  is 
not  suitable  for  asymptotic  expansion.  It  is  desirable  to  derive  an 
integral  representation  of  2^2  via  the  use  of  the  Green's  dyadic  of  the 
circumferentially-propagating  type.  Adopting  the  same  notations  used  in 
Chapter  IV,  the  2-component  of  the  surface  electric  field  due  to  mode  1 
can  be  directly  obtained  from  (4,63)  as: 

Ei\(b,^,2)  =  -  2  I  ds'  J  dh  I 

ankj  s  —  n=0 

1 

e^W(n,h)cosn(<|»-<|>')J^j(b,i)' ,2' )  ^  ^  ,  (5.17) 

where  Ji2(b,i),2)  and  W(n,h)  are  defined  in  (4.76)  and  (4.85), 
respectively,  and 

Sj  =  {r=(b,i.,2):-i)j<(|.<(t.j,-Zj<2<2j}  .  (5.18) 

By  definition, 

^12  “  -  i  E2z^^’.^*z)J22(b,t,z)ds  ,  (5.19) 

^2 
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where  J2z(b.<J».z)  is  again  defined  in  (4.76),  and 


$2  =  {r  =  (b,i).z)  : 


Jh(z-z') 


Combining  (5.17)  -  (5.20)  one  can  explicitly  write 


(ijn  Z  +Z  ^  Z  ^ 

2  2c  2  2c  2  1  1 

^  ^2  /  ‘iz  /  bd(ti  /  dz*  /  bdi|)'  /  dh 

2  Z2(.-Z2  <t>2c-'^2  '^1  ~h 


I  %W(n,h)cosn(t-(t.')J,  (b,t',z')Jp_(b,4.,z) 


Since  W(n,h)  is  an  even  function  of  n,  it  follows  from  (5.4)  and 


(5.5)  that 


jh(z-z') 


o  /  dz  /  bd(|>  /  dz'  /  bdi}>'  /  dh  e 

lirk^ 

^  ^2c"^2  *2c"‘^2  “^1  ‘h 


I  W(n,h)eJ"^‘*’"‘*’  ^0,,(b,<t.',z')Jp  fb,4,z) 


Carrying  out  the  integrations  over  the  variables  <J),iJ>',z  and  z'  yields 


^  -  7  /hW  V  MlHilll  I  ^^(1(120+ 1^2)  Jn('()2c-'t'2)  I 

Zj2  =  8  i  dh  Zj(h)Z2(h)  I  e  -  e 


I  -jfiiii  jn^ii 

•  I  e  -  e 
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where  2=1,2,  and  ^re  specified  in  (4.81)  and  (4.83), 

respectively.  According  to  (5.13),  (5.14)  and  (5.16),  i*’  (5.23)  can 

be  written  as  follows: 


^2="r  I  dh  Z^(h)Z^(h)  I  J  W(v,h) 
-00  n=0 


+  ^^(v)  -  r^(v)  -  r^(v)  I  dv 


Ai2  °®  ~  ~ 

=  -r  /  dh  Z^(h)Z2(h)  I  J  W(v,h)  . 
-«>  n=0  Cy 


+  r2(v)  -  r3(v)  -  r4(v)  dv 


^12  ^  ® 

-T'  /  dhZ^(h)Z2(h)  /W(v,h)  • 


ri[v)  +  r2(v)  -  r3(v)  -  V^{v)  i  dv  ,  (5.24) 


where 


jv(|t2c+<^2-'^ll+?nTi) 

ri(v)  »  e 

jv(2iT-|  ij)2c+i>2-i*l  1  ■*'2mi) 

+  e 

.  (5.25) 

jv(|t2c-'('2+<('il+2nir) 

j  v(  2ir- 1  i>2c- ('2'*’ '<’1 1  ^ ) 

r2(v)  =  e 

+  e 

,  (5.26) 

^  jv(|i'2c-<^2-'t'l|+2nit) 

jv(2Tr-|  4>2c-'}>2-4»il+2mT) 

r3(v)  »  e 

+  e 

,  (5.27) 

J ''( 1  i'2c+<t'2+<l'l!+2nTr) 

jv(2w-l  '|)2c+  <>2+<l>lK2mi) 

/  C  OQ  ^ 

and  Cy  is  the  contour  enclosing  all  the  poles  of  W(u,h)  in  the  upper 
half  v-plane.  The  closed  contour  integral  in  (5.24)  can  be  readily 
evaluated  by  residue  calculus  to  yield 


j2TrAi2  “ 


0»  90 


’12 


8 -  /  dh2^(h)Z2(h)  I  I  Res  [W(v  h)] 


n=0  p=l 


•  ]  ''I'V  "  -  --sIV  - 


j2irAi2  “ 
^  / 


dhZ^{h)Z2(h)  I  Res  [W(Vp,h)]  • 
p=l 


(V 


+  r 


2^^p^  ■  ^3 


(V- 


(5.29) 


As  indicated  in  the  last  section,  the  alternative  representations  of  ziz 
are  derived  on  the  basis  that  W(v,h)-*-0  as  v+»,  and  W(v,h)  is  analytic  on 

the  upper  half  v-plane  except  at  the  poles  Vp,  p=l,2,3 .  Also,  the 

approximations  in  (5.24)  and  (5.29)  are  appropriate  for  large  cylinders. 
It  may  be  noted  that  the  ray-optical  interpretation  of  the  residue 
series  discussed  in  the  last  section  still  applies  here. 

As  far  as  evaluating  the  residue  series  in  (5.29)  is  concerned, 
only  the  first  creeping  wave  mode  will  be  sufficient  for  reasonable 
approximation  due  to  rapid  convergence  of  the  series  when  the  two  modes 
are  well  separated  (one  or  more  free  space  wavelengths  apart).  As  their 
separation  reduces,  the  series  is  becoming  poorly  convergent  and  the 


number  of  creeping  wave  modes  will  hence  increase.  Upon  examining  the 
complexity  of  W(\>,h)  (equations  (4.84)  -  (4.90)  with  v  replacing  n),  it 
appears  that  there  is  no  known  way  to  locate  the  poles  of  W(v,h) 
analytically.  One  needs  to  resort  to  numerical  techniques  such  as  the 
Newton-Rhapson  procedure  to  find  these  poles.  Such  task  is  by  no  means 
trivial.  Thus,  for  small  separations,  it  is  more  appropriate  to 
evaluate  the  integral  representation  of  2^2  in  (5.24). 

To  demonstrate  the  usefulness  and  power  of  the  asymptotic 
solutions,  it  is  best  to  consider  the  mutual  coupling  between  two 
infinite,  z-oriented  conducting  strips  on  the  dielectric  coated  cylinder 
whose  cross-section  is  shown  in  Figure  5.2.  Strip  Z,  centered  at 
t=l,2,  with  tic'Oi  has  constant  and  z-directed  current  J|^,  i.e.. 


Jj^(r)  =  Jj^(b,^,z)  = 


1 

2  2bH  * 

# 

0  ,  elsewhere  (5.30) 


The  mutual  impedance  Z12  between  strip  1  and  strip  2  can  be  expressed  as 


zc  _  _  _  zc  I 

^12  =  “  /  E2(r)-J2(F)dt  =  -  /  E2^(t)'^bd4.  , 

^2c"^2  ^2c‘'^2  (5.31) 

where  Ep{r)  denotes  the  electric  field  at  the  surface  of  the  dielectric 
coated  cylinder  (or  rigorously  at  P=b'*’),  and  its  z-component  is  given  by 
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CIRCULAR  CONDUCTOR 


Figure  5,2,  Cross-section  of  the  dielectric  coated  cylinder  with  two 
Infinite  conducting  strips. 


4(*)  =  - 


- -  J  bd^'  /  dz*  /  dh 

aitk^  , 


"  1 

I  £nW(n,h)cosn( (!)-((.')  2bn  *  (5.32) 

n=0 


Making  use  of  the  results  that 


/  dz'  =  27re'^^^5(h)  ,  (5.33) 


and 


/  f(h)5(h)dh  =  f(o)  ,  (5.34) 


one  can  readily  convert  (5.32)  to 


u>u. 


£*,(♦)  = 


n=— 


«  W(n,o) 

I  -T— 


jn(j»i  -jni))i 
e  -  e 


jntj) 


(5.35) 


Substituting  (5.35)  into  (5.31)  yields 


‘^''2  0.  W(n,o) 

^  2 
n=-“>  ^ 


jn( 'j>2c"'<’2'^4’l)  jn(  i|>2c+'<>2“‘l*l) 

e  +  e 


^jfi(  4’2c-4'2-^'1)  _  g:5n((fi?c+'l>2'*'  h) 


(5.36) 
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A. 


I  ILa  W  Tti  ^iffTI  ^rirtl 


Next  upon  recalling  Identities  (5.13)-(5.15)  one  will  arrive  at  the 
following  representations  of  Z12: 


*^^2  “»  *  W(v,o) 

"1 —  I  f  ~r~ 

16k2^ji^2  n=0  ^ 


ri(v)  +  r2(v)  -  r3(v)  -  r^(v) 


dv 


j  ir'JU2 


^**‘2'**l'**2 


I  Res 
p=l 


W( vp.o) 


where  rji(v),  )i=l,2,3,4,  are  defined  in  (5.25)-(5.28)  and  vp,  p=l,2, 
are  the  poles  of  W(v,o)  in  the  upper  half  v-plane. 

As  just  explained  earlier,  for  large  cylinders,  reasonable 
approximation  can  be  attained  by  ignoring  terms  which  correspond  to 
creeping  waves  making  n>0  encirclements  around  the  coated  cylinder, 
i  .e. , 
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It  has  been  pointed  out  that  the  manner  in  which  z\2  in  (5.48)  is 
being  evaluated  depends  on  the  separation  between  the  two  strips.  The 
quantity  b(ji2c  is  a  measure  of  separation  between  the  strips.  When  b4)2c 
is  small  (typically  less  than  a  quarter  of  free  space  wavelength),  the 
integral  representation  of  zi2  is  most  suitable.  When  bi^^c  ’s  large 
(typically  greater  than  a  free  space  wavelength),  the  residue  series  is 
an  efficient  representation  of  212*  Each  case  is  dealt  with 
specifically  in  the  following: 

Case  (a):  b-ji2c  is  small . 

The  exponential  terms  of  the  form  ej''(2ir-UI)  in  the  functions 
Ta(v),  t  =  1,2, 3, 4,  can  be  dropped  without  sacrificing  significant 
accuracy.  Hence  one  can  rewrite  (5.48)  as 


Jn, 


^12  ~ 


/ 


H^dCjb) 


sinvi^j  sinv(^2 


0  Hv(k2b)+Z^(0)H^(k2b) 


cosv(j»2^  dv  . 


(5.49) 


Letting  a=/T^k2a,  S  =  '^Crk2b,  and  y  =  k2b,  one  may  write 


H  (kb) 

V  2 
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H;  '(3)  - 


(5.50 


Noting  that  Q(v)  is  a  slowly  varying  function,  the  contribution  to  the 
value  of  the  integral  of  (5.49)  due  to  large  values  of  v  will  be 
negligible  because  the  remaining  portion  of  the  integrand,  i.e., 
sinv(|)i  sinvi»2 

V  V  cosv4)2^,  is  highly  oscillatory  (and  diminishing)  for 

large  v  so  that  the  positive  and  negative  parts  of  the  entire  integrand 
tend  to  cancel  each  other.  It  is  then  clear  that  major  contribution  to 
the  value  of  the  integral  arises  from  small  values  of  v  only.  In 
Appendix  J,  the  Debye's  asymptotic  approximations  of  Hankel  functions 
are  employed  to  simplify  Q(v)  yielding 


•k«sinF( v) 


(v/a)^  cosF(v)  -j/k^  _  (v/a)2  sinF(v) 


"1  “1 

+  v(cos  v/a  -  cos  "v/S)  . 


(5.5 


Introducing  new  variable  t  =  v/a  and  using  the  thin  substrate 
approximations  that  bi>^  s  t=l,2,  b^  “ 

rewrite  (5.51)  as 

ko%  „  sinF(at) 


•12 


2  /  2 
C  sinF(at)  +  j 


-  c  cosF(at) 


sinaitj^C  sinaii^C 


a<t.2?  cosa<^2j.c  dc 


(5.53) 


Recognizing  that 


F(v)  =  b  ye^k2-( v/b)^  -  a  J + 


V  cos  I  "^*^2^ 


(b-a)  =  t 


(5.54) 


as  a  -*•  »  while  keeping  t=(b-a),  a(|»j^  a4>2  and  constant,  one  obtains 


tim 


a+<»  z,-  ~ 

IZ  n 


Z  d  « 


sin(t'/k2-t2) 


^  V  k^-c^  sin(t  +  j  v/k^-c^  cos{t  \/  k^-c'^) 


sin(a<|>^c)  sin(a((»2c) 


34.2? 


(5.55) 


which  is  the  mutual  impedance  between  two  infinite  strips  on  a  planar 
dielectric  slab  (grounded),  with  widths  2a4)^,  ^^zc  ^^istance 

apart,  just  as  one  might  have  expected.  (5-55)  is  the  final  form  used 
for  computing  z^^  when  small. 

Case  (b):  b4)2^  is  large 

As  just  commented,  zjp  (5.48)  is  most  appropriately  represented 
as  a  residue  series: 


'IZ 


^ 

P=1 


rrr  (n 

3v  [Hy  (y)+Z^(o)H^  ^y)] 


rj(v)+r2(v)-r3(v)-r^(v) 


V=V 


(5.56) 
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are  the  roots  of 


where  and  v^,  p=l,2,3 . 

+  Z^(o)hJ,^^y)  =  0  .  (5.57 

In  principle,  the  roots  of  (5.57)  can  be  located  by  numerical 
techniques.  However,  the  problem  of  finding  even  the  first  root  would 
seem  to  be  a  formidable  one.  It  is  then  most  desirable  to  be  able  to 
approximate  these  roots  analytically. 

In  applications,  the  thickness  t  of  the  dielectric  coating  is  much 
less  than  the  radius  a  of  the  cylinder.  Thus  it  is  reasonable  to 
assume 


|ub|  >>  1,  |ub-ua|  =  |ut|  <  1,  (5.58 

where  the  second  inequality  means  that  the  dielectric  subtrate  is 
electrically  thin.  Then,  as  shown  in  Appendix  K,  it  is  justifiable  to 
write 


which  is  independent  of  v.  In  particular, 
1  3 


So  it  is  evident  that  the  roots  of  (5.57)  will  be  reasonably 
approximated  by  the  roots  of 


f 


(5.61) 


where 

\  =  l/Zg  (5.62) 

is  considerably  less  than  1  for  thin  substrates.  Hence  one  would  expect 
the  roots  of  (5.61)  are  just  slightly  perturbed  from  the  zeros  of 

^!i^\y).  It  is  known  that  the  zeros  of  Hf.^^(Y)  or  (y)  when  v  -  y, 

|y|  large,  are  given  in  first  approximation  by  the  zeros  of  the 
appropriate  Airy  function  combinations  or  their  derivatives.  This  is 
evident  when  one  appeals  to  the  Watson  approximations  [26]: 


H 


(1) 


(y)  ~ 


jw^(a) 

m/T 


(1)’ 

Hy  (y)  ~ 


jw|(a) 

m^/iT 


where 


Y  1/3  kpb 

=  (2)  =  (“2“) 


(5.63) 


(5.64) 


(5.65) 


v-Y  v-k2b 

0  =  m  “  m  »  (5.66) 

wi(a)  =  (Bi(j)  +  jAi(a))  ,  (6.67) 

with  Ai  and  B-j  denoting  the  Airy  functions. 
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Substituting  (5.63)  and  (5.64)  into  (5.61)  leads  to 
Yg  , 

(  tn]  Wj(a)  -  Wj(a)  =  0  .  (5.68 

It  is  derived  in  Appendix  L  that  the  roots  of  (5.68),  when  (Yg/m)  is 
small,  can  be  expanded  into  fast  converging  series  given  below: 


J.  3  1  4  ‘P  6  ’‘P  6 

»p(d)  =  tp  .  d  .  3  d  *  ^  d  .  ^  d  ♦  yT  “ 


-.3 


•p  5 
r-  +  TE- 


d'. 


1,2,3  ... 


(5.69 


where 


d 


|Ys 

with  j-m 


«  1 


(5.70 


(5.71 


3nd  Tp  are  the  zeros  of  A^(-t).  The  values  of  ip  are  well  tabulated  in 
Logan  [29],  and  only  the  first  few  are  listed  below: 

=  2.3381 
^2  =  4.0879 
^3  =  5.5206 

^4  =  6.7867 

T5  =  7.9441 

T5  =  9.0227 


T7  =  10.0402 


It  is  clear  from  transformation  (5.66)  that  the  roots  of  (5.61)  are 
related  to  the  Op's  as  follows: 

»  P  "  1»2»3,  ...  (5.72 

Next  employing  the  Watson  approximations  (5.63)  and  (5.64)  in  0(v) 
defined  in  (5.50)  yields 

-jw^(a)/(m''iT) 

Q(^)  ; - 

jw^(a)/(m2/ff)  -  Zgjwj((j)/(m/ir) 

Wj(a) 

T  '  T  •  (5.7: 

m  [Wj(a)  -  ^  Wj(a)] 

Since  satisfies  the  Airy  differential  equation,  i.e., 

II 

Wj( a)  -  0Wj( a)  =0 

one  can  readily  verify  that 

3  •  '•  da  ayii(a) 

3v  Wj^(a)  =  Wj^(cf)  3^  =  JJj 

Using  (5.75),  and  the  fact  that 


'  _!!!,  i. 

”l(ap)  *  Y  ”l('^p^  "  d  '^l(‘^p^  »  P  =  •••  • 


one  can  compute 


'  i  f  '  >  i 

I  m  -  d  Wj(a)J 


wi(ap)  1 

iv  ■  <*•' 


or  0=0 
P  P 


It  follows  from  (5.77)  that 


Res  lq(»  ))  =  . -j - P  .  (5.7 

[»p  -  ^1 

Making  use  of  (5.78)  in  (5.48)  results  in  the  following  residue  series 
representation  of 


p=i  [Yd2]  [%  ^ 


*  “  '■  p  d'^J  >■  p  <^2  '  ^ 

■  (>'i(Vp)  +  r2(vp)  -  r3(vp)  -  r,(Vp))  .  (5.; 

Expression  (5.79)  is  the  final  form  used  for  computing 

Employing  (5.46),  (5.50)  in  (5.36),  and  making  use  of  identities 
(5.4)  and  (5.5),  one  can  easily  obtain  the  eigenfunction  expansion  of 
the  mutual  impedance  between  the  two  infinite  strips,  i.e.. 


sinn4)^  sinn(t>^ 


■12  =  2TTbi),  n 

i  ^  n=0 


n 


Expression  (5.80)  will  be  used  to  check  the  results  obtained  by  (5.55) 
or  (5.79). 
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D.  A  NUMERICAL  EXAMPLE 


An  example  of  mutual  coupling  between  two  infinite,  z-oriented 
strips  on  a  dielectric  coated  cylinder  is  presented  to  illustrate  the 
effectiveness  of  asymptotic  solutions.  The  antenna  geometry  and 
pertinent  parameters  under  consideration  are  specified  below: 
current  on  each  strip  -  uniform,  z-directed 
frequency  -  300  (MHz) 

Xq  (free  space  wavelength)  -  1.  (m) 

radius  of  conducting  cylinder  -  5.  (m) 

thickness  of  coating  -  0.005  (m) 

dielectric  constant  -  4. 

loss  tangent  -  0.0001 

width  of  each  strip  -  0.125  (m). 

The  mutual  impedance  ziz  is  a  measure  of  coupling  between  the  two 
strips.  It  is  computed  using  two  different  asymptotic  forms,  which  are 
the  integral  representation  in  (5.55)  and  the  residue  series  in  (5.79), 
and  the  exact  eigenfunction  expansion  in  (5.80)  as  well.  The  magnitude 
of  zi2  obtained  using  each  representation  is  plotted  against  the 
separation  S  between  the  strips  in  Figure  5.3. 

It  is  clear  from  Figure  5.3  that  the  coupling  between  strip  1  and 
strip  2  weakens  as  their  separation  increases.  The  integral 
representation  agrees  exceedingly  well  with  the  eigenfunction  expansion 


for  S<0.25Xq,  and  they  start  drifting  apart  for  S>0.75Xo,  with  the 
former  showing  greater  degree  of  coupling.  This  can  be  attributed  to 
the  fact  that  the  integral  representation  stands  for  the  coupling  in  the 
planar  dielectric  slab  case  in  which  no  energy  is  shed  into  the 
surrounding  space  by  the  surface  field  while  propagating  along  the 
air-dielectric  interface,  whereas  in  the  case  of  the  coated  cylinder 
energy  is  continuously  dissipated  into  the  surrounding  space  while  the 
surface  field  propagates  around  the  curved  surface  of  the  cylinder. 

Hence  coupling  becomes  gradually  weaker  as  strip  2  moves  deeper  into  the 
shadow  region  of  strip  1.  The  residue  series,  on  the  other  hand,  shows 
excellent  agreement  with  the  eigenfunction  expansion  for  S>0.75Xq.  The 
number  of  terms  required  in  the  series  decreases  as  S  increases.  When 
S=2Xq,  for  instance,  3  creeping  wave  modes  are  used,  and  only  1  creeping 
wave  mode  is  needed  when  S=3Xo.  In  the  other  extreme,  a  total  of  56 
modes  are  summed  for  S<0,25Xo.  It  is  quite  apparent  that  a  considerably 
large  number  of  modes  will  be  required  to  provide  a  good  approximation 
of  z\2  when  S=0.  It  may  be  noticed  that  the  overlapping  region  for 
these  two  asymptotic  representations  is  0.25Xo<S<0,75Xo. 
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CHAPTER  VI 


SUMMARY 


Mutual  coupling  analysis  is  important  in  the  design  of  antenna 
arrays  with  low  sidelobes,  or  where  other  tight  control  of  the  radiation 
pattern  is  required.  This  work  presents  an  analysis  of  the  mutual 
coupling  between  conformal  microstrip  antennas  for  the  following 
configurations: 

1)  rectangular  microstrip  patches  printed  on  a  grounded 
planar  dielectric  slab  (the  planar  slab  problem),  and 

2)  rectangular  microstrip  patches  printed  on  a  dielectric 
coated  cylinder  (the  coated  cylinder  problem). 

It  should  be  noted  that  the  method  of  analysis  proposed  here  can  be 
applied  to  microstrip  patches  of  arbitrary  shape;  although  rectangular 
patches  are  treated  for  the  sake  of  simplicity. 

In  analyzing  the  radiation  or  scattering  from  an  array,  one  can 
formulate  the  problem  as  an  integral  equation  in  terms  of  the  unknown 
currents  flowing  on  the  array  elements.  In  the  case  of  the  microstrip 
array,  this  integral  equation  would  contain  in  it's  kernel  the  grounded 
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dielectric  slab  Green's  function  (planar  slab  problem)  or  the  dielectric 
coated  cylinder  Green's  function  (coated  cylinder  problem).  The  moment 
method,  which  transforms  the  integral  equation  into  a  system  of 
simultaneous  equations,  is  the  technique  proposed  to  numerically  solve 
this  integral  equation.  In  setting  up  the  moment  method  solution,  a 
crucial  step  is  the  calculation  of  the  self  and  mutual  impedances 
between  the  expansion  and  test  modes.  The  central  theme  of  this  work  is 
the  efficient  computation  of  these  moment  method  quantities.  A  model 
has  been  provided  to  relate  the  moment  method  modal  impedances  to  the 
port  impedances  of  the  microstrip  antenna  array  under  consideration. 

For  the  planar  slab  problem,  a  moment  method  solution  for 
microstrip  antennas  using  the  rigorous  grounded  dielectric  slab  Green's 
function  has  been  presented.  Piecewise  sinusoids  are  chosen  as 
expansion  and  test  modes.  Good  agreement  between  calculated  and 
measured  values  of  mutual  coupling  has  been  shown.  An  expression  for 
the  far  field  radiation  pattern  has  been  derived  via  the  method  of 
stationary  phase. 

It  has  been  demonstrated  that  the  coated  cylinder  problem  can  be 
solved  using  two  different  approaches.  The  first  approach  utilizes  the 
Green's  function  restricted  to  the  case  where  the  currents  are 
tangential  to  and  also  situated  at  the  surface  of  the  dielectric  coated 
cylinder.  The  fact  that  the  source  is  being  introduced  only  through  the 
boundary  conditions  simplifies  the  calculation.  This  is  similar  to  the 
method  of  solution  for  the  planar  slab  problem.  The  second  approach  is 
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based  on  the  dyadic  Green's  function  for  the  dielectric  coated  cylinder. 
The  use  of  the  Green's  dyadic  enables  one  to  solve  a  larger  class  of 
problems  since  the  source  can  be  arbitrary,  and  is  not  restricted  to  the 
surface  of  the  coated  cylinder.  For  this  reason  the  Green's  dyadic  is 
employed  in  this  mutual  coupling  analysis,  with  the  idea  that  it  may  be 
extended  to  a  more  general  configuration  later. 

The  dyadic  Green's  function  for  the  coated  cylinder  has  been 
constructed  using  the  principle  of  scattering  superposition  which 
requires  the  complete  expansion  of  the  free  space  dyadic  Green's 
function.  The  structure  of  the  free  space  Green's  dyadic  has  been 
characterized  for  a  general  orthogonal  curvilinear  system.  It  is  shown 
that  the  Green's  dyadic  can  be  cast  in  a  form  that  consists  of  two 
solenoidal  components,  and  an  irrotational  component  which  has  the 
simple  form  of  a  dyadic  delta  function.  This  compact  form  permits  one 
to  obtain  the  eigenfunction  expansion  of  the  Green's  dyadic  most 
easily,  since  only  the  solenoidal  vector  wave  functions  are  involved  in 
the  expansion.  This  expansion  process  is  further  facilitated  by  a  new 
set  of  orthogonality  conditions  which  do  not  involve  the  p  coordinate. 

Very  general  expressions  for  the  elements  of  the  impedance  matrix 
and  voltage  vector  (in  the  moment  method  solution)  have  been  presented. 
Self  and  mutual  impedances  of  an  expansion  mode  on  dielectric  coated 
cylinders  of  various  radii  have  been  calculated.  The  results  converge 
to  that  of  the  Infinite  planar  slab  case  as  the  radius  of  cylinder 
increases,  just  as  one  expects.  Mutual  impedances,  between  two 
expansion  modes,  on  coated  cylinders  of  various  radii  have  also  been 
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computed.  As  expected,  the  magnitude  of  the  mutual  Impedance  decreases 
as  the  separation  between  the  modes  increases,  or  as  the  cylinder  radius 
decreases.  Expressions  for  the  far  field  radiation  pattern  have  been 
obtained  using  the  steepest  descent  integration  technique. 

All  of  the  above  calculations  for  the  coated  cylinder  problem  are 
based  upon  eigenfunction  expansions  derived  using  the  dyadic  Green's 
function  which  is  of  the  radially  propagating  type.  These  eigenfunction 
expansions  are  poorly  convergent  when  the  coated  cylinder  radius  becomes 
electrically  large.  For  large  cylinders  asymptotic  high  frequency 
solutions  are  preferable,  because  they  are  computationally  more 
efficient,  and  can  be  interpreted  in  terms  of  ray-optics.  An  essential 
step  in  deriving  asymptotic  solutions  for  the  coated  cylinder  problem 
has  been  the  conversion  of  the  radially  propagating  Green’s  dyadic  into 
a  circumferentially  propagating  type  by  using  the  Poisson  summation 
formula.  Thereby  alternative  expressions  for  the  mutual  impedance 
between  two  expansion  modes  have  been  obtained.  These  include  the 
residue  series  which  lends  itself  to  ray-optical  interpretation,  and  the 
integral  representation  which  is  in  a  form  more  suitable  for  asymptotic 
evaluation.  A  numerical  example  of  coupling  between  two  infinite 
conducting  strips  on  a  cylinder  with  thin  dielectric  coating  has  been 
presented.  This  example  serves  to  illustrate  the  validity  and  power  of 
asymptotic  solutions. 
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This  work  presents  an  effective  approach  for  analyzing  coupling 
between  conformal  microstrip  antennas  by  combining  the  Green's  function 
theory  and  the  moment  method.  Only  a  narrow  aspect  of  the  coupling 
problem  has  been  considered.  It  is  by  no  means  complete.  The  following 
areas  are  recommended  for  future  investigation: 


1.  Extending  the  aysmptotic  solution  of  the  mutual  impedance 
between  infinite  strips  to  the  case  of  finite  length 
patches. 


2.  Rigorous  modeling  of  source  excitation  (which  requires 
the  Green's  function  corresponding  to  source  inside 
the  dielectric). 


3.  Generation  of  data  on  microstrip  antenna  characteristics 
such  as  input  impedance,  radiation  pattern  and  gain  for 
practical  design  purposes. 


APPENDIX  A 


DETERMINATION  OF  SPECTRAL  FUNCTIONS 


Expanding  Equations  (2.2)  and  (2.3)  in  rectangular  coordinates  and 
making  use  of  (2.4)  and  (2.5)  yields  the  following: 


In  region  1  (dielectric). 
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In  region  2  (air) 
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Enforcing  the  boundary  conditions  at  2=0  (dielectric-air  interface) 
specified  by  Equations  (2.8)  and  (2.9)  leads  to 
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Use  of  Equations  (A. 13),  (A.l)  and  (A. 7)  gives 
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Use  of  (A. 15),  (A. 4),  (A. 10)  and  (2.10)  gives 
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Use  of  (A. 16),  (A. 6),  (A. 11)  and  (2.10)  gives 
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(A. 20) 


Adding  (A.17)  multiplied  by  ky  to  (A.18)  multiplied  by  kx  yields 
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Adding  (A. 17)  multiplied  by  kx  to  (A. 18)  multiplied  by  (-ky)  yields 
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Adding  (A. 19)  multiplied  by  kx  to  (A. 20)  multiplied  by  (-ky)  yields 
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Adding  (A. 19)  multiplied  by  ky  to  (A, 20)  multiplied  by  kx  yields 
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Using  (A. 21)  in  (A. 23)  leads  to 
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where 


APPENDIX  B 


MUTUAL  IMPEDANCE  BETWEEN  TWO  EXPANSION  (DIPOLE)  MODES 

In  this  appendix,  the  exact  expressions  for  the  tnutual  impedance 
between  two  surface  expansion  dipole  modes  on  a  lossy  grounded 
dielectric  slab  are  presented.  The  mutual  impedance  between  an 
expansion  dipole  mode  and  an  impressed  source  current,  required  in 
(2.46),  will  also  be  presented. 

Figure  B.l  shows  two  dipole  modes  of  current  densities  Jm  and  Jn, 
located  on  the  surface  of  a  grounded  dielectric  slab  with  parameters  uq 
and  epEg.  The  ambient  medium  is  free  space  with  parameters  wo  and  gq. 
Mode  m  is  centered  with  respect  to  the  (x,y)  coordinate  system.  The 
center  of  mode  n  coincides  with  the  origin  of  the  (x',y')  system  which 
is  displaced  from  the  center  of  mode  m  by  a  position  vector  (xo,yo). 

The  x'-axis  is  at  an  angle  a  with  respect  to  the  x-axis. 
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The  mutual  impedance  between  modes  m  and  n  Is  given  by 
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where  the  integration  is  over  the  surface  of  mode  n  in  the  z*0  plane; 
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Use  of  (B.6),  (B.7),  (B.l'J)  and  (8.20)  in  (8.18)  leads  to 
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The  mutual  impedance  between  expansion  dipole  mode  Em  and  source 


current  Ji  can  be  computed  according  to  (2,46) 


V™  =  /  E^(r)  •  J,(r)dy 


Since  source  Ji  is  modeled  as  a  vertical  filament  of  constant  current 
inside  the  dielectric,  it  can  be  represented  as 
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(B.23) 


where  li  is  the  magnitude  of  the  feed  current  and  (xf.yf)  denotes  the 
filament  location  in  the  x-y  plane,  with  respect  to  the  center  of 
mode  m. 

Making  use  of  (2.13),  one  can  write 
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Inserting  (B.23)  and  {B.24)  into  {B.22)  and  carrying  out  the  volume 
integration  gives 
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In  both  (B.24)  and  (B.25) 
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APPENDIX  C 


EVALUATION  OF  AN  OSCILLATORY  INTEfiRAL  BY  THE 
STATIONARY  PHASE  METHOD 


This  appendix  presents  the  derivation  of  the  leading  term  of  the 
asymptotic  expansion  of  an  oscillatory  integral  of  the  form 

E(f)  =  //  dk^^dk^  (C.i) 

AAAmAAA  »  A 

where  f  =  xx+yy+zz;  k  =  xk  +yk  +zk  with  |kl  «  k  ;  and  f  is  a  slowly 

X  y  z  0 

varying  vector  function  of  k^  and  k^,  and  has  no  singularity  near  the 
stationary  points  of  k*r. 

Equation  (C,l)  is  a  2-0  Fourier  transform  which  can  be  evaluated 
asymptotically  by  the  method  of  stationary  phase.  The  method  of 
stationary  phase,  simply  stated,  is  based  on  the  fact  that,  as  |r|  », 

the  dominant  contribution  to  the  integral  in  (C.I)  arises  from  the 
neighborhood  of  those  values  of  kx  and  ky  which  make  the  phase  k*r 
stationary,  i.e.,  k*r  does  not  change  for  first-order  changes  in  kx  and 


ky.  For  those  values  of  kx  and  ky  that  make  the  phase  kT  vary  rapidly, 
the  integrand  oscillates  rapidly  between  equal  positive  and  negative 
values,  and  the  contribution  to  the  value  of  the  integral  is  small. 


First  it  is  noted  that  k. 


2  2 

and 


/  2  2  2 

k*f  =  r| 

|k^sin9cos<ti  +  k^sin9sin<|)  +  ^ 

/k  -k  -k  cos9 
r  0  X  y 

where  9  and  <t)  are  the  usual  angle  variables  in  spherical  coordinates. 
The  stationary  phase  points  are  determined  by  a  solution  of 


(C.3) 


=  0 
y 


Using  (C.2)  in  (C.3)  and  (C.4)  leads  to 


kzsinscosifi 

k  =  a 

X  cos9 

k^sinesiniji 

k  =  2 

y  cos9 


(C.4) 


(C.5) 

(C.6) 


Since  k^+k^  =  k^-k^,  it  is  seen  that 


2  sin^e  2 
*^z  cos^e  ~  *^o 


2 

k 

z 


and  hence 


k,  =  k  COS9 

z  0 


(C.7) 


Substituting  (C.7)  into  (C.B)  and  (C.6)  yields  the  stationary  values  of 
kx  and  ky  given  below; 


k  =  k  sin9cos<t) 
xo  0 


k  =  k  sinasind) 
yo  0 


(C.8) 


(C.9) 


provided  9  *  ti/2.  In  the  vicinity  of  the  stationary  phase  point 
(kxo»ky^)  k*r  can  be  approximated  by  the  first  few  terms  of  its  Taylor 


senes,  i.e., 


k*r  =  k*r 


(k  -k  )  TT-  +  (k  -k  )  — 
X  xo'  3k  y  yo'  3k 

^  J 


1  (*^xOf*^yo) 

1 

3  3  “1 

2 

^2 

(k  -k  )  +  (k  -k  )  .. 

'  x  xo  3k  '  y  yo  3k 

^  V 

k*r 

— 

k*r 


(•^xo»*^yo) 


(*^xo»*^yo)  •  (C.IO) 


(It  should  be  noted  that  only  the  partial  derivatives  in  (C.IO)  are 
evaluated  at  (kxo.kyo).) 
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It  can  be  shown  that 


_  1  kx 

3k2  kT  =  -  r  COS9  +  ^3) 

^  (k  ,k  )  ^  ^  (k  ,k  ) 

xo  yo  xo  yo 


k^cos^e  (cos20+sin20cos2(|,) 


32  -  _  1  ky 

3k2  k-r  =  -  r  cose  +  ^3) 

y  (k  ,k  )  ^  ^  (k  ,k  ) 

xo  yo  xo  yo 


■  kQCOs29  [cos29+sin2osin2^) 


ak  ak 


=  -  r  COS9  ^,3 


(k  ,k  ) 
xo  yo 


stn2es1n*cost 


(k  ,k  ) 
xo  yo 


Consequently  the  phase  can  be  written  as 


kT  =  k^r  -  A(k  -k  )2-B(k  -k  )2-C(k  -k  )(k  -k  ) 
0  X  xo  y  yo  X  xo  y  yo 


where 


(C.ll) 
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X 


) 


B  =  - 


32  -  - 

kT 


3k2 

y 


(k  ,k  ) 
xo  yo 


> 


3k  3k 
X  y 


32  -  _ 

kT 


In  the  neighborhood  of  the  stationary  phase  point  f(kx,ky)  is 
slowly  varying,  so  it  may  be  replaced  by  f(kxo.kyo).  Hence  the 
asymptotic  value  of  the  integral  for  E{T)  is  given  by 


■Jl'o''  1 


4ii2 


X  y 


£ 


where  Rq  is  a  small  region  centered  on  the  stationary  phase  point  from 
which  comes  the  major  contribution  to  the  integral.  As  noted 
previously,  for  large  r,  the  domain  of  integration  may  be  extended  over 


the  whole  (k^.ky)  plane  without  introducing  significant  error  because  of 
the  rapid  oscillation  of  the  integrand  for  nonzero  values  of  (kx-kxo) 
and  (ky-kyo).  Thus 

E(?)  -  ^  ll 

X  y 

(C.12) 

Next,  applying  a  linear  transformation  T  defined  by 


k 

=  T 

s 

S//A  +  kxo 

t 

t//W  +  kyo 

to  (C.12)  yields 


E(F) 


-jkor 


1  *  j[s2«2-fCst/AB] 

4.2  //  e 


Idet  Jjldsdt 


~  ffk 


JL  “  j[s2+t2+Cst/AB]  d^ 


3kx 

3kx 

1 

3s 

3t 

/A 

0 

3ky 

3ky 

1 

3S 

3t 

0 

/B 

__ 

is  the  Jacobian  matrix  of  T 


Completing  the  square  for  the  s  variable  in  (C.13)  to  obtain 


-jk  r  0^  2  Qi 

e  1  «  j[(s  +  777^)  +  t2[l  -  -m)] 


Using  the  result 


“  JMu-uja  rj- 
/  e  ■  v|l  ® 


in  (C.14)  gives 


Efr)  ~  f(k,o.kyo)  4,2  •  /AB  ’ 


/T  e 


jit/A 


e  1 


““  ^^^*'xo’'^yo^  2n  ''4AB-C2 


1  -  4AB 


It  may  be  observed  that 


y4AB-C  =  j  4  .  4|^2(.os‘*0  (cos20+sin20cos2(^)(cos29+sin20sin2(|») 
I  0 


■  k2cos‘*0 
0 

r  ~  ~ 

"  k  cos20  cos‘*0+cos20sin20(sin2(j)+cos2^) 
0  I 


kpCosS 


Therefore,  the  leading  term  in  the  asymptotic  expansion  of  the  integral 
in  (C.l)  is  given  by 


jkQcos0  e 


T“  .k  ) 

r  '  xo  yo' 


(C.15 


APPENDIX  D 


ORTHOGONALITY  PROPERTIES  OF  CYLINDRICAL  VECTOR  EIGENMODES 


Several  orthogonality  properties  of  the  solenoidal  vector 
eigenmodes  or  wave  functions  in  cylindrical  coordinates  are  stated  and 
proved  in  the  form  of  lemmas.  For  the  sake  of  clarity,  shortand 
notations  are  introduced  first: 


1.  Sn  designates  a  "+"  or  sign,  n=0,l,2,... 


(n.i) 


2. 


^mn 


1  m=n 
0  m*n 


,  m,n  *  0,1,2,.... 


(0.2) 


3. 


1  n=0 

2  n>l 


(D.3) 


4.  X'  =  yk2  -  h'^  ,  ImX'  >  0,  -«<h'<» 

X  =  yk2  -  h2  ,  ImX  >  0,  -«»<h<««» 


(D.4) 
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•  r"T*  V- 


^  mi  ,•  y»r.  J  .  u  '■  J< 


■:  '"j  '-f  L-- 


"  'T*  .*  ■  «ru 


r- 


I  i 


i  r- 


■ 


^em 
0  0 


Lemma  1 


A  r  +  A  r  ,  n>0 
en  en  on  on 

,  n=0 

_  en  en  * 


(D.5) 


/  dh'  /  dz  /  pd(t.  [Ae^  p  x 


-OD  -M  0 

An^ 


0  0 


—  Ae_  r.  px2  (Xp)  7-  Z^"'(Xp)  6 
^n  0”  0"^  n  3p  m  mn 


Proof; 


Substituting  (3.46)  into  the  integrand  furnishes 


L.H.S.  =  /  dh'  /  dz  /  pd^  [Aen{4(x2Z^"(Xp)^"^  nt)Tz(^)}] 


=  /  dh'  /  dz  /  pd(|)([Ae„  5°^  5°n  "i*] 

.CO  -<n  0  0  0 


?13 


( 


^  *-•  *-•  V  '-•  **- 


■  rf..*  ■■  iV  . 


One  can  observe  that 


/  d*  [Ae„  "♦][■■=„  IT„ 
0  0  0 


/  cosn(()Costn(|)  +  A  r  cosni^sinm(|)  +  A  r  sinn(^costnd) 

Q  en  em  ^  ^  en  om  ^  on  em  ^  ^ 


+  A  r  sinn^sinm<|i] 
on  om  •' 


2^  A^„r 
en  en 


n=m=0 


it[a  r  +  A  r  ]  n*m^0 

^  en  en  on  on-* 


?ii 

—  Ap^  rp„  6 
Cn  en  en,  mn 


/  dz  =  2iT5(h+h')  . 


(0.6) 


Thus, 


“  4ir2  Sn  3  Sm 

L.H.S.  .  J  dh'  1  —  »e„''e„P>'Z„  (>»)  T;  Z„ 
-at  n  0  0  ^ 


3p  m 


,  Sp  3  Sm 

^»e„re„px2z„  (Xp)  “  Z„  (Xp)J, 


n  0  0 


V-VvlV.T'V'^-  V-'.  .■■■ 


Lemma  2. 
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Making  use  of  (D.6)  it  is  clear  that  integration  on  z  produces  the  delta 
function  2Tr5(h+h').  Hence, 


Clearly  L 


P.TT 

!  cif 
0 


Simi  1  arly 


2n 

I 


=  2n  /  dh'6(h+h‘)  /  pd^  {[-Ae„  n*][?re„  n*] 

-«  0  0  0 


S„  Sm  ,  , 

k.»  z„ 


*  cos  "♦ltre„  ,To 

0  0 


cos 

sin  k„p 


j(X*)2hn  Sn  Sm 

-  Z  "(Xp)  Z^"’(X-p)  . 

0-  " 


H.S.  =  0  for  n=m=0,  and  n#m.  When  n=m^O, 


[-A 


cos 

®n  sin 
0 


"♦ll'i'en  cos 
0 


77T 

I 


0 


[^onf^^cosn-iksinn*  -  A„„r„„cosntcosn(|>  +  A^„r^„sinn(|)sinn(|) 
en  en  ^  en  on  on  en  ^ 


~  sinn(|icosn(|)]  =  u[-A  r  +  A  r  1 

on  on  ■'  '•  en  on  on  en-* 


it  can  be  verified  that 


[  +  Af 


sin 
n  cos 


ni.][r 


cos  s 
Sn  sin 


n[-A  r  +  A  r 
'■  en  on  on  en- 


+ 


I  ^’''(-Aonr  +A  r  )  /  dh'6(h+h')[' 
'  en  on  on  en  •'  '  '■ 


L.H.S.  = 


Sn  Sn 

.  Z^"(Xp)Z^"(X*p) 


n=ni*0 


=  0  for  al 1  m,n 


Lemma  4. 


j  dh'  /  dz  i  pd<t.  [Aep  P  x  7  x  ]  •  [^en,  Me^x' ^ 

-00  .00  0  0  0  0  0 


Proof: 

After  substituting  (3.46)  and  (3.47)  in  the  integrand  one  will  see 


00  00  2ir 

L.H.S.  =  /  dh'  /  clz  /  oM  {Aen[-z(*)]h{re^[?p(*)4(*)]}  =  0  . 


APPENDIX  E 


PROOF  OF  EQUATION  (3.73) 

The  proof  starts  from  Equations  (3.69)  and  (3.71)  which  are 
relabelled  here  for  easy  reference: 


7x7xGg^(F,F')  -  k^G^^(F,F')  =  0  .  (E.l) 


7x7xG^l(r,F' )  -  kjG^j(r,F')  =  0  ,  (E.^) 

2  2  2 

where  r  e  V^,  r'  t  'Jj  =  to  =  to  u^e^e^.  Taking  the  curl  of 

(E.l)  and  then  subtracting  it  from  (E.2)  leads  to  the  conclusion  that 

"  '7xGg^(r,r')  +  t|»(r,r')  (E.3) 

where  tp  satisfies  the  same  differential  equation  governing  and  G^^. 
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... , 1^1  ,..|,.I .  .  pi.  »  Ml ,  P'  ■  J  ■  j'-  '."■  H'»  1' ■  M  ■■  1  ■  1."  1  ■  ».■  »■  1  .■  »■ 

"T" 

t: 

ft. 

fF- 

i 

Next  letting  and  be  the  electric  and  magnetic  field  in 

“i- 

•  /r 

respectively,  due  to  electric  current  source  J  in  (with 

parameters  \i2  ^nd  e, ),  it  can  be  shown  that  [19] 

•  . 

■ 

E,(r)  =  ju)U2  /  G  (r,r')»J(r‘)dv* 

(E.4) 

^ . 

• ' 

From  Maxwell's  equations,  it  is  clear  that 

-  .  '7xE^(r)  H2  = 

"l"'*  '  j-u,  %.  1  7.G^,(r.r’).0(f)dv'  . 

(E.5) 

■ 

However,  in  order  to  be  consistent  with  the  role  of  the  magnetic 

Green’s 

•k 

dyadic  defined  by  (3.18),  G^i  must  be  sought  such  that 

• 

Hi{r)  =  /  G^i(F,F').j(F')dv' 

''2 

(E.6) 

1 

Substituting  (E.3)  into  (E.6)  yields 

ft 

Hj(r)  =  /  [7xG^^(r,r')+ij;(r,r')]*3(r')dv' 

''2 

{E.7) 

_ . 
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— 

*■  *  .*  ’  ■ 
L  ,  .'1‘ 

V-  V-  0-  n  -  - 

From  (E.S)  and  {E.7)  it  is  seen  that 


'l'(r.r')  =  (—  -  l)vxGg^{r,r*) 


(E.S) 


Using  (E.S)  in  (E.3)  gives 


''-ml"'-'"* 


(E.q) 


(D  II 


APPENDIX  F 


Likewise,  from  (F.l)  and  (3.89)  it  is  clear  that 


=  Gg^(r,F') 


p=b"  or  p=a 


p=b  or  p*a 


0  0  e  e  0 


e  e  0  p=b  or  p-a 


(F.3) 


These  24  unknown  coefficients  (Ae*,  Be-,  Cg-,  Den*  “ten*  ^^n*  ^en*  ‘^en* 

0  0  0  0  0  0  0  0 
ae  ,  be  ,  de  ^en^  determined  from  the  24  linear  equations 

o"  0^  0^  0^ 

derived  from  the  following  set  of  boundary  conditions  (as  defined  by 
Equations  (3.82)  to  (3.84)): 


p  X  Gel^^’^'^1  ^  ° 

P=a 


(F.4) 


P  X  Sg2(r,r')  "  ^  ^  Ggj(r,r') 

p=b  P=b 


(F.5) 
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•V/.. 


I  af: 


p  X  V  X  G^^(r,r') 


p=b* 


=  ;x^7xG^^(r.r-) 


p=b‘ 


(F.6) 


Furthermore,  it  may  be  helpful  to  recall  the  explicit  expressions  of  the 
and  N  functions  which  are  given  below  for  the  sake  of  convenience. 


Menx 

0 


(h.r) 


T  *  -  \  sin  ^ 

^  P  P  cos 


£ 

3p 


\  cos 

♦  Z^(^p)  s^n  n4. 


Jhz 


(F.7) 


Nenx(^r) 

0 


PIT  sjZ„(ip) 


cos 

sin 


*  jhn 
^  ~Tp 


l-„M 


sin 

cos 


nij) 


*  +  co< 


Jhz 


(F.8) 


where  Z]^(Xp)  =  H^^^(Ap),  Z”(Xp)  =  J  (Xp),  and  X  ,  ImX>0.  For 

simplicity,  H^^^(Xp)  will  be  denoted  by  H^(Xp),  and  ^F{?)  by  f'(5) 
throughout  this  appendix. 

Substituting  (F.3)  into  (F.4)  produces,  for  n=0,l,2 . 


-+ 


a  a  A  A  A  A 


J 


Using  (F,7)  and  (F.8)  in  (F.9)  and  in  (F.IO)  yields,  respectively 


zl-u[penH„(ua)  *  Te„J„(iia)l  *  ^7  [8o„H„(lia)  *  t,  J  (i,a)]}“ 

0  0  1  g  e  ^ ' 

1  e  e 


A  1  I  inh 

2l-p[bo„H„(ua)  *  fo„J„(pa)]  *  77  [ae„H^(ua)  +  be„J„(ua)]}'' 
e  6  1  0  0 


**  '-k^tynfbb'^yntPb)]”^ 


Next  substituting  (F.2)  and  (F.3)  into  (F.5)  gives,  for  n=0,l,2,. 


0  00  e  e  p=b 


P  [‘^n”enp(h,^)+SonNJny(h,r)+YenM;  (h,r)+TonNo  (h,r)] 


p  X  [Ne„;^(h.F)+CenNenx(^‘*^)+Vonx(^’'^)l  ,+  = 

0  0  0  e  e  p=b 

P  [aenNJnu(^*'')+^OnM5np^^»^^‘*’^enNinj,(h,r)+fonM5  (h.r)] 

Using  (F,7)  and  (F.8)  in  (F.13)  and  in  (F,14)  yields,  respectively 

j  nh  co^  ^ 

and 

2  l-0o„2H^(xb)  *  ^  [0„(xb)rte„H„(xb)])®^"  nt 

6  2  Q 

♦  ♦|-[Jn(xb)tCe„H^(ib))  b»  =  M-u[b„„H^(ub)+fo„J^(gb)J  J 

0  2  e  e 

2 

•  [aenH^(yb)  +  <lenJn^P^^  1  ^cos  ^ ^s 

Finally,  substituting  (F.2)  and  (F.3)  into  (F.6)  gives,  for 

n=0,l,2,.,., 

P  X  ^{7xMen;^(h,r)  +  Aen7xMenj^(h.F)  +  Bo„7xNj^j^(h,F) } 

2  0  0  0  6  6  p=b 

"  ^  Son^xN^^  (h,F)  +  Yen'^xMe^  (h,F) 

loo  e  e  00 

+  Ton7xNon^^(h,F)}  _ 


and 


P  X  ^{7xNenx(h,0  +  Cen7xNenx(h.r)  +  . 

2  0  0  0  e  e  p=b 

=  P  X  ^{aenVxNepy(h,F)  +  bon7xMony(^r)  +  denVxNef,y(h.F) 
loo  e  e  00 

+  fon^xMon  (h,F)}  .  (F.18) 

e  e  p=b 

Employing  (F.7)  and  (F.8)  in  (F.17)  results  in 

'‘Ml  •  jnh  sin 

2  -Z  (-k2XB„„H„(Xb)  J  —  lJ„(Xt>)*Ae„H„(Xb)ll^„5  H 

e  0 

+  *  ■;;7{-^MJn(>^b)+AenH^(Xb)]}^®^  n^.  =  z{-kiii[eonHn(Mb)+Tg^J^(Mb)] 

jnh  51 n 

+  i{-MM'*enHn(Mb)+YenJn(''‘’^^^s?n  (''•19) 


and 


^  u  1  •  •  jnh  ,  ro!^ 

z  ■;:7  {-k2X[J„(Xb)+CenH^(Xb)]  ±  ^  DonHn( Xb) n^ 

"  z{-'^iM[aenHn(Mb)+denJn(Mb)] 

jbh  1  ,cos 

-  b  ['’0nHn(Mb)+Wn{Mb)]}5^„  n* 

G  G 

+  <fr{-‘'^[^OnHn(Mb)+fonJ„(Mb)]}^Q5  n^  .  (F.?0) 

G  G 


^.: 


%U. 


!r 
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For  n=0,l,2 .  (F.ll)  must  hold  for  all  ^  and  z,  thus  one  can 

that 

SonH  (ua)  +  To  J  (wa)  =  0 
e  e 

oenH  (ua)  +  Ye^J  (ua)  =  0 
0  0 

Using  the  same  argument  one  can  see  from  (F.15)  that 


*  Ye„uO^(bb)  5  („b)  t  („b) 

0  0  el  el 

-  Ae„XHjxb)  t  Bo„^H^(Xb)  =  x/{lb) 

^  "„(“»)  *  'o„  ^  J„(Mb)  -  Bo„  ^  K^(Xb)  =  0  . 

Similarly,  (F,19)  implies 

jnh  jnh  • 

“^n  b  ^®n  b  "  ‘^On'^ii'J 

0  0  "  e  "  e 

^'i  jnh  ui  I  Ui  jnh 

*  ''bni;; “%(''>)  +  Bon^tx'HJxb)  .  >--5-J„(xb) 

<>e„u2H^(ub)  +  Ye„j2j^(ub)  -  ^  x2H^(xb)  •  ^  x2j„(xb) 


conclude 

(F.21) 

(F.22) 

(F.23) 

{F.24) 

(ub) 

9 

(F.25) 

•  {F.26) 
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A  second  group  of  equations  corresponding  to  (F.21)  -  (F.26)  can  be 
similarly  obtained  from  (F.12),  (F.16)  and  (F.20).  For  each 
n=0,l,2,...,  and  for  all  41,  Equation  (F,12)  implies 


aenHn(ua) 

0 

0 

(F.27) 

boX(“») 

1 

+  fonJ^(ua)  =  0 

• 

(F.28) 

Equation  (F.16)  implies 


jnh  '  jnh  < 

*  k  b  ^  ic,b 

0  1  e  0  1  e 


jnh  '  jnh 

-Do,XH^(Xb)  «^]^J„(Xb) 


(F.29) 


«en  ki  +  de,  J^(ub)  -  Ce„  H^(xb)  =  J^(xb) 


(F.30) 


and  (F.20)  leads  to 


'  jnh  '  jnh 

-  aepk^pH  (wb)  ±  bg^  b  H  (ub)  _  dg^k^uJ  (ub)  ±  fgp  b  J-(ub) 

11  ^  i»  A  If  ^  II 


^1  '  jnh  ui  ' 

*  cen  ’  “o„  iT— HJ'M  “  -;::k2u  (xb) 


(F.31) 


U2  b  ^  (J  [ae^H^(wb)  + 


One  can  also  rearrange  (F.25)  such  that 

Bo„  ;^k2»Hr^(Xb)  »  *  [Vn'*-*’'  *  ''n'''”'!  * '’T 


+  YenJn(Mb)]+  kiu[BonH^(pb)  +  To^J^(ub)] 


Combining  (F.38)  and  (F.39)  together  leads  to 

1  ^2  p  ^  I  ^  ^  I  ^  p2 

"  H;(Xb)  *  HI  X  J  (l? 


'2  ''  e 


Equating  (F,36)  and  (F.37),  then  using  (F.33)  and  (F.34)  yield 


1  U2  2  J^(iia) 

";;7(x)  ^n^^^)her^] 

1  '  .  £  '  '  jn^  1^  u^ 

"  H;(xb)  +  X  “  H;(ua)  *  k^b  X  ^2 

Op(ua) 


which  can  further  be  rearranged  to  get 


E 


^  ■  Jn(i‘b)H^(Xb)}  =  J^(Xb)H^(Xb)  -  J^(Xb)H^(Xb) 


jnh  1 

-  kjb  X  ■  H^(ua)  •  (F.4] 


Making  use  of  the  observation  that 


J„(XMH„(Xb)  -  J„(Xb)H^(Xb)  .—  . 


Equation  (F.41)  can  be  used  to  express  Ye^  follows: 

0 


j2  1 

”  •“!>  *0  *  *0  > 


(F.4? 


(*here 


''o  ■  -  J„(-b)H_^(Xb)] 


JWua) 

_  ^2  u 


iilM  t«  I  -  \(wb)H„(Xb)] 


(F.43 


jnh  1  u  2  J^(w3) 

=  k,b  M  -77;rj-H^(ub)]H^(Xb) 


(F.44 


Next  equating  (F.35)  and  (F.40)  and  rearranging  gives 


2 

k2  m2 


[eonHn(>^b)H„(Xb)  +  Tof,J^( ub)H^( Xb)  ] 


e 

U2 


1^1 


{uki[6onH^(ub)H^(Xb)  +  To„J„(ub)H^(Xb)] 


jnh  kiu  M  ^ 

"Tm  1(1)  -  llK„H„(iib)H^(xb)  *  Ye„J„(»b)H„(xb)ll  , 


(F.45) 


which  can  further  be  manipulated  by  using  (F.33)  and  (F.34)  to  give 


J^(pa) 


k2  2 


1*2  I 


iu;:!)  1(^7)  7“„(>‘W„(xb)  .;7H„(„b)H„(xb)i, 


k  -p  2  p 


J'{Ma) 
n 


^2  jnh  M  ^ 


(F.46) 


or 


(’^2-X3)Ton  =  ^  ^4^en 
e  0 


(F.47) 


where 


Ic  2 

■  KkT)  I  Jn< --TT  V 1 


( 


X3  =  [(]^)  I  H^(ub)H^(Xb)  -  -  H^(ub)H^(Xb)] 


‘^2  jnh  1  u  2 

"  yi  kib  U  ■  H;;(ua) 


Finally,  using  (F.42)  in  (F.47)  yields 


j2  1  ^ 

lX2-K3ho„  =  '  xa.b^  Xo  -  Xo  >'  “'■ 


j2  H* 


^On  "  ^  ffbu  Xo 


X^X, 


X2  -  X3  + 


Up  to  this  point,  coefficients  og^,  Sg^,  Yer,»  t 

000 

been  determined.  The  remaining  set  of  undertermined 

consists  of  ag^,  bg^,  dg^^  fg^,  Cg  and  Dg  and  will 
00000  0 

similar  fashion  as  the  previous  set. 


en,  Aen  af’d  B 
0  0 

coefficients 

be  found  in  a 


It  is  clear  that  (F.27)  and  (F.28)  give 

Jn(  v»3i) 

’0"  '  ■  ’’S" 


Jn( Ma) 

'  ■  H:(lia)  '2n 


OIO 


From  (F.32),  one  will  find 


'  *  Wo''*'’’) 

e 


On  ■  H^{xb)  ui  ^"on"n 


J  (ua) 
n 


H^(Xb)  Ui  ^^x^  ■  H;(ua) 


and  from  (F.29) 


'  fb  *  <:e„”„(‘''>l  ■  *  Tb  IVn'*'^’  *  V"'"*'" 

2  0  10  O 


-  a 


[bon^(^b)  + 


Rut  one  gets  from  (F,30)  that 

inh  jnh  vi 

n  (p  [ae„H„(»'>)  *  «eA(bl>)] 

2  o  10  0 


Employing  (F.56)  in  (F.55)  results  in 


=  H;(Xb)  X 


{F.54) 


(F.55) 


(F.56) 


(F.57) 
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Next,  equating  {F.54)  and  (F.67)  and  making  use  of  (F.52)  and  (F.53) 
produce 


1  U2  u  2 


J'  (ua) 


H  (Xb)  -  HVua) 


1  u  . 

hIW  X 


.  .  1  2  J-(ua) 

inh  1  u  ^  f' 

*l^u  1(I>  - 


which  can  be  rewritten  as 


°  ’  ’‘o  H" 


One  can  readily  deduce  from  (F.56)  and  (F,52)  that 


'<2  ,1.,^ 


J  (ua) 
n 


^o"  "  ki  ^X^  '  Hn(ua) 


Also,  one  can  obtain  from  (F.31)  that 


jnh  1 


U2  ki 


Can  ^'^2*H„(^b)  .-k,*l.J^(xb)  *— —  D„„H^(xb) 


I  I  U2  jnh  1 

•  [ae„H„(ub)  ♦  dy^(Bb)]  *  —  7  *  Vn’"*’’!' 


I 


Making  use  of  (F.52),  (F.53)  and  (F.54)  will  reduce  (F.60)  to 


'fin 

0 


H;(Xb) 


^2  jnh  1  , 


{-J  (Xb)  ±  .  ,  . 

^  n  ui  k2b  X 


[(I)  - 1] 


•  [Jn{ub) 


J'(ua) 


U2 

U1  k2X 


•  [J^(ub) 


J^(Ma) 

H„(MaT 


Hn(wb)]de^}  . 


(F.61) 


Now  if  one  equates  (F.59)  and  (F.61),  makes  use  of  (F.58)  and  observes 
•  •  j2 

that  J  (Xb)H  (Xb)  -  J  (Xb)H  (Xb)  =  ,  he  will  find  that 

n'  'n  n  n  uXb 


iL  r  ^1^4  _j2 

k2  X  1^2  -  ^3  +  Xo  '  “  nXb 


or 


irb  ^ki  VJ 


X2  - 


^3 


(F.62) 


So  far  all  coefficients  have  been  determined.  However,  they  can  be 
simplified  significantly.  To  facilitate  the  simplification  process, 
it  is  convenient  to  define 


Pf,  = 


Qn  =  Jn(Mb)  +  qnH^(yb) 


^n  =  On(iib)  +  %\W 
Then  it  can  be  shown  that 

’‘o-  l-'ln  "  '■„T„lH„(Xb) 

jnh  1  M  ^ 

=  -llH„(Xb)S„ 


^2  I  k2  ^ 


^3  =  tui  Pn^(P^)  -  ^ki)  X 


^2  k2  2 


•  l-TT '■n  *  W  ICxSnlHnl^O) 


Hn(Xb)  2 

Xo(X2-X3)  +  =  [-;— J 


li2  u2  P2  1^2  2  y2 

,1  X  Vnl*t-ui  U J  X 


1^2  nh  2  n  2  2 

^  (m)  KD  -  1|  Vn  • 


u  — 

-  ,Tbu  (Xq(X2-X3)+XjX4) 

2nh  ^^2  u  2  Hf)(Xb)  Tp 

'  ^  rtkib2  (^)  [i\)  -  1]  i  y2  ]  [H^{Xb)/y]2Y 


2nh  M2  £  2  Tp 

=  Tikib2  (jJ  [l  -  (x)  ]  H  (Xb)Y^ 


TTby  (X(,(X2-X3)+X^X4) 


[H  (Xb)/M]2Y 


TTb  ^ki%2  [H^(Xb)/y]2 


0)0  o  o  c  o  o 


APPENDIX  G 


ON  THE  FUNCTION 


In  solving  electromagnetic  problems  that  involve  Fourier  integrals, 
one  will  encounter  a  multi-valued  function  of  the  form 

f(h)  =yk2-h2  (G.l) 

where  k  =  k'+jk"  is  the  propagation  constant  and  h  =  h'+jh"  is  the 
integration  variable  going  from  -•  to  •,  In  addition,  one  may  notice 
that  k'>0  and  k"><0  for  the  time  variation.  To  ensure  convergence 

of  the  integrals  one  often  finds  it  necessary  to  enforce  the  condition 
that 


Imf(h)  ><  0  (G.2) 

for  the  e'^j'^t  time  dependence.  It  turns  out  that  condition  (G,2)  will 
define  the  path  of  integration. 


Letting  f(h)  =  t'+jt",  it  is  apparent  that 


f2(h)  =  k2-h2  =  k'’--k"2-h'2+h''2+j2(k‘k"-h'h")  =  T'2-T"2+j2T’ t“  . 

(G.3) 

Equating  the  imaginary  parts  of  (G.3)  gives 

k'k"-h'h"  (G.4) 


Next,  equating  the  ral  parts  of  (G.3)  and  then  employing  (G.4)  yield 
k'k"-h'h"  ^ 

-  ]  =  k'2.k*‘2-h‘2+h"2  .  (G.5) 


which  can  be  readily  solved  to  obtain 


,,  k•2.k«2,^,'2+h'•2  2  2,1/2 

T-  =  {[( - -  )  +  (k'k"-h'h")  ] 


k'2.k"2.h'2+h''2  1/2 

- ; -  } 


(G.6) 


It  is  clear  that  t'  is  always  positive.  Making  use  of  this  fact  in 
conjunction  with  (G.4)  leads  to  the  conclusion  that 

1.  for  t">0  on  the  path  of  integration  (e-jwt  time  dependence), 
one  can  set  h"=0  when  k">0,  and  h'>^0.  h"<xO  when  k"=0: 


h-  PLANE 


Figure  G. 


Figure  G. 


.  Path  of  integration  r  in  the  h-plane  and  the  analytic 
properties  of^k^-h^  when  k"=0  (for  time  dependence) 


Imh 


Path  of  integration  r  in  the  h-plane  and  the  analytic 
properties  of  /k2-h2  when  k"=0  (for  ej“^  time  dependence). 


APPENDIX  H 


ON  THE  RATIOS  Pn/Sn  AND  On/Tn 


From  Equations  (3.11)  -(3.14)  and  (3.116)-(3.117) ,  it  is  apparent 

that 


Pn(i')  i>^(ub)H^(ua)-J^(ua)H^(ub) 


S  (li)  J  (ijb)H  (lia)-J  (ua)H  (ub) 
n  n  n  n  n 


(H.l 


and 


T„(,i) 


jjab)Hjja)-j'(ua)H'(ub) 
n  n  n  n 


J^(iib)fyua)-J^(ua)H^(ub) 


(H.; 


Since  H^(5)  =  =  Jf,(C)  +  where  Y^(5)  is  the  Neumann 

function  of  order  n,  one  can  rewrite  (H.l)  and  (H.2)  as  follows: 


Pn('^)  J^(ub)Y^(ua)-J^(ua)Y^(pb) 

'  S 

Sn(u)  J^(ub)Y^(;ja)-J^(ua)Y^(ub) 


(H.: 


ftl 


On(u)  J^(ub)Y^(iia)-J^(ua)Y^(ub) 


T^(ii)  J^{ub)Y^(iia)-J^(iia)Y^(ub) 


(H.4) 


By  examining  the  ascending  series  for  and  Y^  [2fi],  one  will  note  that 


(H.6) 


and  Y  can  be  put  in  the  form 
n 


Y^(2)  =  ^2)  +|lnf  J^(2) 


(H.6) 


where  Y(2)  has  the  property  that  Y(-2)  =  (-1)  Y(2).  Consequently, 
Y^(-2)  =  Y(-2)  +|ln  (^)  J^(-2)  =  (-1)" 


'  2  2 

Yn(z)  +  7  In  Y  d^(z)  ±  j2J^(2) 


(H.7) 


Formally  differentiating  (H,6),  (H.6)  and  (H,7)  yields 


J^(-2)  =  (-l)"'^j’(2) 


(H.8) 


Y,(z)  = 


- '  2  1  2  '  ^ 

Y„(z)  +7  I  J„(z)  +  In  7  J  (2) 


(H.8) 


n-1  2(1,  2*1  '  I 

Y^(-2)  =  (-1)  Y^(2)  +  7  I  7  (z)  +  In  2  Jn(z)  I  ^  J^JJz)  I  . 


It  follows  from  (H.5)-(H.10)  that 


Sn(-P) 


=  i  [Y„(ua)  +7  In^j  (pa) 

L 

±  j2J^(ua)]}  -  {(-l)"j^(ua)}  {(-1)"’^  [^ub) 


2  ua 


2  ub  1  1  1  ■ 

+  n  (In  2  J„(ub)  +  yb  ~ 


^  2d 

i  {(-l)"jn(ub)}  K-D"  LY^(ua)  +  7  In  1“  J^(wa)  ±  j2J^(ija)]} 


{(-l)"jri(’^a)}  •  ((-l)"[Y^(ub)  +  l  J„(pb)  ±j2J^(ub)]} 


(-1)^  ^{Jf,(‘'‘^)Yn(ua)  -  J  (ua)Y  (ub)}  P  (u) 


(-1)^"  {J„(ub)Y^(ua)  -  J„{ua)Y^(ub)}  S^(u) 


t(-i)"'V  („1,)|  l(-i)"‘'lyM) 


ua  '  '  n-l  ' 

+  In  T-  j,{ua))  ±  j2J^(ua)]}  -  {(-1)"  'j„{ua)} 


•  {(-1)”  ^(Vub)  +  n  Jn(ub)  +  In  2^  Jn(ub)) 


-  j2J^(Mb)]}  I  /  I  ((-l)"j„(ub)}  {(-1)"“My  (ua) 


n  (ua  J^(wa)  +  In  2  0^  (wa))  ±  j2J^(via)]} 

-  {(-1)'^  ^Jn(wa)l  •  K-D”  +  I  In  ^  vl^(pb) 

i  j2J„(ub)]} 

(-l)2(n  ^^(J^(ub)Y^(ua)  -  J^(wa)Y^{ub)}  0^(u) 

(-l)^"‘Mjn(wb)Y’(ua)  -  j'(ua)Y^(ub)}  ‘  T^(y) 

Thus  from  (H.ll)  and  (H.12),  it  is  seen  that  both  P  /S  and  Q  /T 

n  n 

odd  functions  of  u. 


APPENDIX  I 


CALCULATION  OF  CYLINDER  FUNCTIONS: 


I 


H 

n 


It  is  well-known  that  all  cylinder  functions  satisfy  the  following 
recurrence  relations: 


Vi'^) 


2n 

-0„(z) 


(I.l) 


C„.iU)  -  c„^,(z)  .  2C^(z) 


(1.2) 


\  (  1  ) 

where  C  denotes  J,  Y  and  H  =  H  .  Equation  (I.l)  indicates  that  any 
two  known  consecutuve  menbers  of  a  sequence  {Cp(z)}^  can  be  used  to 
compute  the  entire  sequence.  Moreover,  the  corresponding  sequence  of 
derivatives  {C^(2)}^  can  be  subsequently  computed  via  (1.2).  Thus  it  is 
sufficient  to  consider  how  to  evaluate  the  sequences  |J^(z)}^  and 
{Hn(2)lo  only. 

From  experience  it  is  learned  that  the  calculation  of  the  sequence 
iJn(z)lo  is  stable  numerically  only  if  the  recurrence  relations  are 
applied  downward,  and  unstable  otherwise.  In  addition,  only  a  finite 
sequence  can  be  generated,  i.e.,  N  <  «,  because  of  the  limited  dynamic 
range  of  a  digital  computer. 


One  can  verify  that  the  sequence  satisfies  (I.l)  and 

^  N+1 

(1.2),  where  B  is  an  arbitrary  scalar.  Next  letting  {Jp(z)}Q  be  such 
a  sequence  that 


=  0  .  (1-3) 

and 

J^{z)  =  a  .  (1.4) 

where  a  is  some  vanishingly  small  number.  The  remainders  of  the 
sequence  are  calculated  using  the  recurrence  relation  (I.l)  downward. 

By  setting  =  {bJ^(z)}q,  one  finds  J^(z)  =  J^(2)/B,  n=0,l . N, 

and  thus  the  sequence  {^^(z)}^  will  be  determined  as  soon  as  B  is  known. 

A  A 

Noting  3  =  ’  n=0,l,...,N,  it  is  clear  that  B=Jjj(2)/J^(z) 

is  most  easily  determined  since  Jq(z)  can  be  simply  computed  as  follows: 


1.  When  |z|  <8, 

1  2  k 
-  {-Jz  ) 

^  (k!)^ 

k=n 


(1.5) 


where  the  series  converges  rapidly. 
2.  When  |z|  >8, 

^  I  \  -  f  jzcosB^ 

=  it  /  e*'  d9 

0 


where  the  integral  is  evaluated  by  Simpson's  quadrature. 


Once  the  sequence  is  generated,  one  can  immediately  move 

onto  the  calculation  of  the  sequence  which  is  described  below: 

1.  When  |z|  <8,  Lommel's  formula  [27,  pp.  1431  can  be  employed  to 

evaluate  Y  (z): 

0 


2  z  •  (z/2)  J|<(2) 

^0^^^  ^  IT  (lo9  2  +  ^ 

k=l 


(1.7) 


..11  1 

where  y  =  lim  [1  +  7  +  7  +  ...  +  m  -  Inm]  =  .5772157  is  known  as  Euler's 
m>« 

constant.  Next  it  follows  from  the  Wronskian, 

W{Jn(z).Y^(z)}  =  J^^j(z)Y^(z)  -  J^{z)Y^^j(z)  =  fj,  n=0,l,2,...N-l, 

(1.8) 

N 

that  the  sequence  (Y^(z)}q  is  readily  determined.  Finally,  making  use 

of  the  relation  H„(z)  *  J  (z)  +  jY  (z),  the  sequence  {H  (z)}^  is 
n  n  n  n  '  o 

obtained  quite  simply. 

2.  When  \z\  >  8,  the  sequence  {l^f,(z)}Q  is  directly  computed  by 
first  evaluating  H^(z)  via  Hankel's  asymptotic  expansions  [261: 


[P(z)+jQ(z)l  e 


j(z-iT/4) 


(1.9) 


where 


2  2  2  2 
1.3. 5. 7 


l-3*5*7‘9 


then  employing  the  Wronskian, 

It  must  be  mentioned  that,  as  lmz+»,  evaluating  Jo(z)  and  Ho(z) 
according  to  (1.6)  and  (1.9),  respectively,  will  encounter  severe  num¬ 
erical  problems  commonly  known  as  computer  floating-point  overflow  and 

N 

underflow.  Such  will  prevent  one  to  calculate  the  sequences  {Jp(z)} 

N 

^"<1  (H„(z)L  accurately.  A  simple  but  effective  way  to  handle  these 

problems  is  through  the  proper  scaling  of  J^(z)  and  Hq(z),  and  hence  the 

resulting  sequences.  This  means  that  0^(2)  and  H^(z)  are  being  modified 

through  multiplicative  factors  e^  and  e^,  respectively,  where  s  is  some 

appropriate  positive  real  number.  Letting  {^^(z)}^  and  (H  (z)}^  be  the 

N  N 

scaled  sequences  corresponding  to  {^^(z)}^  =  {H^(z)}q,  one  can  write 


where 


T  /  \  f  \  -  T  jzcos9-s 

Jo(z)  =  e  J  (z)  =  If  /  e''  d9 

0 

^n(z)  are  computed  sequentially  according  to  the  Wronskian, 


(1.13) 


(1.14) 


(1.15) 


W{Jn(z),  H^(2)}  =  Jn+i(z)H„(z)  -  n=0.1 . N-l, 


H^(z)  =  e\(z)  [P(2)  +  jO(2)]  ^  ^ 


(1.16) 


(1.17) 


It  is  important  to  observe  that  the  function  W(n,h)  in  (4.83)  is 

♦■c 

independent  of  the  scaling  factors  e~  .  This  is  obvious  if  one  realizes 
that  (from  Equations  (4.88)-(4.92)) 


Hn(Xb)  Hn(Xb) 
Hn(Xb)  Hn(Xb) 


p  J^(ub)H^(ya)-J^(ua)H^(iib)  J^(iib)H^(ija)-J^(pa)H^(iib) 


Jn(yb)Hn(ya)-Jn(ua)H^(ub)  J^(wb)H^(wa)-vJ^(wa)H^(ub) 


1  f 


J  (ub)H_(ua)-J  (ua)H.(ub)  J  (ub)H.(ua)-J  (iia)H  (ub) 


(1.18) 


J_  ( ub )  H  ( ua )  - J  ( uB )  H  ( ub ) 
nnnn  nnnn 


I 


APPENDIX  J 


ASYMPTOTIC  FORM  OF  0(v) 

In  this  appendix,  an  asymptotic  form  of  0(v),  defined  by  (5.50), 
based  on  the  Debye  approximations  is  derived.  From  (4.110)  and  (4.111), 

( M  (^ ) ' 

the  Deybe  approximations  of  Hy2'(5)  and  can  be  written  as: 


1  2  6 


-  V  cos"^ v/5)-it/4) 


(C2-v2)1/4 


hJ,2^'(c)  = hJ^2^c)  ~  ±  j  , 


where  |5|>|v|,  |5-v|>I I .  and  5  large.  It  follows  that 


Hy  (a)  (3)  2 


fjF(v) 


1/4  1/4 


u(2),  , 

Hy  (a)  (3) 


I  (a2-v2)  (32-v2) 


(.1.1) 


(J.2) 


(.1.3) 


F(v) 

-  +  v(cos  (v/a)  -  cos  (v/S)) 

Also, 

y  '  '  ' 

(1)  ~  ^ 
K  iy) 

Y 

Employing  (J.3)  and 

(J.5)  in  Q(v)  yields 

Recalling  that  a  *  fE^  k2a,  6  =  /F^  k2b,  and  f  «  k2b,  one  can  deduce 
from  (J.6)  that 

L-  JF(v)i 


j  c^k2-(v/a)^|  j’’’ '^2-('>/a)^  j  e'"^ 


jF(v)  JF(v) 


where  it  is  assumed  that  a>>(b-a)  such  that  (v/b)=!{v/a) .  It  can  he 
readily  verified  that  (J.7)  leads  to 
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and 


I 


II 


C^(na)  =  C^(ub)  -  utC^(ub) 

Since  Cv  satisfies  the  Bessel  differential  equation,  one  can  write 

Substituting  (K.5)  into  (K.4)  produces 

c|^(..a)  =  c'jM  [1  +  +  ut  [1  -  (j'l  C^(pb)  . 


(K.A) 


(K.5) 


(K.6) 


It  may  be  pointed  out  that  only  (K.3),  but  not  (K,5),  is  needed  in 
approximating  Zv(h).  This  is  clear  when  one  observes  that 


Z„(M  -  - 


or 


i/'^p 

ut 


[H[^^^(ub)-utH(,^^’(ub)]H(^!^)'(ub)  - 
-  [nJ,^^ub)-iitH[,^^'(Mb)]H(,l)'(pb) 
[Hj,l>(„b).utH(^l^'(ub)]H(^=?)(,ib)  - 
-  [H5,^^ub)-utH5,2)'(ub)]HU)(ub) 

(K.7) 

(K.8) 
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One  can  easily  derive  from  (L.4)  the  following  non-linear  differential 
equation  which  governs  the  roots  of  P(a)  =  0: 


do  1 

dy  ~  l-y’-a  ’  ('-.5) 


It  is  readily  observed  that  y=0  if  only  if  Wj(o)  =  0.  Since  w^ 
vanishes  at 


f  P  ~ 


(L.6) 


where  t^'s  are  the  zeros  of  A^(-t)  (Airy  function),  one  sees  that 


o(0)  =  trt  ,  p  -  1,2,3,... 


(L.7) 


It  is  then  natural  to  denote  the  solution  of  (L.5)  corresponding  to  the 
intitial  condition  o{0)  *  tp  by  op(y),  p  =  1,2,3. 

(n)  _  ^ 


'  •  •  •  • 


Next  letting  o 


.(1), 


<Jp(y)»  Of*®  formally  Taylor  expand 


<Jp  '(y)  about  y=0  as  follows: 


(y) 


y 

o 

n=0 


,(n+l) 

!e _ 


(0) 


n! 


n 

y 


(L.8) 


For  notational  simplicity  the  symbol  o^^^  will  mean  Op^^  in  the 
following  computations,  unless  specified  otherwise.  It  is  obvious 
from  (L.5)  that 


7.(1) 


(2^  2y<j  +  f-o 


=  [2yo  +  y2a^^^][a^^^]  ,  (L.IO) 


a^^^(0)=0  . 


(L.ll) 


=  [2<j+2ya^^U2y(j^^Uy^a^^^][o^^^]^+  [2ya+y^a^^^  ][2<j^^^o^^^] 


=  2a[a(l)^  {4y[a(l)]%4yaa(l^a(2)}V[a('^]%('^  . 


(L.12) 


=  2o(0)[n^  =  2t 


(L.13) 


(L.14) 


cr^^^O)  =  6  . 


(L.15) 
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*  22[.(l>w'’'*4,<"(<>'^’r*8j2>o<3'*4[o(l>]%<S> 

.  4,<.'2>c'3>.4j‘>a<‘'>|  *  2yl6<.<"[,''>A3[o<»l'<.''’l 

=  48[a^^^]  a^^^+12o[a^^^ J  +12ao^^^o^^^+y {60a^^^ ] 


^0)  =  24t? 


.  12o<'>»<2)„(3)*l2oa(2>,<^!*4[o("j%('') 

.  4J'>„<’W‘>o'5)}  *  y|6[„<»^(^Ul2[<,(2>f 

15S,<"[o<2'l%92[c<'>)%<3),is,a»)„(«)*48oa<2>,<3Uy|7; 
*  232»<l',<2><,<3Ul2,l,(3)]%16j2>„(‘')*4j'),'5)*,2[, 

,(6) 


(0)  =  280t„ 
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